Effective AC and DC responses of nonlinear composites at higher concentration

For higher concentration of inclusions, an effective medium approximation (EMA) method is used to investigate the bulk effective response of weakly nonlinear composites, which are subject to the constitutive relation of electric displacement and electric field, D-alpha = epsilon E-alpha + chi(alpha)|E|(2) E. As an example of three dimensions, under the external AC and DC electric fields E-app = E-a (1 + sin omega t), we have derived the general effective nonlinear response of composites by the EMA method for higher concentration of spherical inclusions. Furthermore, the effective nonlinear responses at harmonics are predicted.

Dielectric response of spherically anisotropic graded piezoelectric composites

A graded piezoelectric composite consisting of a spherically anisotropic graded piezoelectric inclusion imbedded in an infinite nonpiezoelectric matrix, with the physical properties of the graded spherical inclusion having a power-law profile with respect to the radial variable r, is studied theoretically. Under an external uniform electric field, the electric displacement field and the elastic stress tensor field of this spherically anisotropic graded piezoelectric composite are derived exactly by means of displacement separation technique, based on the governing equations in the dilute limit. A piezoelectric response mechanism, in which the effective piezoelectric response vanishes along the z direction (or x,y directions), is revealed in this kind of graded piezoelectric composites. Furthermore, it is found that the effective dielectric constant decreases (or increases) with the volume fraction p of the inclusions if the exponent parameter k of the grading profile is larger (or smaller) than a critical value. (C) 2007 American Institute of Physics.

基于介电泳的电极阵列电场仿真研究

介电泳方法被广泛地应用于微纳颗粒的分离和操纵中,实现介电泳操作的关键是设计满足所需电场分布的电极阵列。针对目前在微电极阵列设计中尚缺乏简单有效的电场解析方法的现状,提出一种基于格林公式的电极阵列电场的解析方法。首先介绍了传统介电泳和行波介电泳的概念和计算模型,分析了介电泳过程与电极上所施加的交变电压的频率和幅度的关系,然后在确立电极电势的边界条件的基础上,采用基于格林公式的电场解析方法,建立了非均匀电场的解析模型,得出不同条件下的电极阵列电场分布的仿真结果,最后利用FEMLAB有限元仿真软件对解析模型进行了对比仿真,验证了该解析模型的可行性。基于格林公式的电场解析求解方法能够有效地提高电极阵列设计中的针对性以及缩短电极设计的时间。

面向微粒操纵的介电泳芯片系统研究

针对目前在纳米器件及传感器的制造中尚无对大量粒子进行有效操纵的方法,我们利用介电泳方法对大量微粒进行定位和传输操纵,介绍了利用MEMS工艺进行介电泳芯片加工的过程以及整个观测与实验系统的建立,通过有限元软件对传统介电泳和行波介电泳中电极阵列的电场分布进行求解,并在该实验系统下实现了对微通道中的悬浮高度和微粒的运动速度的测量.该实验系统的研究为液体环境下微纳颗粒的装配和分离提供了一条有效的技术路径.

基于机械力与电场的AFM纳米操作/加工研究

纳米技术是在纳米尺度上研究物质的特性和相互作用，以及利用这些特性的多学科交叉的新兴科学技术。纳米加工技术是对纳米技术的重要应用，它使人类在纳米尺度上进行结构和器件的制造成为可能。基于原子力显微镜(AFM)的纳米加工是近20年来发展的技术，由于AFM具有许多新颖的技术特点(如高分辨率，可操控等)，因而基于AFM的纳米加工技术成为纳米科学研究的新热点。 本论文以纳米制造为背景，围绕基于AFM的纳米加工技术所涉及的科学问题，重点开展了基于机械力的AFM纳米操作和基于电场的AFM纳米加工的理论方法和实现技术研究。在这些研究的基础上，开展了基于纳观机械力的刻划、纳米颗粒的排列、碳纳米管的装配实验研究，开展了基于微观电场的氧化点、氧化线、氧化文字的加工实验研究，以及微观电场氧化切割、焊接加工方法的研究。 本文的主要内容包括：交互式纳米操作系统的构建，基于机械力的AFM纳米操作研究，基于电场的AFM纳米加工研究，基于机械力与电场的AFM纳米操作/加工实验研究。具体如下： 在商用AFM的基础上，研究并搭建了具有力觉与视觉反馈的交互式纳米操作系统。通过该系统操作者可以实时感觉到操作中作用在AFM针尖的纳观力并控制针尖的运动，还可以实时观察到模拟出的针尖操作过程，因而显著提高了纳米操作的直观性、可靠性及效率。 在基于机械力的AFM纳米操作研究方面，为了理解纳观环境下普遍存在的黏附力特性，对AFM的力曲线进行了分析并给出了详细解释；通过研究探针悬臂的形变，提出了探针所受三维作用力的建模方法，以此模型获得了操作中的纳观力信息；基于多刚体动力学理论提出了碳纳米管的二连杆动力学建模方法，并对碳纳米管在操作中的弯曲特性进行了仿真分析。这些研究为基于机械力的AFM纳米操作提供了理论指导。 在基于电场的AFM纳米加工方面，研究了微观电场在样本表面的分布，建立了电场氧化加工中电场分布模型，以此模型分析了氧化物特性与电场分布的关系；对氧化加工中氧化物生长过程进行了动力学分析，得到了氧化物生长的理论模型；系统地研究了氧化加工中各因素(偏压，加工速度，针尖－样本距离等)对氧化结构生成的影响，总结了氧化物的生成与各因素的关系，验证了理论分析的结果，提高了电场氧化加工的重复性及可控性。这些研究为基于电场的AFM纳米加工提供了理论及实验依据。 在基于机械力与电场的AFM纳米操作/加工理论研究的基础上，进行了相关的实验研究。通过基于机械力的AFM纳米操作方法实现了在样本表面刻划文字，纳米颗粒的排列，碳纳米管的操作及装配；通过基于电场的AFM纳米加工方法实现了在样本表面加工氧化点，氧化线及氧化文字等纳米结构。利用电场氧化加工方法还实现了对碳纳米管的定点切割与焊接，拓展了AFM纳米加工的应用领域。这些实验研究验证了系统的操作及装配性能，证明了电场加工的有效性及具有潜在的应用前景。 本文的研究工作为制造纳米结构、纳米器件提供了一条可行的技术途径，对AFM纳米加工技术的进一步发展具有一定的理论参考及实践指导意义。

空间电流体系地磁反演及理论研究

The space currents definitely take effects on electromagnetic environment and also are scientific highlight in the space research. Space currents as a momentum and energy provider to Geospace Storm, disturb the varied part of geomagnetic field, distort magnetospheric configuration and furthermore take control of the coupling between magnetosphere and ionosphere. Due to both academic and commercial objectives above, we carry on geomagnetic inverse and theoretical studies about the space currents by using geomagnetic data from INTERMAGNET. At first, we apply a method of Natural Orthogonal Components (NOC) to decomposition the solar daily variation, especially for (solar quiet variation). NOC is just one of eign mode analysis, the most advantage of this method is that the basic functions (BFs) were not previously designated, but naturally came from the original data so that there are several BFs usually corresponding to the process really happened and have more physical meaning than the traditional spectrum analysis with the fixed BFs like Fourier trigonometric functions. The first two eign modes are corresponding to the and daily variation and their amplitudes both have the seasonal and day-to-day trend, that will be useful for evaluating geomagnetic activity indices. Because of the too strict constraints of orthogonality, we try to extend orthogonal contraints to the non-orthogonal ones in order to give more suitable and appropriate decomposition of the real processes when the most components did not satisfy orthogonality. We introduce a mapping matrix which can transform the real physical space to a new mathematical space, after that process, the modified components which associated with the physical processes have satisfied the orthogonality in the new mathematical space, furthermore, we can continue to use the NOC decomposition in the new mathematical space, and then all the components inversely transform back to original physical space, so that we would have finished the non-orthogonal decomposition which more generally in the real world. Secondly, geomagnetic inverse of the ring current’s topology is conducted. Configurational changes of the ring current in the magnetosphere lead to different patterns of disturbed ground field, so that the global configuration of ring current can be inferred from its geomagnetic perturbations. We took advantages of worldwide geomagnetic observatories network to investigate the disturbed geomagnetic field which produced by ring current. It was found that the ring current was not always centered at geomagnetic equator, and significantly deviated off the equator during several intense magnetic storms. The deviation owing to the tilting and latitudinal shifting of the ring current with respect to the earth’s dipole can be estimated from global geomagnetic survey. Furthermore those two configurational factors which gave a quantitative description of the ring current configuration, will be helpful to improve the Dst calibration and understand the dependence of ring current’s configuration on the plasma sheet location relative to the equator when magnetotail field warped. Thirdly, the energization and physical acceleration process of ring current during magnetic storm has been proposed. When IMF Bz component increase, the enhanced convection electric field drive the plasma injection into the inner magnetosphere. During the transport process, a dynamic heating is happened which make the particles more ‘hot’ when the injection is more deeply inward. The energy gradient along the injection path is equivalent to a kind of force, which resist the plasma more earthward injection, as a diamagnetic effect of the magnetosphere anti and repellent action to the exotically injected plasma. The acceleration efficiency has a power law form. We use analytical way to quantitatively describe the dynamical process by introducing a physical parameter: energization index, which will be useful to understand how the particle is heated. At the end, we give a scheme of how to get the from storm time geomagnetic data. During intense magnetic storms, the lognormal trend of geomagnetic Dst decreases depend on the heating dynamic of magnetosphere controlling ring current. The descending pattern of main phase is governed by the magnetospheric configuration, which can be describled by the energization index. The amplitude of Dst correlated with convection electric field or south component of the solar wind. Finally, the Dst index is predicted by upstream solar wind parameter. As we known space weather have posed many chanllenges and impacts on techinal system, the geomagnetic index for evaluating the activity space weather. We review the most popular Dst prediction method and repeat the Dst forecasting model works. A concise and convnient Key Points model of the polar region is also introduced to space weather. In summary, this paper contains some new quantitative and physical description of the space currents with special focus on the ring current. Whatever we do is just to gain a better understanding of the natural world, particularly the space environment around Earth through analytical deduction, algorithm designing and physical analysis, to quantitative interpretation. Applications of theoretical physics in conjunction with data analysis help us to understand the basic physical process govering the universe.

磁层-电离层耦合现象的分析研究

The magnetosphere-ionosphere coupling is mainly manifested by the trans- porting processes of energy into the ionosphere , the energy is carried by solar wind and firstly accumulate at the magnetosphere, and the coupling processes also significantly include the interaction between the magnetosphere and ionosphere for mass and energy. At the quiet condition, energy is delivered by the large-scale convection of the geomagnetic field; the huge energy from solar wind bulk will be injected into and consumed at the near magnetosphere and ionosphere by the geomagnetic storm and substorm activities. Aurorae and FACs (Field-aligned currents) are the important phenomena in the coupling processes. In the present work, firstly, we analyze the activity characteristics of auroral precipitating particle, secondly, we study the distribution characters of large-scale field aligned currents (LS FACs) at storm-time using the observations from different satellites at different altitudes. Finally, we investigate the evolution of the geomagnetic field configuration at the nightside sector on the onset of the expansion phase in a substorm event, the substorm event happened at 0430UT to 0630UT on 8th Nov. 2004. The main results as follows： At the first, the data of the estimated power input (EPI) of auroral particles from NOAA/POES (Polar orbiting environmental satellite) for some 30 years have been analyzed. The variation tendencies of the EPI generally coincide with aa, AE and Dst indices. The annual variation of EPI shows equinox peaks and an asymmetric-activity with a higher peak in the winter-hemisphere than in the summer-hemisphere. The diurnal UT variations are different from north and south hemisphere: for north hemisphere, the peak appears at 1200UT, and the relative deviation is 22% to the daily average of the north hemisphere. For south hemisphere, the maximal deviation is 22% at 2000UT. So the diurnal variation of EPI is more dominant than the annual variation which maximal deviation is 3% to 12% for different seasons. Studies on correlations of the hourly average of EPI, Pa, with AE and Dst indices show a correlation coefficient r=0.74 of Pa and AE, and r=-0.55 of Pa and Dst. The hourly EPIs for north and south polar regions, NPa and SPa, show a north-south asymmetry with a higher correlation of SPa and AE (or Dst). Time delays of EPI with respect to magnetic indices are examined, the maximum correlation coefficient of Pa with AE (r=0.78) occurs when the time delay =0, suggesting a synchronous activity of auroral electrojet and auroral precipitating particles, while =1-2h, the correlation coefficient of Pa with Dst is maximum (r=0.57), suggesting that the activity of auroral particle precipitating may influence the ring current on some extent. Sencondly, we use the high-resolution magnetic field vector data of the CHAMP satellite to investigate the distribution of large-scale FACs during the great magnetic storm on 7th to 8th Nov. 2004. The results show that, whether in the northern or southern hemisphere, the number and density of large-scale FACs during the main-phase are more and bigger than these during the recover-phase, and the number of large-scale FACs in morning sector obviously is more than that in afternoon sector. In terms of the magnetic indices, we find that large-scale FACs in morning sector significantly affected by the substorm activities, while in afternoon sector the large-scale FACs mainly indicate the fluctuations of the ring-current in storm time. Accordingly to the former studies, similarly, we find that in the morning sector, the scale of the large-scale FACs move to the high-latitude region, and in the afternoon sector, large-scale FACs distinctly expand to the low-latitude region. During the time periods that the NOAA/POES auroral precipitating particle power data temporally correspond to the large-scale FACs, the more the power of auroral particle is, the more and bigger the number and density of FACs are. At the same time, we use the magnetic field vector data of POLAR obtain a good form of region 1, region 2, and three pieces of cusp FACs during a single transit at 1930UT-2006UT on 07th. And the characteristics of simultaneous electric field and energy particles observations on Polar are coincide with the five FACs pieces. Finally, by means of the observation of Cluster 4 and Goes 10、 Goes 12, we analyze the evolution process of the change of the magnetic field configuration at night sector at the expansion phase of a substorm event which happened during 0430UT to 0630UT on 8th Nov. 2004, we find that the times of the beginning of the polarizations of magnetic field are observed from Goes 10 to Goes 12 then to Cluster 4. So, at the synchronous orbit ( 6.6 RE) to 10RE distance scale of the neutral sheet, the current disruption spread tailward. Simultaneously, the strengthen of the FACs deduced from these satellites’ magnetic field observations are almost consistent with the times of polarizations, as well as the high energy particles injection and the electric field dominant variation. The onset times determined by the magnetic field polarizations from these satellites are all ahead of the onset time that confirmed from the auroral electrojet indices. So, these characters of different observations can be used as the criterions to determine the onset time for the substorms of such type as we studied.

A collision model for DNA separation by capillary electrophoresis in dilute polymer solution

A theoretical description. based on chemical kinetics and electrochemistry, is given of DNA separation in dilute polymer solution by capillary electrophoresis. A self-consistent model was developed leading to predictions of the DNA electrophoretic velocity as a function of the experimental conditions - polymer concentration, temperature, and electric field strength. The effect of selected experimental variables is discussed. The phenomena discussed are illustrated for the example of 100 bp DNA ladder separation in dilute HPMC solution by capillary electrophoresis. This model is the first single model that can fully explain the dependence of DNA electrophoretic velocity on electrophoretic conditions.

Modeling subauroral polarization streams equatorward of the plasmapause footprints

Pint?r, B.; Thom, S. D.; Balthazor, R.; Vo, H.; Bailey, G. J., Modeling subauroral polarization streams equatorward of the plasmapause footprints, Journal of Geophysical Research, Volume 111, Issue A10, CiteID A10306 RAE2008

Nowe metody wyznaczania parametru asymetrii tensora gradientu pola elektrycznego w spektroskopii NQR i ich zastosowanie do badania układów nieuporządkowanych

Wydział Fizyki

Dynamical Instabilities and Deterministic Chaos in Ballistic Electron Motion in Semiconductor Superlattices

We consider the motion of ballistic electrons within a superlattice miniband under the influence of an alternating electric field. We show that the interaction of electrons with the self-consistent electromagnetic field generated by the electron current may lead to the transition from regular to chaotic dynamics. We estimate the conditions for the experimental observation of this deterministic chaos and discuss the similarities of the superlattice system with the other condensed matter and quantum optical systems.

Symmetry-breaking and chaos in electron transport in semiconductor superlattices

We study the motion of electrons in a single miniband of a semiconductor superlattice driven by THz electric field polarized along the growth direction. We work in the semiclassical balance-equation model, including different elastic and inelastic scattering rates, and incorporating the self-consistent electric field generated by electron motion. We explore regions of complex dynamics, which can include chaotic behaviour and symmetry-breaking. We estimate the magnitudes of dc current and dc voltage that spontaneously appear in regions of broken-symmetry for parameters characteristic of modern semiconductor superlattices. This work complements PRL 80(1998)2669 [ cond-mat/9709026 ].

Dissipative Chaos in Semiconductor Superlattices

We consider the motion of ballistic electrons in a miniband of a semiconductor superlattice (SSL) under the influence of an external, time-periodic electric field. We use the semi-classical balance-equation approach which incorporates elastic and inelastic scattering (as dissipation) and the self-consistent field generated by the electron motion. The coupling of electrons in the miniband to the self-consistent field produces a cooperative nonlinear oscillatory mode which, when interacting with the oscillatory external field and the intrinsic Bloch-type oscillatory mode, can lead to complicated dynamics, including dissipative chaos. For a range of values of the dissipation parameters we determine the regions in the amplitude-frequency plane of the external field in which chaos can occur. Our results suggest that for terahertz external fields of the amplitudes achieved by present-day free electron lasers, chaos may be observable in SSLs. We clarify the nature of this novel nonlinear dynamics in the superlattice-external field system by exploring analogies to the Dicke model of an ensemble of two-level atoms coupled with a resonant cavity field and to Josephson junctions.

A microscopic study of structural and electronic properties of functionalized silicon surfaces based on first-principles

Surface modification of silicon with organic monolayers tethered to the surface by different linkers is an important process in realizing future (opto-)electronic devices. Understanding the role played by the nature of the linking group and the chain length on the adsorption structures and electronic properties of these assemblies is vital to advance this technology. This Thesis is a study of such properties and contributes in particular to a microscopic understanding of induced changes in the work function of experimentally studied functionalized silicon surfaces. Using first-principles density functional theory (DFT), at the first step, we provide predictions for chemical trends in the work function of hydrogenated silicon (111) surfaces modified with various terminations. For nonpolar terminating atomic species such as F, Cl, Br, and I, the change in the work function is directly proportional to the amount of charge transferred from the surface, thus relating to the difference in electronegativity of the adsorbate and silicon atoms. The change is a monotonic function of coverage in this case, and the work function increases with increasing electronegativity. Polar species such as −TeH, −SeH, −SH, −OH, −NH2, −CH3, and −BH2 do not follow this trend due to the interaction of their dipole with the induced electric field at the surface. In this case, the magnitude and sign of the surface dipole moment need to be considered in addition to the bond dipole to generally describe the change in work function. Compared to hydrogenated surfaces, there is slight increase in the work function of H:Si(111)-XH, where X = Te, Se, and S, whereas reduction is observed for surfaces covered with −OH, −CH3, and −NH2. Next, we study the hydrogen passivated Si(111) surface modified with alkyl chains of the general formula H:Si–(CH2)n–CH2 and H:Si–X–(CH2)n–CH3, where X = NH, O, S and n = (0, 1, 3, 5, 7, 9, 11), at half coverage. For (X)–Hexyl and (X)–Dodecyl functionalization, we also examined various coverages up to full monolayer grafting in order to validate the result of half covered surface and the linker effect on the coverage. We find that it is necessary to take into account the van der Waals interaction between the alkyl chains. The strongest binding is for the oxygen linker, followed by S, N, and C, irrespective of chain length. The result revealed that the sequence of the stability is independent of coverage; however, linkers other than carbon can shift the optimum coverage considerably and allow further packing density. For all linkers apart from sulfur, structural properties, in particular, surface-linker-chain angles, saturate to a single value once n > 3. For sulfur, we identify three regimes, namely, n = 0–3, n = 5–7, and n = 9–11, each with its own characteristic adsorption structures. Where possible, our computational results are shown to be consistent with the available experimental data and show how the fundamental structural properties of modified Si surfaces can be controlled by the choice of linking group and chain length. Later we continue by examining the work function tuning of H:Si(111) over a range of 1.73 eV through adsorption of alkyl monolayers with general formula -[Xhead-group]-(CnH2n)-[Xtail-group], X = O(H), S(H), NH(2). The work function is practically converged at 4 carbons (8 for oxygen), for head-group functionalization. For tail-group functionalization and with both head- and tail-groups, there is an odd-even effect in the behavior of the work function, with peak-to-peak amplitudes of up to 1.7 eV in the oscillations. This behavior is explained through the orientation of the terminal-group's dipole. The shift in the work function is largest for NH2-linked and smallest for SH-linked chains and is rationalized in terms of interface dipoles. Our study reveals that the choice of the head- and/or tail-groups effectively changes the impact of the alkyl chain length on the work function tuning using self-assembled monolayers and this is an important advance in utilizing hybrid functionalized Si surfaces. Bringing together the understanding gained from studying single type functionalization of H:Si(111) with different alkyl chains and bearing in mind how to utilize head-group, tail-group or both as well as monolayer coverage, in the final part of this Thesis we study functionalized H:Si(111) with binary SAMs. Aiming at enhancing work function adjustment together with SAM stability and coverage we choose a range of terminations and linker-chains denoted as –X–(Alkyl) with X = CH3, O(H), S(H), NH(2) and investigate the stability and work function of various binary components grafted onto H:Si(111) surface. Using binary functionalization with -[NH(2)/O(H)/S(H)]-[Hexyl/Dodecyl] we show that work function can be tuned within the interval of 3.65-4.94 eV and furthermore, enhance the SAM’s stability. Although direct Si-C grafted SAMs are less favourable compared to their counterparts with O, N or S linkage, regardless of the ratio, binary functionalized alkyl monolayers with X-alkyl (X = NH, O) is always more stable than single type alkyl functionalization with the same coverage. Our results indicate that it is possible to go beyond the optimum coverage of pure alkyl functionalized SAMs (50%) by adding a linker with the correct choice of the linker. This is very important since dense packed monolayers have fewer defects and deliver higher efficiency. Our results indicate that binary anchoring can modify the charge injection and therefore bond stability while preserving the interface electronic structure.

Surface Insulation Resistance Degradation and Electrochemical Migration on Printed Circuit Boards

Widespread adoption of lead-free materials and processing for printed circuit board (PCB) assembly has raised reliability concerns regarding surface insulation resistance (SIR) degradation and electrochemical migration (ECM). As PCB conductor spacings decrease, electronic products become more susceptible to these failures mechanisms, especially in the presence of surface contamination and flux residues which might remain after no-clean processing. Moreover, the probability of failure due to SIR degradation and ECM is affected by the interaction between physical factors (such as temperature, relative humidity, electric field) and chemical factors (such as solder alloy, substrate material, no-clean processing). Current industry standards for assessing SIR reliability are designed to serve as short-term qualification tests, typically lasting 72 to 168 hours, and do not provide a prediction of reliability in long-term applications. The risk of electrochemical migration with lead-free assemblies has not been adequately investigated. Furthermore, the mechanism of electrochemical migration is not completely understood. For example, the role of path formation has not been discussed in previous studies. Another issue is that there are very few studies on development of rapid assessment methodologies for characterizing materials such as solder flux with respect to their potential for promoting ECM. In this dissertation, the following research accomplishments are described: 1). Long-term temp-humidity-bias (THB) testing over 8,000 hours assessing the reliability of printed circuit boards processed with a variety of lead-free solder pastes, solder pad finishes, and substrates. 2). Identification of silver migration from Sn3.5Ag and Sn3.0Ag0.5Cu lead-free solder, which is a completely new finding compared with previous research. 3). Established the role of path formation as a step in the ECM process, and provided clarification of the sequence of individual steps in the mechanism of ECM: path formation, electrodeposition, ion transport, electrodeposition, and filament formation. 4). Developed appropriate accelerated testing conditions for assessing the no-clean processed PCBs' susceptibility to ECM: a). Conductor spacings in test structures should be reduced in order to reflect the trend of higher density electronics and the effect of path formation, independent of electric field, on the time-to-failure. b). THB testing temperatures should be modified according to the material present on the PCB, since testing at 85oC can cause the evaporation of weak organic acids (WOAs) in the flux residues, leading one to underestimate the risk of ECM. 5). Correlated temp-humidity-bias testing with ion chromatography analysis and potentiostat measurement to develop an efficient and effective assessment methodology to characterize the effect of no-clean processing on ECM.