High temperature thermoelectric properties of Zr and Hf based transition metal dichalcogenides: A first principles study

We investigate the electronic and thermal transport properties of bulk MX2 compounds (M = Zr, Hf and X = S, Se) by first-principles calculations and semi-classical Boltzmann transport theory. The band structure shows the confinement of heavy and light bands along the out of plane and in-plane directions, respectively. This results in high electrical conductivity (sigma) and large thermopower leading to a high power factor (S-2 sigma) for moderate n-type doping. The phonon dispersion demonstrates low frequency flat acoustical modes, which results in low group velocities (v(g)). Consequently, lowering the lattice thermal conductivity (kappa(latt)) below 2 W/m K. Low kappa(latt) combined with high power factor results in ZT > 0.8 for all the bulk MX2 compounds at high temperature of 1200 K. In particular, the ZT(max) of HfSe2 exceeds 1 at 1400 K. Our results show that Hf/Zr based dichalcogenides are very promising for high temperature thermoelectric application. (C) 2015 AIP Publishing LLC.

Growth mechanism of the interdiffusion zone between platinum modified bond coats and single crystal superalloys

Pt-modified beta-NiAl bond coats are applied over the superalloys for oxidation protection in jet engine applications. However, as shown in this study, it also enhances the growth of the interdiffusion zone developed between the bond coat and the superalloy along with brittle precipitates. Location of the Kirkendall plane indicates that a precipitate free sublayer grows from the bond coat, whereas another sublayer grows from the superalloy containing very high volume fraction of precipitates. With increasing Pt content, thickness of both the sublayers increases because of an increase in diffusion rates of the components. Quantitative electron probe microanalysis indicates high concentration of refractory components in the precipitates. Transmission electron microscopy shows that Rene N5 superalloy produces TCP phases mu and P, whereas CMSX-4 superalloy produces mu and sigma in the interdiffusion zone. With increasing Pt content in the bond coat, the average size of the precipitates decreases when coupled with Rene N5. Precipitates become much finer when the same bond coats are coupled with CMSX-4. (C) 2015 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Growth mechanism of the interdiffusion zone between platinum modified bond coats and single crystal superalloys

Pt-modified beta-NiAl bond coats are applied over the superalloys for oxidation protection in jet engine applications. However, as shown in this study, it also enhances the growth of the interdiffusion zone developed between the bond coat and the superalloy along with brittle precipitates. Location of the Kirkendall plane indicates that a precipitate free sublayer grows from the bond coat, whereas another sublayer grows from the superalloy containing very high volume fraction of precipitates. With increasing Pt content, thickness of both the sublayers increases because of an increase in diffusion rates of the components. Quantitative electron probe microanalysis indicates high concentration of refractory components in the precipitates. Transmission electron microscopy shows that Rene N5 superalloy produces TCP phases mu and P, whereas CMSX-4 superalloy produces mu and sigma in the interdiffusion zone. With increasing Pt content in the bond coat, the average size of the precipitates decreases when coupled with Rene N5. Precipitates become much finer when the same bond coats are coupled with CMSX-4. (C) 2015 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Localized Charge Carrier Transport Properties of Zn1-x Ni (x) O/NiO Two-Phase Composites

We report the localized charge carrier transport of two-phase composite Zn1-x Ni (x) O/NiO (0 a parts per thousand currency sign x a parts per thousand currency sign 1) using the temperature dependence of ac-resistivity rho (ac)(T) across the N,el temperature T (N) (= 523 K) of nickel oxide. Our results provide strong evidence to the variable range hopping of charge carriers between the localized states through a mechanism involving spin-dependent activation energies. The temperature variation of carrier hopping energy epsilon (h)(T) and nearest-neighbor exchange-coupling parameter J (ij)(T) evaluated from the small poleron model exhibits a well-defined anomaly across T (N). For all the composite systems, the average exchange-coupling parameter (J (ij))(AVG) nearly equals to 70 meV which is slightly greater than the 60-meV exciton binding energy of pure zinc oxide. The magnitudes of epsilon (h) (similar to 0.17 eV) and J (ij) (similar to 11 meV) of pure NiO synthesized under oxygen-rich conditions are consistent with the previously reported theoretical estimation based on Green's function analysis. A systematic correlation between the oxygen stoichiometry and, epsilon (h)(T) and J (ij)(T) is discussed.

Structure-rheology relationship in a sheared lamellar fluid

The structure-rheology relationship in the shear alignment of a lamellar fluid is studied using a mesoscale model which provides access to the lamellar configurations and the rheology. Based on the equations and free energy functional, the complete set of dimensionless groups that characterize the system are the Reynolds number (rho gamma L-2/mu), the Schmidt number (mu/rho D), the Ericksen number (mu(gamma)/B), the interface sharpness parameter r, the ratio of the viscosities of the hydrophilic and hydrophobic parts mu(r), and the ratio of the system size and layer spacing (L/lambda). Here, rho and mu are the fluid density and average viscosity, (gamma) over dot is the applied strain rate, D is the coefficient of diffusion, B is the compression modulus, mu(r) is the maximum difference in the viscosity of the hydrophilic and hydrophobic parts divided by the average viscosity, and L is the system size in the cross-stream direction. The lattice Boltzmann method is used to solve the concentration and momentum equations for a two dimensional system of moderate size (L/lambda = 32) and for a low Reynolds number, and the other parameters are systematically varied to examine the qualitative features of the structure and viscosity evolution in different regimes. At low Schmidt numbers where mass diffusion is faster than momentum diffusion, there is fast local formation of randomly aligned domains with ``grain boundaries,'' which are rotated by the shear flow to align along the extensional axis as time increases. This configuration offers a high resistance to flow, and the layers do not align in the flow direction even after 1000 strain units, resulting in a viscosity higher than that for an aligned lamellar phase. At high Schmidt numbers where momentum diffusion is fast, the shear flow disrupts layers before they are fully formed by diffusion, and alignment takes place by the breakage and reformation of layers by shear, resulting in defects (edge dislocations) embedded in a background of nearly aligned layers. At high Ericksen number where the viscous forces are large compared to the restoring forces due to layer compression and bending, shear tends to homogenize the concentration field, and the viscosity decreases significantly. At very high Ericksen number, shear even disrupts the layering of the lamellar phase. At low Ericksen number, shear results in the formation of well aligned layers with edge dislocations. However, these edge dislocations take a long time to anneal; the relatively small misalignment due to the defects results in a large increase in viscosity due to high layer stiffness and due to shear localization, because the layers between defects get pinned and move as a plug with no shear. An increase in the viscosity contrast between the hydrophilic and hydrophobic parts does not alter the structural characteristics during alignment. However, there is a significant increase in the viscosity, due to pinning of the layers between defects, which results in a plug flow between defects and a localization of the shear to a part of the domain.

Status of the 98-125 GeV Higgs bosons scenario with updated LHC-8 data

In the context of the minimal supersymmetric standard model (MSSM), we discuss the possibility of the lightest Higgs boson with mass M-h = 98 GeV to be consistent with the 2.3 sigma excess observed at the LEP in the decay mode e(+)e(-) -> Zh, with h -> b (b) over bar. In the same region of the MSSM parameter space, the heavier Higgs boson (H) with mass M-H similar to 125 GeV is required to be consistent with the latest data on Higgs coupling measurements at the end of the 7 + 8 TeV LHC run with 25 fb(-1) of data. While scanning the MSSM parameter space, we impose constraints coming from flavor physics, relic density of the cold dark matter as well as direct dark matter searches. We study the possibility of observing this light Higgs boson in vector boson fusion process and associated production with W/Z-boson at the high luminosity (3000 fb(-1)) run of the 14 TeV LHC. Our analysis shows that this scenario can hardly be ruled out even at the high luminosity run of the LHC. However, the precise measurement of the Higgs signal strength ratios can play a major role to distinguish this scenario from the canonical MSSM one.

Electron-hole asymmetry in the electron-phonon coupling in top-gated phosphorene transistor

Using in situ Raman scattering from phosphorene channel in an electrochemically top-gated field effect transistor, we show that phonons with A(g) symmetry depend much more strongly on concentration of electrons than that of holes, wheras phonons with B-g symmetry are insensitive to doping. With first-principles theoretical analysis, we show that the observed electon-hole asymmetry arises from the radically different constitution of its conduction and valence bands involving pi and sigma bonding states respectively, whose symmetry permits coupling with only the phonons that preserve the lattice symmetry. Thus, Raman spectroscopy is a non-invasive tool for measuring electron concentration in phosphorene-based nanoelectronic devices.

FAST AND ACCURATE BILATERAL FILTERING USING GAUSS-POLYNOMIAL DECOMPOSITION

The bilateral filter is a versatile non-linear filter that has found diverse applications in image processing, computer vision, computer graphics, and computational photography. A common form of the filter is the Gaussian bilateral filter in which both the spatial and range kernels are Gaussian. A direct implementation of this filter requires O(sigma(2)) operations per pixel, where sigma is the standard deviation of the spatial Gaussian. In this paper, we propose an accurate approximation algorithm that can cut down the computational complexity to O(1) per pixel for any arbitrary sigma (constant-time implementation). This is based on the observation that the range kernel operates via the translations of a fixed Gaussian over the range space, and that these translated Gaussians can be accurately approximated using the so-called Gauss-polynomials. The overall algorithm emerging from this approximation involves a series of spatial Gaussian filtering, which can be efficiently implemented (in parallel) using separability and recursion. We present some preliminary results to demonstrate that the proposed algorithm compares favorably with some of the existing fast algorithms in terms of speed and accuracy.

Evaluation of the Effectiveness of Representative Methods for Determining Young's Modulus and Hardness from Instrumented Indentation Data

The effectiveness of Oliver & Pharr's (O&P's) method, Cheng & Cheng's (C&C's) method, and a new method developed by our group for estimating Young's modulus and hardness based on instrumented indentation was evaluated for the case of yield stress to reduced Young's modulus ratio (sigma(y)/E-r) >= 4.55 x 10(-4) and hardening coefficient (n) <= 0.45. Dimensional theorem and finite element simulations were applied to produce reference results for this purpose. Both O&P's and C&C's methods overestimated the Young's modulus under some conditions, whereas the error can be controlled within +/- 16% if the formulation was modified with appropriate correction functions. Similar modification was not introduced to our method for determining Young's modulus, while the maximum error of results was around +/- 13%. The errors of hardness values obtained from all the three methods could be even larger and were irreducible with any correction scheme. It is therefore suggested that when hardness values of different materials are concerned, relative comparison of the data obtained from a single standard measurement technique would be more practically useful. It is noted that the ranges of error derived from the analysis could be different if different ranges of material parameters sigma(y)/E-r and n are considered.

Strain gradient theory with couple stress for crystalline solids

A new phenomenological strain gradient theory for crystalline solid is proposed. It fits within the framework of general couple stress theory and involves a single material length scale Ics. In the present theory three rotational degrees of freedom omega (i) are introduced, which denote part of the material angular displacement theta (i) and are induced accompanying the plastic deformation. omega (i) has no direct dependence upon u(i) while theta = (1 /2) curl u. The strain energy density omega is assumed to consist of two parts: one is a function of the strain tensor epsilon (ij) and the curvature tensor chi (ij), where chi (ij) = omega (i,j); the other is a function of the relative rotation tensor alpha (ij). alpha (ij) = e(ijk) (omega (k) - theta (k)) plays the role of elastic rotation reason The anti-symmetric part of Cauchy stress tau (ij) is only the function of alpha (ij) and alpha (ij) has no effect on the symmetric part of Cauchy stress sigma (ij) and the couple stress m(ij). A minimum potential principle is developed for the strain gradient deformation theory. In the limit of vanishing l(cs), it reduces to the conventional counterparts: J(2) deformation theory. Equilibrium equations, constitutive relations and boundary conditions are given in detail. For simplicity, the elastic relation between the anti-symmetric part of Cauchy stress tau (ij), and alpha (ij) is established and only one elastic constant exists between the two tensors. Combining the same hardening law as that used in previously by other groups, the present theory is used to investigate two typical examples, i.e., thin metallic wire torsion and ultra-thin metallic beam bend, the analytical results agree well with the experiment results. While considering the, stretching gradient, a new hardening law is presented and used to analyze the two typical problems. The flow theory version of the present theory is also given.

A crack perpendicular to and terminating at a bimaterial interface

Using dislocation simulation approach, the basic equation for a finite crack perpendicular to and terminating at a bimaterial interface is formulated. A novel expansion method is proposed for solving the problem. The complete solution to the problem, including the explicit formulae for the T stresses ahead of the crack tip and the stress intensity factors are presented. The stress held characteristics are analysed in detail. It is found that normal stresses sigma(x) and sigma(y) ahead of the crack tip, are characterised by Q fields if the crack is within a stiff material and the parameters \p(T)\ and \q(T)\ are very small, where Q is a generalised stress intensity factor for a crack normal to and terminating at the interface. If the crack is within a weak material, the normal stresses sigma(x) and sigma(y) are dominated by the Q field plus T stress.

Simulations of oil spread on a water surface

Numerical simulation of an oil slick spreading on still and wavy surfaces is described in this paper. The so-called sigma transformation is used to transform the time-varying physical domain into a fixed calculation domain for the water wave motions and, at the same time, the continuity equation is changed into an advection equation of wave elevation. This evolution equation is discretized by the forward time and central space scheme, and the momentum equations by the projection method. A damping zone is set up in front of the outlet boundary coupled with a Sommerfeld-Orlanski condition at that boundary to minimize the wave reflection. The equations for the oil slick are depth-averaged and coupled with the water motions when solving numerically. As examples, sinusoidal and solitary water waves, the oil spread on a smooth plane and on still and wavy water surfaces are calculated to examine the accuracy of simulating water waves by Navier-Stokes equations, the effect of damping zone on wave reflection and the precise structures of oil spread on waves.

Asymptotic analysis for a crack on interface of damaged materials

An asymptotic analysis for a crack lying on the interface of a damaged plastic material and a linear elastic material is presented in this paper. The present results show that the stress distributions along the crack tip are quite similar to those with HRR singularity field and the crack faces open obviously. Material constants n, mu and mo are varied to examine their effects on the resulting stress distributions and displacement distributions in the damaged plastic region. It is found that the stress components sigma(rr), sigma(theta theta), sigma(r theta) and sigma(e) are slightly affected by the changes of material constants n, mu and m(0), but the damaged plastic region are greatly disturbed by these material parameters.

The Stability of a Vertical Single-Walled Carbon Nanotube Under Its Own Weight

It has been reported recently that single carbon nanotubes were attached to AFM tips to act as nanotweezers. In order to investigate its stability, a vertical single-walled carbon nanotube (SWCNT) under its own weight is studied in this paper. The lower end of the carbon nanotube is clamped. Firstly the governing dimensionless numbers are derived by dimensional analysis. Then the theoretical analysis based on an elastic column model is carried out. Two ratios, I.e., the ratio of half wall thickness to radius (t=R) and the ratio of gravity to elastic resilience ($\rho$gR=E), and their influences on the ratio of critical length to radius are discussed. It is found that the relationship between the critical ratio of altitude to radius and ratio of half thickness to radius is approximately linear. As the dimensionless number $\rho$gR=E increases, the compressive force per unit length (weight) becomes larger, thus critical ratio of altitude to radius must become smaller to maintain stability. At last the critical length of SWCNT is calculated. The results of this paper will be helpful for the stability design of nanotweezers-like nanostructures.

Revelation of a Functional Dependence of the Sum of Two Uniaxial Strengths/Hardness on Elastic Work/Total Work of Indentation

Dimensional and finite element analyses were used to analyze the relationship between the mechanical properties and instrumented indentation response of materials. Results revealed the existence of a functional dependence of (engineering yield strength sigma(E,y) + engineering tensile strength sigma(E,b))/Oliver & Pharr hardness on the ratio of reversible elastic work to total work obtained from an indentation test. The relationship links up the Oliver & Pharr hardness with the material strengths, although the Oliver & Pharr hardness may deviate from the true hardness when sinking in or piling up occurs. The functional relationship can further be used to estimate the SUM sigma(E,y) + sigma(E,b) according to the data of an instrumented indentation test. The sigma(E,y) + sigma(E,b) value better reflects the strength of a material compared to the hardness value alone. The method was shown to be effective when applied to aluminum alloys. The relationship can further be used to estimate the fatigue limits, which are usually obtained from macroscopic fatigue tests in different modes.