Molecular dynamics simulation of nanoindentation of Fe3C and tetrahedral Fe4C

Study of nanomechanical response of iron carbides is important because presence of iron carbides greatly influences the performance and longevity of steel components. This work contributes to the literature by exploring nanoindentation of nanocrystalline Fe3C and tetrahedral-Fe4C using molecular dynamics simulation. The chemical interactions of iron and carbon were described through an analytical bond order inter-atomic potential (ABOP) energy function. The indentations were performed at an indentation speed of 50 m/sec and a repeat trial was performed at 5 m/sec. Load-displacement (P-h) curve for both these carbides showed residual indentation depth and maximum indentation depth (hf/hmax) ratio to be higher than 0.7 i.e. a circumstance where Oliver and Pharr method was not appropriate to be applied to evaluate the material properties. Alternate evaluation revealed Fe3C to be much harder than Fe4C. Gibbs free energy of formation and radial distribution function, coupled with state of the average local temperature and von Mises stresses indicate the formation of a new phase of iron-carbide. Formation of this newer phase was found to be due to deviatoric strain rather than the high temperature induced in the substrate during nanoindentation

Solid-liquid interfacial free energy out of equilibrium

The properties of the interface between solid and melt are key to solidification and melting, as the interfacial free energy introduces a kinetic barrier to phase transitions. This makes solidification happen below the melting temperature, in out-of-equilibrium conditions at which the interfacial free energy is ill defined. Here we draw a connection between the atomistic description of a diffuse solid-liquid interface and its thermodynamic characterization. This framework resolves the ambiguities in defining the solid-liquid interfacial free energy above and below the melting temperature. In addition, we introduce a simulation protocol that allows solid-liquid interfaces to be reversibly created and destroyed at conditions relevant for experiments. We directly evaluate the value of the interfacial free energy away from the melting point for a simple but realistic atomic potential, and find a more complex temperature dependence than the constant positive slope that has been generally assumed based on phenomenological considerations and that has been used to interpret experiments. This methodology could be easily extended to the study of other phase transitions, from condensation to precipitation. Our analysis can help reconcile the textbook picture of classical nucleation theory with the growing body of atomistic studies and mesoscale models of solidification.

SrTiO3(001)(2x1) reconstructions: First-principles calculations of surface energy and atomic structure compared with scanning tunneling microscopy images

(1x1) and (2x1) reconstructions of the (001) SrTiO3 surface were studied using the first-principles full-potential linear muffin-tin orbital method. Surface energies were calculated as a function of TiO2 chemical potential, oxygen partial pressure and temperature. The (1x1) unreconstructed surfaces were found to be energetically stable for many of the conditions considered. Under conditions of very low oxygen partial pressure the (2x1) Ti2O3 reconstruction [Martin R. Castell, Surf. Sci. 505, 1 (2002)] is stable. The question as to why STM images of the (1x1) surfaces have not been obtained was addressed by calculating charge densities for each surface. These suggest that the (2x1) reconstructions would be easier to image than the (1x1) surfaces. The possibility that the presence of oxygen vacancies would destabilise the (1x1) surfaces was also investigated. If the (1x1) surfaces are unstable then there exists the further possibility that the (2x1) DL-TiO2 reconstruction [Natasha Erdman Nature (London) 419, 55 (2002)] is stable in a TiO2-rich environment and for p(O2)>10(-18) atm.

Atomic excitation during recollision-free ultrafast multi-electron tunnel ionization

Modern intense ultrafast pulsed lasers generate an electric field of sufficient strength to permit tunnel ionization of the valence electrons in atoms(1). This process is usually treated as a rapid succession of isolated events, in which the states of the remaining electrons are neglected(2). Such electronic interactions are predicted to be weak, the exception being recollision excitation and ionization caused by linearly polarized radiation(3). In contrast, it has recently been suggested that intense field ionization may be accompanied by a two-stage 'shake-up' reaction(4). Here we report a unique combination of experimental techniques(5-8) that allows us to accurately measure the tunnel ionization probability for argon exposed to 50-fs laser pulses. Most significantly for the current study, this measurement is independent of the optical focal geometry(7,8), equivalent to a homogenous electric field. Furthermore, circularly polarized radiation negates recollision. The present measurements indicate that tunnel ionization results in simultaneous excitation of one or more remaining electrons through shake-up(9). From an atomic physics standpoint, it may be possible to induce ionization from specific states, and will influence the development of coherent attosecond extreme-ultraviolet-radiation sources(10). Such pulses have vital scientific and economic potential in areas such as high-resolution imaging of in vivo cells and nanoscale extreme-ultraviolet lithography.

Ab initio simulation of charged slabs at constant chemical potential

We present a practical scheme for performing ab initio supercell calculations of charged slabs at constant electron chemical potential mu, rather than at constant number of electrons N-e. To this end, we define the chemical potential relative to a plane (or "reference electrode") at a finite distance from the slab (the distance should reflect the particular geometry of the situation being modeled). To avoid a net charge in the supercell, and thus make possible a standard supercell calculation, we restore the electroneutrality of the periodically repeated unit by means of a compensating charge, whose contribution to the total energy and potential is subtracted afterwards. The "constant mu" mode enables one to perform supercell calculation on slabs, where the slab is kept at a fixed potential relative to the reference electrode. We expect this to be useful in modeling many experimental situations, especially in electro-chemistry. (C) 2001 American Institute of Physics.

A simple model of atomic interactions in noble metals based explicitly on electronic structure

A total energy tight-binding model with a basis of just one s state per atom is introduced. It is argued that this simplest of all tight-binding models provides a surprisingly good description of the structural stability and elastic constants of noble metals. By assuming inverse power scaling laws for the hopping integrals and the repulsive pair potential, it is shown that the density matrix in a perfect primitive crystal is independent of volume, and structural energy differences and equations of state are then derived analytically. The model is most likely to be of use when one wishes to consider explicitly and self-consistently the electronic and atomic structures of a generic metallic system, with the minium of computation expense. The relationship to the free-electron jellium model is described. The applicability of the model to other metals is also considered briefly.

Current-induced forces in atomic-scale conductors

We present a self-consistent tight-binding formalism to calculate the forces on individual atoms due to the flow of electrical current in atomic-scale conductors. Simultaneously with the forces, the method yields the local current density and the local potential in the presence of current flow, allowing a direct comparison between these quantities. The method is applicable to structures of arbitrary atomic geometry and can be used to model current-induced mechanical effects in realistic nanoscale junctions and wires. The formalism is implemented within a simple Is tight-binding model and is applied to two model structures; atomic chains and a nanoscale wire containing a vacancy.

The importance of hydrogen's potential-energy surface and the strength of the forming R-H bond in surface hydrogenation reactions

An understanding of surface hydrogenation reactivity is a prevailing issue in chemistry and vital to the rational design of future catalysts. In this density-functional theory study, we address hydrogenation reactivity by examining the reaction pathways for N+H -> NH and NH+H -> NH2 over the close-packed surfaces of the 4d transition metals from Zr-Pd. It is found that the minimum-energy reaction pathway is dictated by the ease with which H can relocate between hollow-site and top-site adsorption geometries. A transition state where H is close to a top site reduces the instability associated with bond sharing of metal atoms by H and N (NH) (bonding competition). However, if the energy difference between hollow-site and top-site adsorption energies (Delta E-H) is large this type of transition state is unfavorable. Thus we have determined that hydrogenation reactivity is primarily controlled by the potential-energy surface of H on the metal, which is approximated by Delta E-H, and that the strength of N (NH) chemisorption energy is of less importance. Delta E-H has also enabled us to make predictions regarding the structure sensitivity of these reactions. Furthermore, we have found that the degree of bonding competition at the transition state is responsible for the trend in reaction barriers (E-a) across the transition series. When this effect is quantified a very good linear correlation is found with E-a. In addition, we find that when considering a particular type of reaction pathway, a good linear correlation is found between the destabilizing effects of bonding competition at the transition state and the strength of the forming N-H (HN-H) bond. (c) 2006 American Institute of Physics.

Laser-Induced Fluorescence Measurements of Absolute Atomic Densities: Concepts and Limitations

The potential of laser-induced fluorescence spectroscopy of atoms is reviewed with emphasis on the determination of absolute densities. Examples of experiments with single-photon and two-photon excitation are presented. Calibration methods applicable with the different schemes are discussed. A new method is presented that has the potential to allow absolute measurement in plasmas of elevated pressure where collisional depletion of the excited state is present.

Potential and kinetic sputtering of alkanethiol self-assembled monolayers by impact of highly charged ions

Highly charged ions have been used to study the sputtering of positive molecular fragments from mercaptoundecanoic acid and dodecanethiol self-assembled monolayers on gold surfaces. The samples were bombarded with Arq+ (42n⁺, and Cn+1O2H_{2n + 1}⁺ from mercaptoundecanoic and H+, CnH_2n⁺, and Cn+1H_{2n + 3}⁺ from dodecanethiol. The proton yields were increased with larger charge state q of the highly charged ion (HCI) in both samples, scaling as qgamma, with gamma~5. The charge state dependence is discussed in terms of electron transfer to the HCI. The final yield of protons depends on molecular functional group characteristics, orientation on the surface, and reneutralization phenomena.

Melt-compounded nanocomposites of titanium dioxide atomic-layer-deposition-coated polyamide and polystyrene powders

Polyamide and polystyrene particles were coated with titanium dioxide films by atomic layer deposition (ALD) and then melt-compounded to form polymer nanocomposites. The rheological properties of the ALD-created nanocomposite materials were characterized with a melt flow indexer, a melt flow spiral mould, and a rotational rheometer. The results suggest that the melt flow properties of polyamide nanocomposites were markedly better than those of pure polyamide and polystyrene nanocomposites. Such behavior was shown to originate in an uncontrollable decrease in the polyamide molecular weight, likely affected by a high thin-film impurity content, as shown in gel permeation chromatography (GPC) and scanning electron microscope (SEM) equipped with an energy-dispersive spectrometer. Transmission electron microscope image showed that a thin film grew on both studied polymer particles, and that subsequent melt-compounding was successful, producing well dispersed ribbon-like titanium dioxide with the titanium dioxide filler content ranging from 0.06 to 1.12wt%. Even though we used nanofillers with a high aspect ratio, they had only a minor effect on the tensile and flexural properties of the polystyrene nanocomposites. The mechanical behavior of polyamide nanocomposites was more complex because of the molecular weight degradation. Our approach here to form polymeric nanocomposites is one way to tailor ceramic nanofillers and form homogenous polymer nanocomposites with minimal work-related risks in handling powder form nanofillers. However, further research is needed to gauge the commercial potential of ALD-created nanocomposite materials. Copyright (C) 2011 John Wiley & Sons, Ltd.

Density-potential mapping in time-dependent density-functional theory

The key questions of uniqueness and existence in time-dependent density-functional theory are usually formulated only for potentials and densities that are analytic in time. Simple examples, standard in quantum mechanics, lead, however, to nonanalyticities. We reformulate these questions in terms of a nonlinear Schroedinger equation with a potential that depends nonlocally on the wave function.

Characterization and parametrization in terms of atomic number of x-ray emission from K-shell filling during ion-surface interactions

When highly charged ions are incident on a surface, part of their potential energy is emitted as characteristic radiation. The energies and yields of these characteristic x rays have been measured for a series of elements at the Tokyo electron-beam ion trap. These data have been used to develop a simple model of the relaxation of the hollow atoms which are formed as the ion approaches the surface, as well as a set of semiempirical scaling laws, which allow for the ready calculation of the K-shell x-ray spectrum which would be produced by an arbitrary slow bare or hydrogenlike ion on a surface. These semiempirical scaling laws can be used to assess the merit of highly charged ion fluorescence x-ray generation in a wide range of applications.

The potential of electron beam radiation for simultaneous surface modification and bioresorption control of PLLA

Bioresorbable polymers have been widely investigated as materials exhibiting significant potential for successful application in the fields of tissue engineering and drug delivery. Further to the ability to control degradation, surface engineering of polymers has been highlighted as a key method central to their development. Previous work has demonstrated the ability of electron beam (e-beam) technology to control the degradation profiles and bioresorption of a number of commercially relevant bioresorbable polymers (poly-l-lactic acid (PLLA), Llactide/DL-lactide co-polymer (PLDL) and poly(lactic-co-glycolic acid (PLGA)). This work investigates the further potential of ebeam technology to impart added biofunctionality through the manipulation of polymer (PLLA) surface properties. PLLA samples were subjected to e-beam treatments in air, with varying beam energies and doses. Surface characterization was then performed using contact angle analysis, X-ray photoelectron spectroscopy (XPS), Raman spectroscopy, and atomic force microscopy. Results demonstrated a significant increase in surface wettability post e-beam treatment. In correlation with this, XPS data showed the introduction of oxygen-containing functional groups to the surface of PLLA. Raman spectroscopy indicated chain scission in the near surface region of PLLA (as predicted). However, e-beam effects on surface properties were not shown to be dependent on beam energy or dose. E-beam irradiation did not seem to affect the surface roughness of PLLA as a direct consequence of the treatment.

Exchange potential from the common energy denominator approximation for the Kohn-Sham Green's function: Application to (hyper)polarizabilities of molecular chains

An approximate Kohn-Sham (KS) exchange potential v(xsigma)(CEDA) is developed, based on the common energy denominator approximation (CEDA) for the static orbital Green's function, which preserves the essential structure of the density response function. v(xsigma)(CEDA) is an explicit functional of the occupied KS orbitals, which has the Slater v(Ssigma) and response v(respsigma)(CEDA) potentials as its components. The latter exhibits the characteristic step structure with "diagonal" contributions from the orbital densities \psi(isigma)\(2), as well as "off-diagonal" ones from the occupied-occupied orbital products psi(isigma)psi(j(not equal1)sigma). Comparison of the results of atomic and molecular ground-state CEDA calculations with those of the Krieger-Li-Iafrate (KLI), exact exchange (EXX), and Hartree-Fock (HF) methods show, that both KLI and CEDA potentials can be considered as very good analytical "closure approximations" to the exact KS exchange potential. The total CEDA and KLI energies nearly coincide with the EXX ones and the corresponding orbital energies epsilon(isigma) are rather close to each other for the light atoms and small molecules considered. The CEDA, KLI, EXX-epsilon(isigma) values provide the qualitatively correct order of ionizations and they give an estimate of VIPs comparable to that of the HF Koopmans' theorem. However, the additional off-diagonal orbital structure of v(xsigma)(CEDA) appears to be essential for the calculated response properties of molecular chains. KLI already considerably improves the calculated (hyper)polarizabilities of the prototype hydrogen chains H-n over local density approximation (LDA) and standard generalized gradient approximations (GGAs), while the CEDA results are definitely an improvement over the KLI ones. The reasons of this success are the specific orbital structures of the CEDA and KLI response potentials, which produce in an external field an ultranonlocal field-counteracting exchange potential. (C) 2002 American Institute of Physics.