旌节花科的分类和系统学研究

旌节花科是东亚特有科，仅有旌节花属1属。本文研究了该科植物的生活习性、外部形态、花粉形态和外壁超微结构以及花器官发生，并根据多个DNA区段的序列探讨了科的系统位置和属下关系。在此基础上，综合各方面的证据，对旌节花科进行了全面的分类学修订，完善了该科的分类系统。主要内容包括： 　　1、习性和外部形态学 　　结合室内标本和野外考察，对旌节花科植物的生活习性、叶、花序、花和果实各个性状的变异式样、变异幅度、变异规律及其相关性进行了深入研究，并讨论了它们的分类学价值。结果表明：叶片形状、叶缘齿、花序下的叶片是否宿存以及花序是否具花序梗等性状在种间变异非常大，单个性状对旌节花属属下分类意义不大;而叶片毛被、脉序和结构，花序长度，花的大小和颜色以及果实大小等性状通常较稳定，可作为旌节花属划分种和变种的重要依据。 　　2、花粉形态学和外壁超微结构 　　在透射电镜下研究了旌节花属5个代表种和1个外类群（Crossosoma californicum）的花粉外壁超微结构，并在扫描电镜下研究了外类群Crossosomataceae 1个代表种（Crossosoma californicum）的花粉形态学特征，结合前人的研究，探讨了上述性状的分类学意义。结果表明：旌节花属植物花粉形态的共同特征是：花粉粒球形或近球形，三孔沟;外壁具穴状纹饰，穴较细小;外壁内壁分层明显，容易区分（柳叶旌节花除外），覆盖层光滑，几无附属物，柱状层明显，穿孔较少，不形成半覆盖层。各种植物花粉外壁超微结构的特征在旌节花属两个组中较一致，但其差异可作为探讨属下种的分类和系统学关系的依据。此外，对旌节花科和Crossosomataceae花粉表面纹饰和外壁超微结构的研究结果提示：这两个科与省沽油科和山茶科厚皮香亚科有较近缘的关系，而不支持这两个科的姐妹群关系。 　　3、花器官发生 　　首次对旌节花属常绿组和落叶组共2个代表种的花器官发生进行了研究，结果表明：该科植物花为功能性单性花，雌雄异株，花发育过程长达10个月，其中2-3个月是冬季休眠期。在花发育早期，雄花和雌花的发育没有差异，花性别的分化发生在花发育后期，即每年春天花开前后的2-3天，雌花中雌蕊的子房开始膨大，柱头伸长，胚珠多数，能育，而雄蕊逐渐干瘪收缩直至最后花粉不育;此时雄花的雄蕊发育良好，花丝粗壮发达，花药饱满，花粉能育，而雌蕊发育逐渐停滞直至败育。花器官发生方式为向心式，发生顺序为：苞片—外轮萼片—内轮萼片—花瓣—外轮雄蕊—内轮雄蕊—雌蕊。花萼4枚，两轮，交互对生，花瓣和雄蕊均4枚一轮几乎同时发生，花瓣覆瓦状排列。旌节花属常绿组和落叶组植物的花器官发生发育过程没有差别。 　　4、基于多个DNA区段序列的系统发育分析 　　首次选取进化速率较快的5个DNA区段：核核糖体ITS区和叶绿体rps16, ndhF，trnL-F及trnS-G-G区段，对旌节花属分别代表常绿组和落叶组的10种2变种及外类群Crossosomataceae和省沽油科3种植物的相关序列进行测定，并将上述序列合并构建系统树。研究结果表明：旌节花科作为一个单系群得到很好的支持。旌节花属落叶组和常绿组的划分在基于ITS、叶绿体4个区段单独分析及联合分析得到的系统树中没有得到支持，表现在该联合分析得到的树的两个分支与该属植物的习性和外部形态特征无关，而与它们的地理分布密切相关：其中一支主要分布在东亚的偏北部，另一支主要分布在东亚的南部地区，常绿和落叶的种类则混合在这两支中。此外，分子序列分析的结果对该属植物的分类、系统关系、起源和该科的系统位置均有重要意义。系统树显示出分布于台湾的S. sigeyosii与分布于日本的早春旌节花聚成稳定的一支，说明前者与后者，而不是与西域旌节花关系密切。叶绿体基因树与核基因树存在冲突，表明该属经历过网状进化。推测变种披针叶旌节花与云南旌节花可能是杂交起源。根据叶绿体rbcL序列得出的系统树显示旌节花科与产北美的Crossosomataceae的姊妹群关系支持率低于50%。基于最大似然树，运用罚分似然法推算出这两个科分化时间为68.25 ± 10.36 MYA，这两者共同与省沽油科近缘。 　　5、分类学研究 　　综合上述研究结果，对旌节花属进行了分类学修订。属下不设立组;确认该属植物有6种6变种，其中中国特有种3个，变种3个。5种5变种被作为新异名，新组合3个，其中2个改级新组合。对西域旌节花、云南旌节花和滇缅旌节花等进行了模式指定。本文还给出了分种和变种检索表、各类群的形态学描述、地理分布图、标本引证和讨论等。 　　

蕨属的系统发育及谱系地理学研究

蕨属（Pteridium Gled. ex Scop）分布范围极广，几乎遍布于全世界。关于此属的属下分类一直争论不休，主要是因为此属的分布广，形态变异比较大，变异类型之间形态差异不显著，缺乏明确的鉴别性状，或者变异类型之间存在过渡。 本研究利用叶绿体rps4-trnS DNA序列和trnS-G DNA序列构建了世界不同地区蕨属植物的系统关系及中国区域蕨属植物的系统关系，同时构建了世界不同地区蕨属植物的单倍型关系。我们发现，（1）南美洲和澳洲的蕨属植物分化最早;（2）北美洲的蕨属植物遗传变异最丰富，有多条相互独立的进化线（evolutionary lineages）;（3）欧亚大陆和非洲的蕨属植物都与北美洲具有密切关系;（4）大洋洲的蕨属植物同时受南美洲和亚洲蕨属植物的影响;（5）东亚的蕨属植物直接受北美洲蕨属植物的影响，与欧洲和非洲没有直接的关系;（6）非洲和欧洲的蕨属植物关系密切，有着相同的起源。 根据不同地点的蕨属植物的系统关系和单倍型Network关系，我们认为，蕨属植物起源于南美洲，而后向北美洲和澳洲扩散。在北美洲东部和东南部蕨属植物得到繁荣并分别向三个方向扩散与分化。第一个方向是朝北美西部扩散，第二个方向是两次独立扩散到亚洲，第三个方向是扩散到非洲，进而由非洲向欧洲扩散。 根据叶绿体基因所揭示的世界不同地区蕨属植物的发生和发展规律结合前人的研究，我们认为蕨属包含四个物种。Pteridium aquilinum (L.) Kuhn，分布于北美洲、亚洲、非洲和欧洲。该种分布广泛，存在一定的地理分化，可进一步划分8个亚种。Pteridium esculentum (G. Forst.) Cokayne，分布于南美洲，并扩散到澳大利亚，甚至东南亚国家。Pteridium caudatum (L.) Maxon，分布于中美洲和南美洲，是个异源四倍体。Pteridium semihastatum (Wall. ex Ag.) S.B.Andrews，分布于澳洲和东南亚，是个异源四倍体。 通过分子数据和形态分析相结合，作者认为中国蕨属包括两个亚种，即蕨[Pteridium aquilinum (L.) Kuhn subsp.japonicum (Nakai) Á.Löve & D.Löve]和毛轴蕨[Pteridium aquilinum (L.) Kuhn subsp. revolutum (Bl.) X. Q. Chen, stat. nov.]。作者认为，《中国植物志》中记载的食蕨，特有种云南蕨、糙轴蕨、镰羽蕨和长羽蕨并不存在或只是蕨或毛轴蕨的变形，或者毛轴蕨与蕨的杂交个体。 为了今后更进一步研究蕨属植物的系统与进化关系、居群的遗传结构等问题，作者还分离了8个具有多态性的微卫星位点，探讨了微卫星引物开发设计过程中如何提高效率等问题。

山茶属茶组植物分子系统学初步研究

山茶科植物的叶绿体基因和nrDNA基因分别为母系和双亲遗传，二者的联合分析为系统发育重建，特别是揭示网状进化过程提供了便利条件。茶组(Sect. Thea (L.) Dyer )属山茶科(Theaceae)山茶属(Camellia)，除茶（C. sinensis (L.) O. Kuntze）为广布种外，其余种主要分布在我国云南、广西、贵州和四川等地。自从W. T. Dyer（1874）建立茶组以来，该组作为一个自然类群，不存在大的争议，但是组下系统分类至今在国内外学术界中争议颇多。 本文选择来自于叶绿体的2个DNA片段，通过序列分析对茶组的28个种进行了系统发育重建。主要结果如下： (1) 叶绿体DNA rpS12-rpL20和trnS-G的序列分析 2个叶绿体DNA片段联合分析，将茶组植物分为两大支，部分种之间的系统发育关系得到了分辨。该叶绿体基因的分析结果支持茶组植物为一个单系类群，但与以往将茶组分为子房五室类群和三室类群的观点不相吻合。 (2) 茶组nrDNA ITS区研究结果 选用14个种茶组植物，对其nrDNA ITS区进行序列分析。茶组植物nrDNA ITS区扩增比较困难，且存在种内多态性，因而在测序之前需要对PCR产物进行克隆。

粗茎秦艽与西藏秦艽(龙胆科)种间的形态学、染色体与分子界定

粗茎秦艽Gentiana crassicaulis Duthie ex Burk和西藏秦艽G. tibetica King ex Hook. f.的主要区别是前者花冠小、檐部蓝色, 后者花大、檐部紫褐色; 但二者在西藏东南部地区的形态界限不清楚, 具有较多的中间形态过渡类型.本文以该地区居群作为主要研究对象, 并以云南和西藏其他地区的典型形态类群为对照, 比较研究了两个物种间的形态学、染色体和分子标记界定.共研究了10个居群69个个体的trnS-G和ITS序列、6个居群26株植株的染色体倍性和54份标本的形态变异.结果表明, 花冠大小与染色体倍性密切相关.云南等地分布的粗茎秦艽均为二倍体(2n=26), 花冠较小; 西藏西部的西藏秦艽均为四倍体(2n=52), 花冠较大.而西藏东南部的形态过渡类型也为四倍体, 花冠较大, 且其中的部分植株花颜色变浅, 类似西藏秦艽.但母系与双亲遗传背景的分子标记都表明西藏东南部居群与典型形态特征上的粗茎秦艽聚为一支, 应属于粗茎秦艽; ITS序列上也未发现杂合位点.这些分子证据都不支持近期杂交渗入造成形态过渡性的假设.综合分析表明多倍化造成了粗茎秦艽种内形态特征的变异, 其四倍体个体与四倍体西藏秦艽在形态特征上的趋同, 最终导致了二者在西藏东南部地区形态界限的模糊.但是, 分子证据可以很好地区分这种多倍化造成的形态趋同.我们的研究结果还表明, 综合利用分子、染色体和广义形态证据来解决近缘种之间种范围的界定具有更大的优势.此外, 还讨论了粗茎秦艽不同染色体倍性的分布以及该高山物种可能的起源与扩散途径.

Direct numerical test of the B-G-K model equation by the DSMC method

In the present paper the rarefied gas how caused by the sudden change of the wall temperature and the Rayleigh problem are simulated by the DSMC method which has been validated by experiments both in global flour field and velocity distribution function level. The comparison of the simulated results with the accurate numerical solutions of the B-G-K model equation shows that near equilibrium the BG-K equation with corrected collision frequency can give accurate result but as farther away from equilibrium the B-G-K equation is not accurate. This is for the first time that the error caused by the B-G-K model equation has been revealed.

Immunosensor Interface Based on Physical and Chemical Immunogylobulin G Adsorption onto Mixed Self-Assembled Monolayers

An immunosensor interface based on mixed hydrophobic self-assembled monolayers (SAMs) of methyl and carboxylic acid terminated thiols with covalently attached human Immunoglobulin G (hIgG), is investigated. The densely packed and organised SAMs were characterised by contact angle measurements and cyclic voltammetry. The effect of the non-ionic surfactant, Tween 20, in preventing nonspecific adsorption is addressed by ellipsometry during physical and covalent hIgG immobilization on pure and mixed SAMs, respectively. It is clearly demonstrated that nonspecific adsorption due to hydrophobic interactions of hIgG on methyl ended groups is totally inhibited, whereas electrostatic/hydrogen bonding interactions with the exposed carboxylic groups prevail in the presence of surfactant. Results of ellipsometry and Atomic Force Microscopy, reveal that the surface concentration of covalently immobilized hIgG is determined by the ratio of COOH/CH3-terminated thiols in SAM forming solution. Moreover, the ellipsometric data demonstrates that the ratio of bound anti-hIgG/hIgG depends on the density of hIgG on the surface and that the highest ratio is close to three. We also report the selectivity and high sensitivity achieved by chronoamperometry in the detection of adsorbed hIgG and the reaction with its antibody.

Numerical simulation of g-jitter effects on half floating zone convection under microgravity environment

The g-jitter effects on the thermocapillary convection in liquid bridge of floating half zone were studied by numerical simulation for unsteady and axi-symmetric model in the cylindrical coordinate system. The g-jitter field was given by a steady microgravity field in addition to an oscillatory low-gravity field, and the effects on the flow field, temperature distribution and free surface deformation were analyzed numerically.

G-jitter effects on half floating-zone convection in intermediate-frequency range

The g-jitter influence on thermocapillary convection and critical Marangoni number in a liquid bridge of half-floating rone was discussed in the low frequency range of 0.4 to 1.5 Hz in a previous paper. This paper extended the experiments to the intermediate frequency range of 2 to 18 Hz, which htrs often been recorded as vibration environment of spacecrafts. The experiment was completed on the deck of a vibration machine, which gave a periodical applied acceleration to simulate the effects of g-jitter. The experimental results in the intermediate frequency range are different from that in the low frequency range. The velocity field and the shape of the free surface have periodical fluctuations in response to g-jitter. The amplitude of the periodical varying part of the temperature response decreases obviously with increasing frequency of g-jitter and vanishes almost when the frequency of g-jitter is high enough. The critical Marangoni number is defined to describe the transition from a periodical convection in response to g-jitter to an oscillatory convection due to internal instability, and will increase with increasing g-jitter frequency. According to the spectral analysis, it can be found that the oscillatory part of temperature is a superposition of two harmonic waves if the Marangoni number is larger than a critical value.

Effects of g-jitter on the critical marangoni number

A half floating zone is fixed on a vibrational deck, which supports a periodical applied acceleration to simulate the effect of g-jitter. This paper deals with the effects of g-jitter on the fluid fields and the critical Marangoni number, which describes the transition from a forced oscillation of thermocapillary convection into an instability oscillatory convection in a liquid bridge of half floating zone with top rod heated. The responses of g-jitter field on the temperature profiles and flow pattern in the liquid bridge were obtained experimentally. The results indicated that the critical Marangoni number decreases with the increasing of g-jitter effect and is slightly smaller for higher frequency of g-jitter with fixed strength of applied gravity.

Quantum and variational monte-carlo interaction potentials for li2 (x1-sigma-g+)

Optimized trial functions are used in quantum Monte Carlo and variational Monte Carlo calculations of the Li₂(X ¹Σ⁺_g) potential curve. The trial functions used are a product of a Slater determinant of molecular orbitals multiplied by correlation functions of electron—nuclear and electron—electron separation. The parameters of the determinant and correlation functions are optimized simultaneously by reducing the deviations of the local energy E_L (E_L  Ψ⁻¹_THΨ_T, where Ψ_T denotes a trial function) over a fixed sample. At the equilibrium separation, the variational Monte Carlo and quantum Monte Carlo methods recover 68% and 98% of the correlation energy, respectively. At other points on the curves, these methods yield similar accuracies.

源于G.I.Taylor工作的现代经典物理学

一位科学家的工作提供了一门完整课程的素材，涉及从流体动力学稳定性、湍流到流体电动力学、微生物的运动．

Starting Nozzle Flow Simulation Using K-G Two-Equation Turbulence Model

The starting process of two-dimensional nozzle flows has been simulated with Euler, laminar and k - g two-equation turbulence Navier-Stokes equations. The flow solver is based on a combination of LUSGS subiteration implicit method and five spatial discretized schemes, which are Roe, HLLE, MHLLE upwind schemes and AUSM+, AUSMPW schemes. In the paper, special attention is for the flow differences of the nozzle starting process obtained from different governing equations and different schemes. Two nozzle flows, previously investigated experimentally and numerically by other researchers, are chosen as our examples. The calculated results indicate the carbuncle phenomenon and unphysical oscillations appear more or less near a wall or behind strong shock wave except using HLLE scheme, and these unphysical phenomena become more seriously with the increase of Mach number. Comparing the turbulence calculation, inviscid solution cannot simulate the wall flow separation and the laminar solution shows some different flow characteristics in the regions of flow separation and near wall.