一种用于多叶准直器适形放疗的剂量算法研究

描述了一种用于多叶准直器适形放疗的剂量计算方法。 在不同大小的方野和非规则野照射情况下, 将此方法得到的剂量分布与微分卷积法计算得到的剂量分布进行了精度验证。 结果表明, 这种算法具有较高的剂量精度和较快的计算速度, 因此可在适形放疗中用作MLC适形野的自动生成。 A dose calculation algorithm for MLC based conformal radiotherapy is described in this paper. The algorithm is formulated by the coordinate of MLC leaves. Verification on the algorithm is made by comparing the dose distributions generated by this algorithm with that generated by a Differential Convolution Superposition algorithm for various regular and irregular fields. The results demonstrate that the present algorithm has suitable accuracy and high computational efficiency, thus it could be useful for the treatment planning process in MLC based conformal radiotherapy, where the workload for interactively or automatically designing the shapes of MLC is heavy.     

Imaging photothermal microscopy for absorption measurements of optical coatings

For absorption measurement of large-aperture optical coatings, a novel method of imaging photothermal microscopy based on image lock-in technique is presented. Detailed theoretical analysis and numerical calculation are made based on the image photothermal technique. The feasibility of this imaging method is proved through the coincidence between the theoretical results of single spot method and multi-channel method. The measuring speed of this imaging method can be increased hundreds of times compared with that of the raster scanning. This technique can expand the applications of photothermal technique.     

Specific absorbed fractions of electrons and photons forRad-HUMAN phantom using Monte Carlo method

The specific absorbed fractions (SAF) for self- and cross-irradiation are effective tools for the internal dose estimation of inhalation and ingestion intakes of radionuclides. A set of SAFs of photons and electrons were calculated using the Rad-HUMAN phantom, which is a computational voxel phantom of a Chinese adult female that was created using the color photographic image of the Chinese Visible Human (CVH) data set by the FDS Team. The model can represent most Chinese adult female anatomical characteristics and can be taken as an individual phantom to investigate the difference of internal dose with Caucasians. In this study, the emission of mono-energetic photons and electrons of 10 keV to 4 MeV energy were calculated using the Monte Carlo particle transport calculation code MCNP. Results were compared with the values from ICRP reference and ORNL models. The results showed that SAF from the Rad-HUMAN have similar trends but are larger than those from the other two models. The differences were due to the racial and anatomical differences in organ mass and inter-organ distance. The SAFs based on the Rad-HUMAN phantom provide an accurate and reliable data for internal radiation dose calculations for Chinese females.     

弧形调强放射治疗对剂量计算方法的要求

在弧形调强放射治疗的治疗计划设计中, 由于包含有很多照射方向, 调强最优化的射束元矩阵计算需要很大的计算量和存储量, 为提高计算效率常使用简化剂量计算模型计算射束元矩阵, 因此有必要研究简化模型对治疗计划质量产生影响。 对一个模拟例子和一个临床实例, 使用没考虑散射效应的原射线剂量计算模型计算射束元矩阵, 由此进行最优化计算。 在得到最优化强度分布后, 通过比较原射线剂量计算模型和微分卷积剂量计算模型得到的剂量分布, 研究了不同射束数目条件下, 使用简化剂量计算模型计算射束元剂量矩阵对最终的剂量分布质量的影响。 结果表明, 在射线束很多的情况下(对应弧形调强照射）, 用简化的剂量计算模型, 即不考虑散射来计算射束元剂量矩阵, 会导致靶区剂量分布的质量大大低于预期的剂量分布质量, 因此, 弧形调强放射治疗的最优化计算中, 有效考虑散射的影响是必要的。 In the treatment planning for arc intensity modulated radiation therapy, because many irradiation directions are involved, the computing time and storage space needed for calculating beamlet dose matrices in optimization is quite heavy. In order to improve the computation efficiency, the simplified dose calculation is often used for the calculation of the dose matrices. Thus, it is deserved to study how this simplification could influence the quality of the treatment plan. In this paper, a simulation and a clinical case are adopted. Using the primary dose calculation model without taking into account the scattering effect to generate the dose matrices of beamlets, the optimization for beam intensity profile are firstly carried out. Then, based on the obtained intensity profile, the dose distributions are recalculated by using the primary dose calculation model and the differential convolution superposition dose calculation model which is more accurate but more time consuming. By comparing dose distributions obtained by this two models, the influence of using simplified model for dose matrix calculation on beam profile optimization is studied. The results demonstrate that when the beam number is large(corresponding to the arc modulated radiation), using the simplified model for the calculation of dose matrix of beamlets will reduce the quality of dose distribution greatly comparing with the expected dose distribution quality. Thus it is very necessary to correctly take into account the scattering effect in beam profile optimization for the arc intensity modulated radiation therapy.     

ANGULAR DISTRIBUTION OF THE SOUND FIELD IN SHALLOW WATER

The field of a harmonic point source in a layered medium has a certain angular distribution.Although the angular distribution of the sound field can be calculated for a given receiving point bythe ray theory,the effect of the frequency on the angular distribution function and the sound field inshadow zones can not be calculated.In this paper,the angular distribution of the sound field in layeredmedia is discussed based on the normal-mode theory,and the general expression of the angular dis-tribution function and the attenuation of the sound field with distance for directional reception areobtained.The results given in the paper can be used for the calculation of the sound field in shadowzones and they represent the influence of the frequency on the sound field.When the underwater acous-tic channel and the performance of signal detection in it are investigated,the optimum receiving direc-tion can be chosen by utilizing the obtained results.For a given medium model the numerical cal-culations are conducted     

The application of arbitrary incidence laser beams heat treatment temperature field calculation formulas

Based on the calculation formulas of heat treatment temperature field for an arbitrary incident laser intensity distribution, the transformation intensity distribution of CO2 laser beam passing an integrating mirror is studied theoretically and experimentally. The derived formulas are applied in laser heat treatment research which is transformed by optical system, and the theoretical calculation results are compared with experimental results. It is shown that the formulas can be used to calculate the laser heat treatment temperature field accurately, and the calculation speed is obviously faster than the numerical calculation methods with the same precision. The calculation software can be used to select proper experiment parameters.     

The application of arbitrary incidence laser beams heat treatment temperature field calculation formulas

Based on the calculation formulas of heat treatment temperature field for an arbitrary incident laser intensity distribution, the transformation intensity distribution of CO2 laser beam passing an integrating mirror is studied theoretically and experimentally. The derived formulas are applied in laser heat treatment research which is transformed by optical system, and the theoretical calculation results are compared with experimental results. It is shown that the formulas can be used to calculate the laser heat treatment temperature field accurately, and the calculation speed is obviously faster than the numerical calculation methods with the same precision. The calculation software can be used to select proper experiment parameters.     

高Q2电子原子核单举准弹性散射截面的计算方法

电子与原子核碰撞实验是通过中高能探针探测原子核结构的方法。本文提供了一个高Q2单举准弹性电子原子核散射截面的计算方法,此方法是基于核子-核子短程关联的经验公式与弱束缚近似下的氘核散射截面模型。在弱束缚近似下,氘核可以看成是由近似自由质子与中子组成,质子与中子的短程关联可以忽略,氘核结构函数可以写成质子与中子结构函数线性组合,从而可以得到氘核的散射截面。根据氘核散射截面以及短程关联的经验公式,可以得到考虑短程关联的原子核A > 2的散射截面。我们将计算得到的散射截面与已有的实验和及Bosted拟合方法的结果比较,发现本文的方法在大x和Q2 > 2 GeV2区域得到的结果与实验结果符合得较好并且对于一些重原子核,尤其是4He核,明显优于Bosted方法的结果。The electron nucleus collision experiments are approaches measuring the structure of nuclei by using intermediate and high energy probe. This paper shows a calculation method of inclusive electron nucleus quasielastic scattering cross section at high Q2 which based on a empirical formula of Nucleon-Nucleon Short Range Correlation (NN-SRC) and a model of electron deuteron quasi-elastic cross section in Weak-Binding Approximation (WBA). In WBA, the deuteron can be regarded as the combination of quasi-free proton and neutron and the short range correlation between them can be ignored. Therefore the structure function of deuteron can be written as the linear combination of that of proton and neutron, then one can get the cross section of deuteron. According to the cross section of deuteron and the empirical formula, one can obtain the cross section of nuclei A > 2 which considers NN-SRC effect. We compare our calculation results with existing experiments and the results calculated by Bosted' fit method, then find that our results match the experiments at high x and Q2 > 2 GeV2 and better significantly than the Bosted's results for some heavier nucleus, especially 4He.     

Two methods to deal with an inhomogenous Rabi frequency in microwave transition

We analyse the influence of an inhomogenous microwave field on the coherence of atom ensembles. Two methods are proposed to suppress the dephasing generated by the inhomogenous Rabi frequency. One of them is realized by using a spin echo, and the other one is based on the identical spin rotation effect. The calculation results show that the contrast of a signal acquired in experiment can be improved by using the two methods. Their advantages and drawbacks are discussed. We hope they can be used to improve the contrast of experimental signals in situations where microwave fields are very inhomogenous. Finally, we discuss the case of a continuous working microwave field and show that the dipole force raised with the inhomogeneity can be eased by spin flip.     

Effect of Media on the Electric Field of a Rhombic Nanostructure Array

Electric field distribution is an important parameter for nanostructure arrays in nanobiosensing appfications. It can influence the sensitivity and the resolution of nanobiosensors. We focus on the effect of media on the electric field distribution of a rhombic silver nanostructure array. The finite-difference time-domain algorithm- based numerical calculation method is used to monitor the electric field distribution of the silver nanostructures when the refractive index of the medium around the nanostructure array is changed. The calculated results show that tuning the refractive index of the medium around silver can have a considerable influence on the electric field distribution in the reflection and transmission directions. This effect can be used to increase the extinction efficiency and to improve the resolution of the spectra for nanobiosensing.     

基于GPU加速的质子调强放疗鲁棒优化器

开发一种基于图形处理器（GPU）加速的质子调强放疗鲁棒优化器,用于减小质子束射程不确定性和靶区定位偏差对质子放疗的影响。建立的鲁棒优化模型使用的目标函数包括9种边界剂量目标,分别是:无偏差情况、2种射程偏差（偏长与偏短）、6种摆位不确定性（前后、侧向、上下入射方向各2种正负偏差）。首先靶区和危及器官的剂量贡献矩阵使用笔形束算法计算得到,然后使用共轭梯度法优化目标函数让其满足约束条件,这两部分均采用GPU加速。头颈部、肺部和前列腺三个临床病例被用来检测本优化器的性能表现。与传统基于计划靶区（PTV）的质子调强放疗计划相比,鲁棒优化器能够优化出对射程不确定性和摆位误差更加不敏感的治疗计划,让靶区实现了高剂量均匀性的同时危及器官（OARs）也得到了更好的保护。经过100次迭代,三个病例的优化时间均在10 s左右。该结果证明了基于GPU加速的质子调强放疗鲁棒优化器能够在短时间内设计出高鲁棒性的质子治疗计划,从而提高质子放射治疗的可靠性。This paper describes the development of a fast robust optimization tool that takes advantage of the GPU technologies. The objective function of the robust optimization model considered nine boundary dose distributions--two for ±range uncertainties, six for ±set-up uncertainties along anteroposterior (A-P), lateral (R-L) and superior{inferior (S-I) directions, and one for nominal situation. The nine boundary influence matrices were calculated using an in-house dose engine for proton pencil beams of a finite size, while the conjugate gradient method was applied to minimize the objective function. The GPU platform was adopted to accelerate both the proton dose calculation algorithm and the conjugate gradient method. Three clinical cases-one head and neck cancer case, one lung cancer case and one prostate cancer case-were investigated to demonstrate the clinical significance of the proposed robust optimizer. Compared with conventional planning target volume (PTV) based IMPT plans, the proposed method was found to be conducive in designing robust treatment plans that were less sensitive to range and setup uncertainties. The three cases showed that targets could achieve high dose uniformity while organs at risks (OARs) were under better protection against setup and range errors. The run times for the three cases were around 10 s for 100 iterations. The GPU-based fast robust optimizer developed in this study can serve to improve the reliability of traditional proton treatment planning by achieving a high level of robustness in a much shorter time.     

Restricted Boltzmann machine: Recent advances and mean-field theory

This review deals with restricted Boltzmann machine(RBM) under the light of statistical physics.The RBM is a classical family of machine learning(ML) models which played a central role in the development of deep learning.Viewing it as a spin glass model and exhibiting various links with other models of statistical physics,we gather recent results dealing with mean-field theory in this context.First the functioning of the RBM can be analyzed via the phase diagrams obtained for various statistical ensembles of RBM,leading in particular to identify a compositional phase where a small number of features or modes are combined to form complex patterns.Then we discuss recent works either able to devise mean-field based learning algorithms;either able to reproduce generic aspects of the learning process from some ensemble dynamics equations or/and from linear stability arguments.     

Green's function Monte Carlo method combined with restricted Boltzmann machine approach to the frustrated J1-J2 Heisenberg model

Restricted Boltzmann machine (RBM) has been proposed as a powerful variational ansatz to represent the ground state of a given quantum many-body system. On the other hand, as a shallow neural network, it is found that the RBM is still hardly able to capture the characteristics of systems with large sizes or complicated interactions. In order to find a way out of the dilemma, here, we propose to adopt the Green's function Monte Carlo (GFMC) method for which the RBM is used as a guiding wave function. To demonstrate the implementation and effectiveness of the proposal, we have applied the proposal to study the frustrated J₁-J₂ Heisenberg model on a square lattice, which is considered as a typical model with sign problem for quantum Monte Carlo simulations. The calculation results demonstrate that the GFMC method can significantly further reduce the relative error of the ground-state energy on the basis of the RBM variational results. This encourages to combine the GFMC method with other neural networks like convolutional neural networks for dealing with more models with sign problem in the future.     

重离子治疗计划的生物有效剂量及其验证

介绍了重离子治疗计划中的生物有效剂量，以及基于局部效应模型的相对生物学效应的计算. 报告了用生物测试系统测得的实验结果，并与治疗计划的计算结果做了比较. 靶子内外结果的一致性显示了重离子治疗计划的可靠性及保护紧挨着靶的重要器官的能力.     

Anisotropic effects on the radial breathing mode of silicon nanowires: An ab initio study

The effect of orientation on the frequency of the radial breathing mode (RBM) of silicon nanowires (SiNWs) is investigated by means of the first principles Density Functional Theory approach through the generalized gradient approximation. We compare the RBM frequency of SiNWs orientated in three different directions, [0 0 1], [1 1 1], and [1 1 0]. The RBM is observed by the calculation of the phonon band structure and density of states of the SiNWs through the supercell finite displacement method. Results show that the SiNWs are stable in the three chosen directions since there are no negative frequencies in their phonon band structure and density of states. A clear dependence of the RBM frequency with respect to the growth direction of the nanowires and the phonon confinement was observed as the RBM frequency decreased with an inverse power law in each nanowire direction, with the fitting parameters dependent on the growth direction. These results are important since they could be used as a fingerprint to identify them within different spectroscopy techniques such as Raman.     

Radial breathing vibration of double-walled carbon nanotubes subjected to pressure

A theoretical vibrational analysis of the radial breathing mode (RBM) of double-walled carbon nanotubes (DWCNTs) subjected to pressure is presented based on an elastic continuum model. The results agree with reported experimental results obtained under different conditions. Frequencies of the RBM in DWCNTs subjected to increasing pressure depend strongly on circumferential wave numbers, but weakly on the aspect ratio and axial half-wave numbers. For the inner and outer tubes of DWCNTs, the frequency of the RBM increases obviously as the pressure increases under different conditions. The range of variation is smaller for the inner tube than the outer tube.     

气氧煤油发动机尾焰红外辐射特性研究

提出了一种考虑碳烟颗粒的气氧煤油发动机尾焰红外辐射特性计算方法,首先对气氧煤油发动机纯气相内流场进行计算,然后以喷管喉部作为气体和固体碳烟颗粒的入口边界计算发动机尾焰流场,最后以发动机流场参数分布为基础,采用有限体积法和伪气体理论对发动机尾焰红外辐射特性进行计算。进行了气氧煤油发动机点火实验,并将计算结果与实验结果进行对比分析。结果表明,燃烧室内两个压力测量点的测量与计算误差分别为1.4%和3.4%,燃烧室内计算温度与热力学计算误差为2.16%,证明了燃烧室流场计算模型的准确性。含有碳烟颗粒的尾焰流场计算结果与热像仪测量结果比较吻合,证明了尾焰流场计算方法和模型的准确性。4.3 μm波段尾焰红外成像计算结果与工作在4.3 μm波段的红外热像仪测量结果吻合比较一致,证明了尾焰红外辐射特性计算方法和模型的准确性。     

目标函数设置对放疗逆向计划多目标优化过程的影响

目前放疗逆向计划中常用的目标函数有两种:基于剂量分布的目标函数和基于剂量．体积直方图（DVH）的目标函数。实际系统都基于单目标优化算法进行,不考虑逆向计划的多目标性。在Pareto多目标优化理论的基础上,研究逆向计划的目标函数设置问题,比较基于剂量分布的目标函数和基于DVH的目标函数对多目标优化过程的影响,包括优化时间、收敛性和存在的问题等,为逆向计划过程中多目标优化目标函数的设置提供依据。There are two kinds of objective functions in radiotherapy inverse planning: dose distribution-based and Dose-Volume Histogram （DVH）-based functions. The treatment planning in our days is still a trial and error process because the multi-objective problem is solved by transforming it into a single objective problem using a specific set of weights for each object. This work investigates the problem of objective function setting based on Pareto multi-optimization theory, and compares the effect on multi-objective inverse planning of those two kinds of objective functions including calculation time, converge speed, etc. The basis of objective function setting on inverse planning is discussed.