Mechanisms of femtosecond laser-induced breakdown and damage in MgO

Single-shot laser damage threshold of MgO for 40-986 fs, 800 nm laser pulses is reported. The pump-probe measurements with femtosecond pulses were carried out to investigate the time-resolved electronic excitation processes. A theoretical model including conduction band electrons (CBE) production and laser energy deposition was applied to discuss the roles of multiphoton ionization (MPI) and avalanche ionization in femtosecond laser-induced dielectric breakdown. The results indicate that avalanche ionization plays the dominant role in the femtosecond laser-induced breakdown in MgO near the damage threshold. (c) 2005 Elsevier B.V. All rights reserved.

飞秒激光作用下透明介质材料的反射率演化

利用800nm抽运和400nm探针技术，测量了CaF2和MgO的时间分辨反射率，研究了材料的电子激发和弛豫超快动力学过程。采用耦合动力学模型，探讨了飞秒激光对透明介质材料的激发，以及材料的激发对抽运激光在材料中的传输、分布和反射特性的影响。根据这个理论模型计算了时间分辨反射率的演化，计算结果和实验结果相吻合。研究表明，多光子电离（MPI）和碰撞电离（II）在介质材料的导带电子激发中起着重要的作用。

不同功率下的飞秒激光在水中形成的脉冲压缩和超连续谱

本文研究了飞秒激光脉冲在水中的传输情况.通过改变不同的激光输入功率进行模拟，我们发现从输入功率略高于到远远高于发生自聚焦的临界功率,分别是群速度色散和多光子电离多光子吸收阻止了自聚焦导致的脉冲塌陷，当多光子电离和多光子吸收主导传输时，脉冲能被压缩到几个光学周期.在频域，多光子电离能引起很强的蓝移,而多光子吸收能对这种蓝移起到抑制作用。

Interplay of martensitic phase transformation and plastic slip in polycrystals

Inspired by key experimental and analytical results regarding Shape Memory Alloys (SMAs), we propose a modelling framework to explore the interplay between martensitic phase transformations and plastic slip in polycrystalline materials, with an eye towards computational efficiency. The resulting framework uses a convexified potential for the internal energy density to capture the stored energy associated with transformation at the meso-scale, and introduces kinetic potentials to govern the evolution of transformation and plastic slip. The framework is novel in the way it treats plasticity on par with transformation.

We implement the framework in the setting of anti-plane shear, using a staggered implicit/explict update: we first use a Fast-Fourier Transform (FFT) solver based on an Augmented Lagrangian formulation to implicitly solve for the full-field displacements of a simulated polycrystal, then explicitly update the volume fraction of martensite and plastic slip using their respective stick-slip type kinetic laws. We observe that, even in this simple setting with an idealized material comprising four martensitic variants and four slip systems, the model recovers a rich variety of SMA type behaviors. We use this model to gain insight into the isothermal behavior of stress-stabilized martensite, looking at the effects of the relative plastic yield strength, the memory of deformation history under non-proportional loading, and several others.

We extend the framework to the generalized 3-D setting, for which the convexified potential is a lower bound on the actual internal energy, and show that the fully implicit discrete time formulation of the framework is governed by a variational principle for mechanical equilibrium. We further propose an extension of the method to finite deformations via an exponential mapping. We implement the generalized framework using an existing Optimal Transport Mesh-free (OTM) solver. We then model the $\alpha$--$\gamma$ and $\alpha$--$\varepsilon$ transformations in pure iron, with an initial attempt in the latter to account for twinning in the parent phase. We demonstrate the scalability of the framework to large scale computing by simulating Taylor impact experiments, observing nearly linear (ideal) speed-up through 256 MPI tasks. Finally, we present preliminary results of a simulated Split-Hopkinson Pressure Bar (SHPB) experiment using the $\alpha$--$\varepsilon$ model.

Self-compression of femtosecond pulses in argon with a power close to the self-focusing threshold

Self-compression of femtosecond pulses in noble gases with an input power close to the self-focusing threshold has been investigated experimentally and theoretically. It is demonstrated that either multiphoton ionization (MPI) or space time focusing and self-steepening effects can induce pulse shortening, but they predominate at different beam intensities during the propagation. The latter effects play a key role in the final pulse self-compression. By choosing an appropriate focusing parameter, action distance of the space time focusing and self-steepening effects can be lengthened, which can promote a shock pulse structure with a duration as short as two optical cycles. It is also found that, for our calculation cases in which an input pulse power is close to the self-focusing threshold, either group velocity dispersion (GVD) or multiphoton absorption (MPA) has a negligible influence on pulse characteristics in the propagation process.

Nonlinear X-wave formation and conical emission at different powers of a femtosecond laser pulse in water

Nonlinear X-wave formation at different pulse powers in water is simulated using the standard model of nonlinear Schrodinger equation (NLSE). It is shown that in near field X-shape originally emerges from the interplay between radial diffraction and optical Kerr effect. At relatively low power group-velocity dispersion (GVD) arrests the collapse and leads to pulse splitting on axis. With high enough power, multi-photon ionization (NIPI) and multi-photon absorption (MPA) play great importance in arresting the collapse. The tailing part of pulse is first defocused by MPI and then refocuses. Pulse splitting on axis is a manifestation of this process. Double X-wave forms when the split sub-pulses are self-focusing. In the far field, the character of the central X structure of conical emission (CE) is directly related to the single or double X-shape in the near field. (c) 2007 Elsevier B.V. All rights reserved.

Implementação paralela do algoritmo iterativo de busca do parâmetro de regularização ótimo para o funcional de Tikhonov no problema de restauração de imagens

O uso de técnicas com o funcional de Tikhonov em processamento de imagens tem sido amplamente usado nos últimos anos. A ideia básica nesse processo é modificar uma imagem inicial via equação de convolução e encontrar um parâmetro que minimize esse funcional afim de obter uma aproximação da imagem original. Porém, um problema típico neste método consiste na seleção do parâmetro de regularização adequado para o compromisso entre a acurácia e a estabilidade da solução. Um método desenvolvido por pesquisadores do IPRJ e UFRJ, atuantes na área de problemas inversos, consiste em minimizar um funcional de resíduos através do parâmetro de regularização de Tikhonov. Uma estratégia que emprega a busca iterativa deste parâmetro visando obter um valor mínimo para o funcional na iteração seguinte foi adotada recentemente em um algoritmo serial de restauração. Porém, o custo computacional é um fator problema encontrado ao empregar o método iterativo de busca. Com esta abordagem, neste trabalho é feita uma implementação em linguagem C++ que emprega técnicas de computação paralela usando MPI (Message Passing Interface) para a estratégia de minimização do funcional com o método de busca iterativa, reduzindo assim, o tempo de execução requerido pelo algoritmo. Uma versão modificada do método de Jacobi é considerada em duas versões do algoritmo, uma serial e outra em paralelo. Este algoritmo é adequado para implementação paralela por não possuir dependências de dados como de Gauss-Seidel que também é mostrado a convergir. Como indicador de desempenho para avaliação do algoritmo de restauração, além das medidas tradicionais, uma nova métrica que se baseia em critérios subjetivos denominada IWMSE (Information Weighted Mean Square Error) é empregada. Essas métricas foram introduzidas no programa serial de processamento de imagens e permitem fazer a análise da restauração a cada passo de iteração. Os resultados obtidos através das duas versões possibilitou verificar a aceleração e a eficiência da implementação paralela. A método de paralelismo apresentou resultados satisfatórios em um menor tempo de processamento e com desempenho aceitável.

短脉冲激光辐照下熔融石英薄膜的损伤机理

激光诱导薄膜损伤过程中，雪崩离化（AI）和多光子离化（MPI）的性质和作用到目前仍然存在争议。基于STUART等人的电子密度演化方程，运用数值模拟方法，研究了脉宽为τ∈[0．01，5]ps范围内单脉冲激光作用下熔融石英薄膜中电子密度演化过程；讨论了初始电子密度、激光脉冲宽度对阈值功率密度和阈值能量的影响；分析了初始电子密度、激光脉冲宽度对多光子离化及雪崩离化的影响。研究结果表明，在所研究的脉宽范围内，对于熔融石英光学薄膜、飞秒激光诱导损伤以雪崩离化为主导，多光子离化的影响随着脉宽的降低而增强，雪崩离化所需种子电子主要来源于多光子离化。

Um metodo numerico com paralelismo no tempo para aproximar solucoes de EDPs

Este trabalho de pesquisa tem por objetivo apresentar e investigar a viabilidade de um método numérico que contempla o paralelismo no tempo. Este método numérico está associado a problemas de condição inicial e de contorno para equações diferenciais parciais (evolutivas). Diferentemente do método proposto neste trabalho, a maioria dos métodos numéricos associados a equações diferencias parciais evolutivas e tradicionalmente encontrados, contemplam apenas o paralelismo no espaço. Daí, a motivação em realizar o presente trabalho de pesquisa, buscando não somente um método com paralelismo no tempo mas, sobretudo, um método viável do ponto de vista computacional. Para isso, a implementação do esquema numérico proposto está por conta de um algoritmo paralelo escrito na linguagem C e que utiliza a biblioteca MPI. A análise dos resultados obtidos com os testes de desempenho revelam um método numérico escalável e que exige pouco nível de comunicação entre processadores.

Sea level rise adaptation: emerging lessons for local policy development

Many coastal communities across the United States are beginning to plan for climate-related sea level rise. While impacts and solutions will vary with local conditions, jurisdictions which have begun this process seem to pass through three common stages when developing policy for local sea level rise adaptation: l) building awareness about local sea level rise threats, 2) undertaking analyses of local vulnerabilities, and 3) developing plans and policies to deal with these vulnerabilities. The purpose of this paper is to help advance community dialogue and further inform local decision-makers about key elements and steps for addressing climate-related sea level rise. It summarizes the results of a project the Marine Policy Institute (MPI) undertook during 2011-12 to review experiences from fourteen U.S. coastal jurisdictions representing a variety of city, county, and state efforts with sea level adaptation. There are many more initiatives underway than those reflected in this sample, but the “focus jurisdictions” were selected because of the extensive information publically available on their experiences and lessons being learned that could provide insights for coastal communities, especially in Southwest Florida.

An examination of spatial and temporal genetic variation in walleye pollock (Theragra chalcogramma) using allozyme, mitochondrial DNA, and microsatellite data

We used allozyme, microsatellite, and mitochondrial DNA (mtDNA) data to test for spatial and interannual genetic diversity in wall-eye pollock (Theragra chalcogramma) from six spawning aggregations representing three geographic regions: Gulf of Alaska, eastern Bering Sea, and eastern Kamchatka. Interpopulation genetic diversity was evident primarily from the mtDNA and two allozyme loci (SOD-2*, MPI*). Permutation tests ˆindicated that FST values for most allozyme and microsatellite loci were not significantly greater than zero. The microsatellite results suggested that high locus polymorphism may not be a reliable indicator of power for detecting population differentiation in walleye pollock. The fact that mtDNA revealed population structure and most nuclear loci did not suggests that the effective size of most walleye pollock populations is large (genetic drift is weak) and migration is a relatively strong homogenizing force. The allozymes and mtDNA provided mostly concordant estimates of patterns of spatial genetic variation. These data showed significant genetic variation between North American and Asian populations. In addition, two spawning aggregations in the Gulf of Alaska, in Prince William Sound, and off Middleton Island, appeared genetically distinct from walleye pollock spawning in the Shelikof Strait and may merit management as a distinct stock. Finally, we found evidence of interannual genetic variation in two of three North American spawning aggregations, similar in magnitude to the spatial variation among North American walleye pol-lock. We suggest that interannual genetic variation in walleye pollock may be indicative of one or more of the following factors: highly variable reproductive success, adult philopatry, source-sink metapopulation structure, and intraannual variation (days) in spawning timing among genetically distinct but spatially identical spawning aggregates.

全球变化下中国东北样带植物多样性梯度分析及气候-NDVI-植被的关系研究

中国东北样带（NorthEast China Transect, NECT）是位于中纬度温带以降水量作为主要驱动因素的陆地样带。本文的工作以此作为研究平台，利用生态信息系统（Ecological Information System, EIS）以及Microsoft Excel 7.0软件包建立了样带的地理数据库和植物多样性数据库，包括气候数据库、植被数据库、遥感数据库和内蒙内C_4植物数据库以及样带内生态系统特征数据库。在此基础上，主要研究了以下四个方面的内容： 1. 利用Holdridge的生命地带方法对NECT内的生物群区进行了划分。 主要是确定了生物群区间过渡带的位置与宽度，并预测了在全球变化三种模式下NECT内生物群区，尤其是过渡带的变化图景。湿度升高2 ℃后，过渡带的面积都呈扩大化的趋势。森林区对于降水量的变化反应很敏感。荒漠灌丛（即荒漠草原类型）由于其水热条件处于样带内较极端类型，因而对于全球气候变化反应也比较敏感。 2. 研究了NECT内的α、β多样性以及包括生活型、水分生态型、区系地理成分等在内的植物群落特征多样性的梯度变化规律。 研究了样带内的多样性梯度，提出了在样带内存在的α多样性测度问题以及β多样性沿样带的变化规律：样带内由东到西，β多样性逐渐升高，群落内物种被替代的速率变慢;两种植被类型边界上的两个样地之间的相似程度由东到西呈上升趋势;同一类型群落之间的物种周转率比不同类型群落间的物种周转率相对要低。同时将各个环境因子与α、β多样性作了回归分析，找出样带内决定α、β多样性的主要环境因子指标。 样带内沿43.5°N一线附近植物群落的生活型共有17类，水分生态型8类，区系地理成分包括17类，以此为基础分析了群落特征沿样带的变化规律。并探讨了生活型分布的历史地理原因。 3. 对样带气候－NDVI间的关系以及植被－NDVI的关系进行了探讨。 利用来自气象卫星的遥感数据一归一化植被指数（NDVI），和数值化后的样带1:100万植被图进行叠加，找到NECT内每种植被类型对应的NDVI值。样带内共有植被类型147种，反映在NDVI变化上的植被类型有106类。其中，自然植被101种。 影响年均NDVI分布的因子主要有经度、辐射日照百分率及7月温度，与经度呈正相关，与辐射日照时数及7月温度呈负相关。回归方程如下： NDVI = -220.426 + 3.273Lon - 80.338Ratio - 1.962T_7 (R~2 = 0.9714, F = 521.52, p < 0.001) 4. 研究了NECT内的光合功能型。 主要包括内蒙古地区的C_4植物及其生态地理特性。揭标C_4植物的分类群特性、生活型、水分生态型与区系地理成分等生态学特性。C_4植物分布的科属极其集中。C_4光合型为维管植物某些分类群（科、属、种）的特性，为它们固有的遗传特性。推断C_4起源于草本的某些科属。C_4植物为喜热、耐旱的类群。世界种、泛热带种、泛地中海种C_4植物较集中。 样带内的C_3、C_4功能型及其与环境因子的相关性。样带内C_4和C_3光合型植物组成比例由东到西表现出两高两低的趋势。分布主要与年均温和降水量呈显著相关。 提出了一种新的C_3、C_4鉴别方法。即根据野外测定的光合数据建立了C_3、C_4的判别模型： f_1(x) = -1.5493 + 0.1427Pn + 0.1035Tr + 0.3768ΔT + 0.1000Gs f_2(x) = -15.6142 + 1.0542Pn - 0.2503Tr - 0.2957ΔT + 0.6491Gs 最后，综合7个GCMs模型（GFDL，GISS，LLNL，MPI，OSU，UKMOH，UKMOL）的输出结果，利用此结果和本文建立的回归模型，模拟了样带内生物多样性的窨分布格局，并预测了末来全球变化下归一化植被指数NDVI的空间分布格局的变化。

基于机群架构的并行数据库中间件系统改进研究

中国计算机学会

FM-index分块并行算法及其实现

查询海量数据有压缩和索引两种方法来提高速度,该文结合这两种方法提出了压缩查询的方法.FM-index是一种自索引的全文查询算法,存在内存占用过大的问题,对于复杂的查询效率也不理想.该文提出分块FM-index算法,在分块的基础上采用MPI对算法进行并行化,解决了内存占用过多的问题,达到了较好的并行效率.

并行最优路径算法及并行层次通信模型LogGP(h)研究

随着以数据为中心的超级计算时代的到来，在各种以图为数据结构的应用中数据规模日益增大，数据量的急剧增加使得串行最优路径算法成为应用的性能瓶颈，已不能满足大规模最优路径求解的实时性要求。然而，最优路径求解算法的并行化却对图的类型有高度的敏感性，尤其是对一些稀疏图，如交通网络图，目前并没有很好的算法提供足够的并行性[1]。对此，本文提出、实现和优化了一种新的基于网络划分和迭代更新的并行最优路径算法，并用真实的交通网络数据在IBM JS21集群和曙光5000A两种平台上对算法进行了评测。测试结果表明，该算法具有较好的性能，在IBM JS21集群上，使用16个进程会出现约15倍的加速比；而在曙光5000A上，算法使用同一节点内的16个核获得了约20倍的超线性加速比。 本文的另一部分研究工作集中在新的并行计算模型方面。研究者所提出的并行计算模型有很多，从基于共享存储器结构的第一代并行计算模型，到针对于分布式存储结构的第二代并行模型，再到考虑存储层次性的第三代并行计算模型，却较少有模型考虑到通信的层次性。 大规模集群系统往往由几百到几千、甚至上万个节点构成，网络拓扑结构复杂。在这样的大规模集群系统中，点对点通信的性能并不总是一致的。这种非一致的通信性能随着网络拓扑的复杂化呈现层次性的变化，本文首先提出了层次化非一致性带宽和延迟(Hierarchical Communication Latency and Bandwidth：HCLB)的概念，并以进程映射对MPI Allgather的四种算法性能影响说明了由多核结构引起的通信层次性；之后，本文将大规模集群系统中的HCLB现象模型化，提出LogGP(h)并行通信模型。LogGP(h)在LogGP模型基础上，通过参数h将网络拓扑结构抽象到模型中，很好的表达了通信层次性。本文在具有8个刀片的IBM JS21集群上对点对点通信和集合通信开销进行模型分析和实际测试，实验结果表明，LogGP(h)较未引入参数h的LogGP模型更精确。