超短脉冲激光在二维非均匀介质内的瞬态辐射传输

随着超短脉冲激光的快速发展,吸收散射性介质内的瞬态辐射传输引起了广泛的关注.本文采用最小二乘有限元法模拟了超短脉冲激光局部入射条件下,具有高散射核的二维非均匀介质内的瞬态辐射传输.研究了不同边界位置上反射和透射信号随时间的变化情况.结果表明,对于具有高散射核的非均匀介质,能够揭示散射核位置的双峰现象可能在早期的反射信号中发生,早期的反射信号比后期的、经过介质衰减后信号更重要.因此,在利用短脉冲激光进行光谱分析和成像等技术中,人们应该重视早期反射信号的测量.     

NMR 中非线性梯度场下扩散行为的理论描述和蒙特卡罗模拟

利用传播子方法研究了在一般非线性场梯度下NMR信号的扩散衰减. 在自由扩散和平板间的限制扩散情况下获得了扩散衰减因子的理论表达式. 该表达式适用范围宽，且具有较简单的数学形式和明确的物理意义. 文中还将理论预测与蒙特卡罗模拟结果进行了比较. 结果表明：文中所采用的理论方法适合于表述自由扩散和短脉冲近似下的受限扩散；蒙特卡罗模拟提供了一种定性研究MRI和NMR中非均匀场梯度扩散衰减的方法.     

用于同时多层MRI的选择性射频脉冲的优化设计

近年来,为提高磁共振成像（MRI）信号信噪比（SNR）、缩短成像时间,同时多层成像技术受到了极大的关注.为了实现同时多层的选择性激发,现有的多层成像序列大多使用组合射频（RF）脉冲,该脉冲可包含多个独立的幅值相同相位不同的简单脉冲,由于其采用简单的线性叠加方法,该类脉冲射频功率随脉冲数量呈现平方增长,因而应用受限.针对这一问题,基于自旋动力学和优化控制原理,本文提出了一种针对同时多层MRI的选择性射频脉冲的数值优化方法,该方法充分运用射频脉冲的调控机制,获得优化脉冲,并配合层选梯度,可实现任意层厚、层间距、层数的同时高效选择性激发.最后,通过数字模体的同时多层模拟成像实验验证了优化脉冲的有效性.     

Bloch方程的精细时程积分及其在射频脉冲设计中的应用

射频脉冲的频率选择性直接影响磁共振成像的质量，而射频脉冲的优化设计又归结为对Bloch方程的求解.尽管在某些情况下Bloch方程存在解析解，但由于其缺乏通用性而且形式上过于复杂而难于得到实用.本文提出一种Bloch方程的精细时程积分算法，并结合全局优化算法给出一个完整的射频脉冲设计方案.精细积分算法具有高效、高精度的特点，对于射频脉冲的设计很有裨益.数值算例表明，设计所得的射频脉冲具有较好的频率选择性. 关键词： 磁共振成像 射频脉冲 Bloch方程 精细时程积分     

基于非均匀螺旋线数据和布雷格曼迭代的快速磁共振成像方法

数据采集时间长是制约磁共振成像技术发展的重要瓶颈.为了解决这一问题, 本文基于压缩感知成像理论, 提出了一种结合非均匀螺旋线磁共振数据采集序列和布雷格曼迭代重建的快速磁共振成像方法, 通过欠采样缩短数据采集时间.欠采样引起混迭伪影则通过非均匀螺旋线欠采样特性和布雷格曼迭代重建去除.水模磁共振成像实验和在体磁共振成像实验结果表明: 欠采样情况下, 所提出的方法能有效去除欠采样导致的混迭伪影, 获得的图像结构信息完整的成像结果, 在缩短采样时间的同时, 具有较高的准确度. 关键词： 磁共振成像 非均匀螺旋线 全变分 布雷格曼迭代     

超短THz脉冲在随机散射介质中的传播

用Monte-Carlo方法模拟了超短THz脉冲在随机散射介质中的传播.根据Mie理论计算出随机散射介质的散射系数和各向异性因子，研究了随机散射颗粒及介质厚度的大小对透射脉冲的影响.结果表明在Mie散射范围内，在相同的浓度下，颗粒尺寸越小，散射介质越厚，THz散射越严重，对透射脉冲的影响越大.散射同时降低了THz脉冲在随机散射介质中的成像分辨率. 关键词： 超短THz脉冲 随机散射介质 Mie理论 Monte-Carlo方法     

磁逆多光子非线性Compton散射的能量转换

应用电子与光子非弹性碰撞模型，分析、计算了磁逆多光子非线性COmpton散射中电子与光子间的能量转换。结果表明：磁逆多光子非线性Compton散射前，做Landau运动的电子不是热电子，而是处于激发态的电子，该散射产生高γ射线的几率要比磁逆单光子Compton散射时的几率小；当磁逆多光子非线性Compton散射中的电子被光场俘获时，电子可以从激光场中获得巨大的加速能量，但比多光子非线性Compton散射时获得的能量低。     

低场永磁体磁共振射频场映像

射频场映像是通过一定算法对磁共振射频线圈的发射场进行重建的方法.高场下的射频场经过生物组织时会发生明显变化,在其基础上可以反演生物组织体内电特性,进而对癌症等疾病进行早期诊断,是对生物组织的磁共振结构成像的有力补充.目前为止,射频场映像和电特性研究都以高场鸟笼线圈为主,对低场下的相控阵研究较少.本文主要研究了低场永磁体磁共振射频场的均匀度.有限元仿真和实验验证了在17.8 MHz激励下,射频场在空载和负载下均匀度均发生较大变化.射频场均匀度在负载下的改变在一定程度上可以反映负载生物组织的电特性,对磁共振电特性实用化研究提供了一定的参考价值.     

超宽频带THz脉冲在随机散射介质中传播的理论研究

研究了超宽频带THz脉冲入射到散射介质中时其透射脉冲的时空特性.根据Mie理论计算出随机散射介质的散射系数和各向异性因子,用时间分辨Monte Carlo方法模拟了超宽频带的THz脉冲在随机散射介质中的传播,研究了在不同散射颗粒半径,不同频宽的THz入射脉冲对透射脉冲的影响和散射对成像分辨率的影响.结果表明:散射会降低THz脉冲在随机散射介质中的成像分辨率,散射颗粒越小,散射介质越厚,其成像分辨率越差.     

一种用于高场MRI的多源射频发射机

介绍了高场磁共振成像（MRI）多源发射技术的原理,提出了一种用于高场MRI系统的多源射频信号发射机.它能并行输出多路频率、相位、幅度,可快速独立调节的射频脉冲信号.该射频发射机的实现基于单片现场可编程门阵列（FPGA）和多通道数模转换器（DAC）芯片,FPGA读取预存于双端口随机存取存储器（RAM）中的射频信号参数,并利用读取的参数分别实现每路信号的直接数字频率合成（DDS）和信号调制等核心功能,获得多路数字射频信号;FPGA输出的数字信号经过高性能DAC转化为模拟信号,即所需要的射频信号.该射频发射机在设计中大量采用软件无线电技术,即利用Xilinx提供的IP核实现DDS和信号调制等主要功能,具有集成度高、体积小、灵活度高的优点,同时,该设计可以大大缩短开发时间,有效降低实现的难度和成本,为高场MRI谱仪的多源射频发射机的设计研制提供了一种低成本、高效、高性价比的方案.     

针对限制幅值磁共振脉冲的优化设计

在磁共振脉冲优化领域,优化脉冲普遍存在幅值过大的问题,这极大地限制了优化脉冲的使用范围。为了限制优化脉冲的幅值,扩大其应用范围,提出了一种基于L-BFGS-B数值算法的脉冲优化设计方法。首先基于Liouville-von Neuman方程,使用最优控制思想构建优化模型;然后使用L-BFGS-B算法,在限制幅值的条件下对优化模型进行数值迭代求解;最后以脉冲的激发效率以及激发轮廓的均匀性作为衡量优化脉冲优劣的标准对该方法进行仿真和实验验证。结果表明,采用该方法获得的优化脉冲在幅值被限制的前提下,仍能获得较传统磁共振脉冲更好的共振激发效果,进而增强信号的灵敏度,提高图像的质量。     

Practical pulse synthesis via the discrete inverse scattering transform

This paper provides the practical details required to use the inverse scattering (IST) approach to design selective RF-pulses. As in the Shinnar-Le Roux (SLR) approach, we use a hard pulse approximation to actually design the pulse. Unlike SLR, the pulse is designed using the full inverse scattering data (the reflection coefficient and the bound states) rather than the flip angle profile. We explain how to approximate the reflection coefficient to obtain a pulse with a prescribed rephasing time. In contrast to the SLR approach, we retain direct control on the phase of the magnetization profile throughout the design process. We give explicit recursive algorithms for computing the hard pulse from the inverse scattering data. These algorithms are quite different from the SLR recursion, being essentially discretizations of the Marchenko equations. We call our approach the discrete inverse scattering transform or DIST. Overall, it is as fast as the SLR approach. When bound states are present, we use both the left and right Marchenko equations to improve the numerical stability of the algorithm. We compute a variety of examples and consider the effect of amplitude errors on the magnetization profile.     

Semi-adiabatic Shinnar–Le Roux pulses and their application to diffusion tensor imaging of humans at 7T

The adiabatic Shinnar–Le Roux (SLR) algorithm for radiofrequency (RF) pulse design enables systematic control of pulse parameters such as bandwidth, RF energy distribution and duration. Some applications, such as diffusion-weighted imaging (DWI) at high magnetic fields, would benefit from RF pulses that can provide greater B₁ insensitivity while adhering to echo time and specific absorption rate (SAR) limits. In this study, the adiabatic SLR algorithm was employed to generate 6-ms and 4-ms 180° semi-adiabatic RF pulses which were used to replace the refocusing pulses in a twice-refocused spin echo (TRSE) diffusion-weighted echo planar imaging (DW-EPI) sequence to create two versions of a twice-refocused adiabatic spin echo (TRASE) sequence. The two versions were designed for different trade-offs between adiabaticity and echo time. Since a pair of identical refocusing pulses is applied, the quadratic phase imposed by the first is unwound by the second, preserving the linear phase created by the excitation pulse. In vivo images of the human brain obtained at 7 Testa (7 T) demonstrate that both versions of the TRASE sequence developed in this study achieve more homogeneous signal in the diffusion-weighted images than the conventional TRSE sequence. Semi-adiabatic SLR pulses offer a more B₁-insensitive solution for diffusion preparation at 7 T, while operating within SAR constraints. This method may be coupled with any EPI readout trajectory and parallel imaging scheme to provide more uniform coverage for diffusion tensor imaging at 7 T and 3 T.     

相位可控的核四极矩共振激励脉冲发生器设计

在核四极矩共振（NQR）领域,射频激励脉冲信号的优劣对NQR响应信号有重要影响.针对常规方法中射频激励脉冲参数不可控的问题,本文基于32位闪存微型控制器STM32和直接数字频率合成（DDS）芯片AD9910设计了一种相位可控激励脉冲发生器.采用STM32控制AD9910产生波形参数（脉冲宽度、脉冲间隔、脉冲个数和共振频率等）可控的射频激励脉冲,利用LabVIEW软件平台设计脉冲参数设置界面,并建立计算机与微控制器通信,实现波形参数的精确优化控制.实验结果表明,该方法实现了相位可控的NQR激励脉冲序列,可为后续NQR信号检测提供有效激励源.     

An integrated program for amplitude-modulated RF pulse generation and re-mapping with shaped gradients

Efficient generation of amplitude modulated, frequency selective RF pulses has been demonstrated by the Shinnar-Le Roux (SLR) algorithm. In the present article, we provide an overview of a relatively comprehensive computer program that includes a version of the SLR algorithm and also incorporates an algorithm for re-mapping a selective RF pulse onto a new dwell time with modulated gradients. The re-mapping may be used to reduce SAR, or to shorten the RF pulse time by increasing the gradient and RF strength in regions where the original RF pulse amplitude was low. The program includes additional useful features including a Bloch equations algorithm, and pulse scaling, to enable examination of pulse profiles under a variety of conditions such as RF inhomogeneity and even nuclear relaxation. The program, MATPULSE, was developed with the MATLAB for Windows programming language and makes extensive use of the MATLAB graphical user interface (GUI) features to generate a userfriendly interface. A number of examples are provided to illustrate the capabilities of the MATPULSE program.     

Reducing SAR in parallel excitation using variable-density spirals: a simulation-based study

Parallel excitation using multiple transmit channels has emerged as an effective method to shorten multidimensional spatially selective radiofrequency (RF) pulses, which have a number of important applications, including B1 field inhomogeneity correction in high-field MRI. The specific absorption rate (SAR) is a primary concern in high-field MRI, where wavelength effects can lead to local peaks in SAR. In parallel excitation, the subjects are exposed to RF pulses from multiple coils, which makes the SAR problem more complex to analyze, yet potentially enables greater freedom in designing RF pulses with lower SAR. Parallel-excitation techniques typically employ either Cartesian or constant-density (CD) spiral trajectories. In this article, variable-density (VD) spiral trajectories are explored as a means for SAR reduction in parallel-excitation pulse design. Numerical simulations were conducted to study the effects of CD and VD spirals on parallel excitation. Specifically, the electromagnetic fields of a four-channel transmit head coil with a three-dimensional head model at 4.7 T were simulated using a finite-difference time domain method. The parallel RF pulses were designed and the resulting excitation patterns were generated using a Bloch simulator. The SAR distributions due to CD and VD spirals were evaluated quantitatively. The simulation results show that, for the same pulse duration, parallel excitation with VD spirals can achieve a lower SAR compared to CD spirals for parallel excitation. VD spirals also resulted in reduced artifact power in the excitation patterns. This gain came with slight, but noticeable, degrading of the spatial resolution of the resulting excitation patterns.     

Improving excitation and inversion accuracy by optimized RF pulse using genetic algorithm

In this study, a Genetic Algorithm (GA) is introduced to optimize the multidimensional spatial selective RF pulse to reduce the passband and stopband errors of excitation profile while limiting the transition width. This method is also used to diminish the nonlinearity effect of the Bloch equation for large tip angle excitation pulse design. The RF pulse is first designed by the k-space method and then coded into float strings to form an initial population. GA operators are then applied to this population to perform evolution, which is an optimization process. In this process, an evaluation function defined as the sum of the reciprocal of passband and stopband errors is used to assess the fitness value of each individual, so as to find the best individual in current generation. It is possible to optimize the RF pulse after a number of iterations. Simulation results of the Bloch equation show that in a 90 degrees excitation pulse design, compared with the k-space method, a GA-optimized RF pulse can reduce the passband and stopband error by 12% and 3%, respectively, while maintaining the transition width within 2 cm (about 12% of the whole 32 cm FOV). In a 180 degrees inversion pulse design, the passband error can be reduced by 43%, while the transition is also kept at 2 cm in a whole 32 cm FOV.     

Design and optimization for variable rate selective excitation using an analytic RF scaling function

At higher B(0) fields, specific absorption rate (SAR) deposition increases. Due to maximum SAR limitation, slice coverage decreases and/or scan time increases. Conventional selective RF pulses are played out in conjunction with a time independent field gradient. Variable rate selective excitation (VERSE) is a technique that modifies the original RF and gradient waveforms such that slice profile is unchanged. The drawback is that the slice profile for off-resonance spins is distorted. A new VERSE algorithm based on modeling the scaled waveforms as a Fermi function is introduced. It ensures that system related constraints of maximum gradient amplitude and slew rate are not exceeded. The algorithm can be used to preserve the original RF pulse duration while minimizing SAR and peak b1 or to minimize the RF pulse duration. The design is general and can be applied to any symmetrical or asymmetrical RF waveform. The algorithm is demonstrated by using it to (a) minimize the SAR of a linear phase RF pulse, (b) minimize SAR of a hyperbolic secant RF pulse, and (c) minimize the duration of a linear phase RF pulse. Images with a T1-FLAIR (T1 FLuid Attenuated Inversion Recovery) sequence using a conventional and VERSE adiabatic inversion RF pulse are presented. Comparison of images and scan parameters for different anatomies and coils shows increased scan coverage and decreased SAR with the VERSE inversion RF pulse, while image quality is preserved.     

2H T2rho relaxation dynamics and double-quantum filtered NMR studies

In this study 2H T2rho DQF NMR spectra of water in MCM-41 were measured. The T2rho double-quantum filtered (DQF) NMR signal is generated by applying a radio frequency (RF) field for various durations and then observed after a monitor RF pulse. It was found that the transfer between different quantum coherences by the couplings during long-duration RF fields (i.e., soft pulses) and that residual quadrupolar interaction dominates the signal decay. Knowledge of coherence transfer during long-RF pulses has special significance for the development of sophisticated multi-quantum NMR experiments especially multi-quantum MRI applications.