Multi-scale modelling of Suzuki segregation in gamma' precipitates in Ni and Co-base superalloys

The high temperature strength of alloys with (gamma +gamma') microstructure is primarily due to the resistance of the ordered precipitate to cutting by matrix dislocations. Such shearing requires higher stresses since it involves the creation of a planar fault. Planar fault energy is known to be dependent on composition. This implies that the composition on the fault may be different from that in the bulk for energetic reasons. Such segregation (or desegregation) of specific alloying elements to the fault may result in Suzuki strengthening which has not been explored extensively in these systems. In this work, segregation (or desegregation) of alloying elements to planar faults was studied computationally in Ni-3(Al, Ti) and Co-3(W, Al) type gamma' precipitates. The composition dependence of APB energy and heat of mixing were evaluated from first principle electronic structure calculations. A phase field model incorporating the first principles results, was used to simulate the motion of an extended superdislocation under stress concurrently with composition evolution. Results reveal that in both systems, significant (de) segregation occurs on equilibration. On application of stress, solutes were dragged along with the APB in some cases. Additionally, it was also noted the velocity of the superdislocation under an applied stress is strongly dependent on atomic mobility (i. e. diffusivity).

DESIGN STRATEGIES FOR CIRCULAR ECONOMY

A Circular Economy (CE) values material, technical or biological, as nutrient. CE thinking seeks to accelerate the conversion of technical nutrient cycles along the lines of biological nutrient cycles by re-designing systems till the scale of the economy. Though the notion of products being technical nutrient exists, its situation as an outcome of design intent is not contextually made. One objective of this article is to situate design and nutrient cycles of the earth system as and within natural cycles. This situation emphasizes the mechanism by which design affects nutrient availability to vital earth systems and draws attention to the functions that nutrients afford and serve by default before being embodied in products by human intent. The first principle of CE seeks to eliminate waste and re-purpose nutrients with minimal energy. Towards this, the historic trend of perceiving waste is drawn and Gestalts identified to arrive at the concept of tenancy and inform design. Tenancy is defined as the duration for which the nutrient embodied serves some purpose. Identifying the 6R scenarios as nutrient re-purposing functions, corresponding design strategies are stated.

Length scales in glass-forming liquids and related systems: a review

The central problem in the study of glass-forming liquids and other glassy systems is the understanding of the complex structural relaxation and rapid growth of relaxation times seen on approaching the glass transition. A central conceptual question is whether one can identify one or more growing length scale(s) associated with this behavior. Given the diversity of molecular glass-formers and a vast body of experimental, computational and theoretical work addressing glassy behavior, a number of ideas and observations pertaining to growing length scales have been presented over the past few decades, but there is as yet no consensus view on this question. In this review, we will summarize the salient results and the state of our understanding of length scales associated with dynamical slow down. After a review of slow dynamics and the glass transition, pertinent theories of the glass transition will be summarized and a survey of ideas relating to length scales in glassy systems will be presented. A number of studies have focused on the emergence of preferred packing arrangements and discussed their role in glassy dynamics. More recently, a central object of attention has been the study of spatially correlated, heterogeneous dynamics and the associated length scale, studied in computer simulations and theoretical analysis such as inhomogeneous mode coupling theory. A number of static length scales have been proposed and studied recently, such as the mosaic length scale discussed in the random first-order transition theory and the related point-to-set correlation length. We will discuss these, elaborating on key results, along with a critical appraisal of the state of the art. Finally we will discuss length scales in driven soft matter, granular fluids and amorphous solids, and give a brief description of length scales in aging systems. Possible relations of these length scales with those in glass-forming liquids will be discussed.

I. Rigid body penetration into brittle materials. II. Phase change effect on shock wave propagation

Part I.

We have developed a technique for measuring the depth time history of rigid body penetration into brittle materials (hard rocks and concretes) under a deceleration of ~ 10⁵ g. The technique includes bar-coded projectile, sabot-projectile separation, detection and recording systems. Because the technique can give very dense data on penetration depth time history, penetration velocity can be deduced. Error analysis shows that the technique has a small intrinsic error of ~ 3-4 % in time during penetration, and 0.3 to 0.7 mm in penetration depth. A series of 4140 steel projectile penetration into G-mixture mortar targets have been conducted using the Caltech 40 mm gas/ powder gun in the velocity range of 100 to 500 m/s.

We report, for the first time, the whole depth-time history of rigid body penetration into brittle materials (the G-mixture mortar) under 10⁵ g deceleration. Based on the experimental results, including penetration depth time history, damage of recovered target and projectile materials and theoretical analysis, we find:

1. Target materials are damaged via compacting in the region in front of a projectile and via brittle radial and lateral crack propagation in the region surrounding the penetration path. The results suggest that expected cracks in front of penetrators may be stopped by a comminuted region that is induced by wave propagation. Aggregate erosion on the projectile lateral surface is < 20% of the final penetration depth. This result suggests that the effect of lateral friction on the penetration process can be ignored.

2. Final penetration depth, P_max, is linearly scaled with initial projectile energy per unit cross-section area, e_s , when targets are intact after impact. Based on the experimental data on the mortar targets, the relation is P_max(mm) 1.15e_s (J/mm² ) + 16.39.

3. Estimation of the energy needed to create an unit penetration volume suggests that the average pressure acting on the target material during penetration is ~ 10 to 20 times higher than the unconfined strength of target materials under quasi-static loading, and 3 to 4 times higher than the possible highest pressure due to friction and material strength and its rate dependence. In addition, the experimental data show that the interaction between cracks and the target free surface significantly affects the penetration process.

4. Based on the fact that the penetration duration, t_max, increases slowly with e_s and does not depend on projectile radius approximately, the dependence of t_max on projectile length is suggested to be described by t_max(μs) = 2.08e_s (J/mm² + 349.0 x m/(πR²), in which m is the projectile mass in grams and R is the projectile radius in mm. The prediction from this relation is in reasonable agreement with the experimental data for different projectile lengths.

5. Deduced penetration velocity time histories suggest that whole penetration history is divided into three stages: (1) An initial stage in which the projectile velocity change is small due to very small contact area between the projectile and target materials; (2) A steady penetration stage in which projectile velocity continues to decrease smoothly; (3) A penetration stop stage in which projectile deceleration jumps up when velocities are close to a critical value of ~ 35 m/s.

6. Deduced averaged deceleration, a, in the steady penetration stage for projectiles with same dimensions is found to be a(g) = 192.4v + 1.89 x 10⁴, where v is initial projectile velocity in m/s. The average pressure acting on target materials during penetration is estimated to be very comparable to shock wave pressure.

7. A similarity of penetration process is found to be described by a relation between normalized penetration depth, P/P_max, and normalized penetration time, t/t_max, as P/P_max = f(t/t_max, where f is a function of t/t_max. After f(t/t_max is determined using experimental data for projectiles with 150 mm length, the penetration depth time history for projectiles with 100 mm length predicted by this relation is in good agreement with experimental data. This similarity also predicts that average deceleration increases with decreasing projectile length, that is verified by the experimental data.

8. Based on the penetration process analysis and the present data, a first principle model for rigid body penetration is suggested. The model incorporates the models for contact area between projectile and target materials, friction coefficient, penetration stop criterion, and normal stress on the projectile surface. The most important assumptions used in the model are: (1) The penetration process can be treated as a series of impact events, therefore, pressure normal to projectile surface is estimated using the Hugoniot relation of target material; (2) The necessary condition for penetration is that the pressure acting on target materials is not lower than the Hugoniot elastic limit; (3) The friction force on projectile lateral surface can be ignored due to cavitation during penetration. All the parameters involved in the model are determined based on independent experimental data. The penetration depth time histories predicted from the model are in good agreement with the experimental data.

9. Based on planar impact and previous quasi-static experimental data, the strain rate dependence of the mortar compressive strength is described by σ_f/σ⁰_f = exp(0.0905(log(έ/έ_0) ^1.14, in the strain rate range of 10^-7/s to 10³/s (σ⁰_f and έ are reference compressive strength and strain rate, respectively). The non-dispersive Hugoniot elastic wave in the G-mixture has an amplitude of ~ 0.14 GPa and a velocity of ~ 4.3 km/s.

Part II.

Stress wave profiles in vitreous GeO₂ were measured using piezoresistance gauges in the pressure range of 5 to 18 GPa under planar plate and spherical projectile impact. Experimental data show that the response of vitreous GeO₂ to planar shock loading can be divided into three stages: (1) A ramp elastic precursor has peak amplitude of 4 GPa and peak particle velocity of 333 m/s. Wave velocity decreases from initial longitudinal elastic wave velocity of 3.5 km/s to 2.9 km/s at 4 GPa; (2) A ramp wave with amplitude of 2.11 GPa follows the precursor when peak loading pressure is 8.4 GPa. Wave velocity drops to the value below bulk wave velocity in this stage; (3) A shock wave achieving final shock state forms when peak pressure is > 6 GPa. The Hugoniot relation is D = 0.917 + 1.711u (km/s) using present data and the data of Jackson and Ahrens [1979] when shock wave pressure is between 6 and 40 GPa for ρ₀ = 3.655 gj cm³ . Based on the present data, the phase change from 4-fold to 6-fold coordination of Ge⁺⁴ with O^-2 in vitreous GeO₂ occurs in the pressure range of 4 to 15 ± 1 GPa under planar shock loading. Comparison of the shock loading data for fused SiO₂ to that on vitreous GeO₂ demonstrates that transformation to the rutile structure in both media are similar. The Hugoniots of vitreous GeO₂ and fused SiO₂ are found to coincide approximately if pressure in fused SiO₂ is scaled by the ratio of fused SiO₂to vitreous GeO₂ density. This result, as well as the same structure, provides the basis for considering vitreous Ge0₂ as an analogous material to fused SiO₂ under shock loading. Experimental results from the spherical projectile impact demonstrate: (1) The supported elastic shock in fused SiO₂ decays less rapidly than a linear elastic wave when elastic wave stress amplitude is higher than 4 GPa. The supported elastic shock in vitreous GeO₂ decays faster than a linear elastic wave; (2) In vitreous GeO₂ , unsupported shock waves decays with peak pressure in the phase transition range (4-15 GPa) with propagation distance, x, as α 1/x^-3.35 , close to the prediction of Chen et al. [1998]. Based on a simple analysis on spherical wave propagation, we find that the different decay rates of a spherical elastic wave in fused SiO₂ and vitreous GeO₂ is predictable on the base of the compressibility variation with stress under one-dimensional strain condition in the two materials.

Band structure and transport properties of simple cubic tellurium-base alloys

The electrical transport properties and lattice spacings of simple cubic Te-Au, Te-Au-Fe, and Te-Au-Mn alloys, prepared by rapid quenching from the liquid state, hove been measured and correlated with a proposed bond structure. The variations of superconducting transition temperature, absolute thermoelectric power, and lattice spacing with Te concentration all showed related anomalies in the binary Te-Au alloys. The unusual behavior of these properties has been interpreted by using nearly free electron theory to predict the effect of the second Brillouin zone boundary on the area of the Fermi surface, and the electronic density of states. The behavior of the superconducting transition temperature and the lattice parameter as Fe and Mn ore added further supports the proposed interpretation as well as providing information on the existence of localized magnetic states in the ternary alloys. In addition, it was found that a very distinct bond structure effect on the transition temperatures of the Te-Au-Fe alloys could be identified.

The analytical description and experiments of the optical bottle generated by an axicon and a lens

Based on the interferential theory, we deduce a new type of analytic expression suitable for describing the evolutions of the optical bottle beam generated from the axicon-lens optical system illuminated by the Gaussian beam for the first time. The theory does not use much approximation in the process of mathematical analysis and can better illustrate the optical bottle beam evolutions at any positions. With the derived expression, the three-dimensional (3D) longitudinal and transverse intensity profiles of the optical bottle beam are simulated numerically. The numerical calculations have been confirmed by the experimental results.

Enhanced electroluminescence intensity of InGaN/GaN multi-quantum-wells based on Mg-doped GaN annealed in O-2

InGaN/GaN multi-quantum-well blue (461 +/- 4 nm) light emitting diodes with higher electroluminescence intensity are obtained by postgrowth thermal annealing at 720 C in O-2-ambient. Based on our first-principle total-energy calculations, we conclude that besides dissociating the Mg-H complex by forming H2O, annealing in O-2 has another positive effect on the activation of acceptor Mg in GaN. Mg can be further activated by the formation of an impurity band above the valence band maximum of host GaN from the passivated Mg-Ga-O-N complex. Our calculated ionization energy for acceptor Mg in the passivated system is about 30 meV shallower than that in pure GaN, in good agreement with previous experimental measurement. Our model can explain that the enhanced electroluminescence intensity of InGaN/GaN MQWs based on Mg-doped p-type GaN is due to a decrease in the ionization energy of Mg acceptor with the presence of oxygen. (C) 2008 American Institute of Physics.

The bipolar doping of ZnS via native defects and external dopants

By employing first-principle total-energy calculations, a systematic study of the dopability of ZnS to be both n- and p-types compared with that of ZnO is carried out. We find that all the attempted acceptor dopants, group V substituting on the S lattice site and group I and IB on the Zn sites in ZnS, have lower ionization energies than the corresponding ones in ZnO. This can be accounted for by the fact that ZnS has relative higher valence band maximum than ZnO. Native ZnS is weak p-type under S-rich condition, as the abundant acceptor V-Zn has rather large ionization energy. Self-compensations by the formation of interstitial donors in group I and IB-doped p-type ZnS can be avoided when sample is prepared under S-rich condition. In terms of ionization energies, Li-Zn and N-S are the preferred acceptors in ZnS. Native n- type doping of ZnS is limited by the spontaneous formation of intrinsic V-Zn(2-); high efficient n-type doping with dopants is harder to achieve than in ZnO because of the readiness of forming native compensating centers and higher ionization energy of donors in ZnS. (C) 2009 American Institute of Physics. [DOI 10.1063/1.3103585]

CALCULATED ELECTRONIC AND MAGNETIC-STRUCTURES OF THE NEW TERNARY RARE-EARTH IRON NITRIDE ND2FE17N3

The electronic and magnetic structures of Nd2Fe17 and Nd2Fe17N3 have been calculated using the first-principle, spin-polarized orthogonalized linear combination of atomic orbitals method. Comparative studies of the two materials reveal important effects of the nitrogen atoms (at 9e site) on the electronic and magnetic structures. Results are presented for the total density of states, site-projected partial density of states and the spin magnetic moments on four nonequivalent Fe sites. The highest magnetic moments are found to be located on the 6c site for Nd2Fe17 and on the 9d site for Nd2Fe17N3, in agreement with the neutron and Mossbauer experiments. The variation trends of the magnetic moments on different Fe sites are discussed in terms of the separation between Fe and N atoms. Compared with Nd2Fe17, an increase in the exchange splitting of the Fe d band is found in Nd2Fe17N3, which accounts for its higher Curie temperature as observed in experiments. The calculated results show that the nitrogen atoms are charge acceptors in these compounds.

CALCULATION OF DIELECTRIC-CONSTANTS FOR SEMICONDUCTORS - AN APPLICATION OF ABINITIO PSEUDOPOTENTIAL METHOD

We have used ab initio pseudopotential method to generate basis wavefunctions and eigen energies to carry out first principle calculations of the static macroscopic dielectric constant for GaAs and GaP. The resulted converged random phase approximation (RPA) value is 12.55 and 10.71, in excellent agreement to the experimental value of 12.4 and 10.86, respectively. The inclusion of the exchange correlation contribution makes the calculated result slightly worsen. A convergence test with respect to the number of k points in Brillouin zone (BZ) integration was carried out. Sixty irreducible BZ k points were used to achieve the converged results. Integration with only 10 special k points increased the RPA value by 15%.

Photoluminescence of GaInP under high pressure

Photoluminescence of GaInP under hydrostatic pressure is investigated. The Gamma valley of disordered GaInP shifts sublinearly upwards with respect to the top of the valence band with increasing pressure and this sublinearity is caused by the nonlinear relationship between lattice constant and hydrostatic pressure. The Gamma valleys of ordered GaInP rise more slowly than that of the disordered one and the relationship between the band gap and the pressure can not be explained in the same way. Taking into account the interactions between the Gamma valley and the folded L valleys, as well as, the X valleys, the experimental pressure dependences of the band gap of ordered GaInP epilayers are calculated and fitted quite well using first order perturbation theory. The results indicate that simultaneous ordering along [111] and [100] directions can occur in ordered GaInP. (C) 1996 American Institute of Physics.

ELECTROCHEMICAL AND IN-SITU SPECTROELECTROCHEMICAL STUDY OF (TETRAPHENYLPORPHINATO)MANGANESE(III) PERCHLORATE IN THE PRESENCE OF CHLORIDE ANIONS .2. OXIDATION PROCESSES

An investigation of electrode oxidation processes of (tetra-phenylporphinato) manganese (III) Perchlorate, (TPS)Mn(III)ClO4, was carried out during the titration of chloride anions by conventional cyclic voltammetry, thin-layer cyclic voltammetry and spectroelectrochemistry. It was demonstrated that in the presence of one equivalent amount of Cl-, the first one electron oxidation reaction corresponds to the Mn(III)I cation radical oxidation, and the second one electron oxidation corresponds to the cation radical/dication generation followed by an iso-porphyrin formation reaction, however in the presence of two equivalent amount of Cl-, the first one electron oxidation of Mn(III) gives Mn(IV) porphyrin and the second one electron oxidation generates cation radicals of Mn(IV) followed by an iso-porphyrin formation reactions. Mechanisms of these redox processes are postulated.

Genetic mapping of size-related quantitative trait loci (QTL) in the bay scallop (Argopecten irradians) using AFLP and microsatellite markers

We constructed genetic linkage maps for the bay scallop Argopecten irradians using AFLP and microsatellite markers and conducted composite interval mapping (CIM) of body-size-related traits. Three hundred seventeen AFLP and 10 microsatellite markers were used for map construction. The female parent map contained 120 markers in 15 linkage groups, spanning 479.6 cM with an average interval of 7.0 cM. The male parent map had 190 markers in 17 linkage groups, covering 883.8 cM at 7.2 cM per marker. The observed coverage was 70.4% for the female and 81.1% for the male map. Markers that were distorted toward the same direction were closely linked to each other on the genetic maps, suggesting the presence of genes important for survival. Six size-related traits, shell length, shell height, shell width, total weight, soft tissue weight, and shell weight, were measured for QTL mapping. The size data were significantly correlated with each other. We subjected the data, log transformed firstly, to a principle component analysis and use the first principle component for QTL mapping. CIM analysis revealed one significant QTL (LOD=2.69, 1000 permutation, P<0.05) in linkage group 3 on the female parent map. The mapping of size-related QTL in this study raises the possibility of improving the growth of bay scallops through marker-assisted selection. (c) 2007 Published by Elsevier B.V.

通过第一性原理计算研究内地核的结构与成分

The determination of the composition and structure of the Earth’s inner core has long been the major subject in the study of the Earth’s deep interior. It’s widely believed that the Earth’s core is formed by iron with a fraction of nickel. However, light elements must exist in the inner core because the earth core is less dense than pure iron-nickel alloy (~2-3% in the solid inner core and ~6-7% in the liquid outer core). The questions are what and how much light element is there in the iron-nickel alloy. Besides the composition, the crystal structure of the iron with or without light element is also not well known. According to the seismological observations, the sound waves propagate 3-4% faster along the spin axis than in the equatorial plane. That means the inner core is anisotropic. The densest structure of iron-nickel alloy should be h.c.p structure under the very high pressures. However, the h,c,p structure does not propagate waves anisotropic ally. Then what is the structure of the iron-nickel alloy or the iron-nickle-light element alloy. In this study, we tried to predict the composition and the structure of the inner core through ab initio calculation of the Gibbs free energy, which is a function of internal energy, density and entropy. We conclude that the h.c.p structure is more stable than the b.c.c structure under high pressure and 0 K, but with the increase of temperature, the free energy of the b.c.c structure is decreasing much faster than the h.c.p structure caused by the vibration of the atomics, so the b.c.c structure is more stable at high temperatures. With the addition of light elements (S or Si or both), the free energy of b.c.c. decreases even faster, about 3at% of Si not only explains why the inner core is about 2-3 % lighter than the iron-nickle alloy, but also reasons why the inner core is anisotropic, since the b.c.c. structure becomes more stable than the h.c.p structure at 5500-6000K and b.c.c. is anisotropic in propagating seismic waves. Therefore, we infer that the inner core of the earth is formed by b.c.c iron and a fraction of nickel plus ~3at.% Si, with a temperature higher than 5500K, which is consistent with the studies from other approaches.

分子水平上研究地质流体的物理化学性质

Geological fluids exist in every geosphere of the Earth and play important roles in many processes of material transformations, energetic interchanges and geochemical interactions. To study the physicochemical properties and geochemical behaviors of geological fluids turn Girt to be one of the challenging issues in geosciences. Compared with conventional approaches of experiments and semi-theoretical modeling, computer simulation on molecular level shows its advantages on quantitative predictions of the physicochemical properties of geological fluids under extreme conditions and emerges as a promising approach to find the characteristics of geological fluids and their interactions in different geospheres of the Earth interior.This dissertation systematically discusses the physicochemical properties of typical geological fluids with state-of-the-art computer simulation techniques. The main results can be summarized as follows: (1) The experimental phase behaviors of the systems CH4-C2H6 and. CO2 have been successfully reproduced with Monte Carlo simulations. (2) Through comprehensive isothermal-isobaric molecular dynamics simulations, the PVT data of water hia^e been extended beyond experimental range to about 2000 K and 20 GPa and an improved equation of state for water has been established. (3) Based on extensive computer simulations, am optimized molecular potential for carbon dioxide have been proposed, this model is expected to predict different properties of carbon dioxide (volumetric properties, phase equilibria, heat of vaporization, structural and dynamical properties) with improved accuracies. (4) On the basis of the above researches of the end-members, a set of parameters for unlike interactions has been proposed by non-linear fitting to the ab initio potential surface of CO2-H2O and is superior to the common used mixing rule and the results of prior workers vs/Ith remarkable accuracies, then a number of simulations of the mixture have been carried out to generate data under high temperatures and pressures as an important complement to the limited experiments. (5) With molecular dynamics simulations, various structural, dynamical and thermodynamical properties of ionic solvations and associations have been oomprehensively analyzed, these results not only agree well with experimental data and first principle calculation results, but also reveal some new insights into the microscopic ionic solvation and association processes.