Early Detection of Fretting Corrosion in Hip Replacement by Acoustic Emission Non-Invasive Technique.

The present study aimed to investigate the feasibility of using acoustic emission (AE) as a detection method for identifying failure mechanisms at the modular junction interface in total hip replacements (THRs) subjected to fretting corrosion. The experimental setup involved simulating fretting corrosion using a Ti6Al4V disc representing the femoral neck and a ZrO2 pin representing the femoral head. Mechanical testing provided insights into the wear and frictional behavior occurring at the modular junction interface. The results revealed that for all three potential conditions, a fretting condition of partial slip was observed. These findings highlight the importance of understanding the mechanical interactions and their influence on the overall performance and longevity of THRs. Electrochemical analysis shed light on the corrosion behavior under different potentiostatic conditions. High potentials in the anodic condition led to increased corrosion and ion transfer due to the breakdown of the passive oxide layer. Conversely, the cathodic potential condition exhibited a regrowth of the passive oxide layer, protecting the Ti6Al4V surface from further corrosion. The mid-range corrosion potential condition showed a dynamic equilibrium between corrosion and passivation processes. These electrochemical insights enhance our understanding of the mechanisms involved in fretting corrosion. The AE data proved to be promising in detecting and monitoring the onset and progression of failure mechanisms. The AE signals exhibited distinctive patterns that correlated with the severity of fretting corrosion. Notably, the hit driven data results, derived from AE signals, demonstrated the ability to differentiate between different levels of fretting conditions. This suggests that AE can serve as a valuable diagnostic tool for early detection and continuous monitoring of implant failure in THRs.

Decoding ionic conductivity and reordering in cation-disordered pyrochlores.

The ordered structure A2B2O6O' in pyrochlores engenders twin rows of inequivalent anion sublattices each centred on alternating cations. While it is known that cation antisite disorder augments the ionic conductivity by several orders of magnitude, the local cation environment around the anions and the dynamic anion reordering during the cation disordering are not well-elucidated. Using atomistic simulations on Gd2Zr2O7, we first show that the anions engage in concerted hops to the neighbouring tetrahedral sites mostly along with the 〈1 0 0〉 direction while completely avoiding the octahedral sites. While the initially vacant 8a sites start accommodating oxygen ions with increasing cation disorder, they show noticeable reluctance even at significant levels of disorder. We have also tracked both the distribution of available oxygen sites following random cation disorder, which is dependent only on cation disordering, and the probability of occupation of these sites. Interestingly, the available oxygen sites show a non-monotonic dependence on the number of B ions in the nearest neighbouring shell while the occupation probability of all the available oxygen sites increases monotonically. A tetrahedral oxygen site thus has a better probability of being occupied when it has a greater number of second neighbour B ions. This article is part of the Theo Murphy meeting issue 'Understanding fast-ion conduction in solid electrolytes'.

Exact diagonal representation of normal mode energy, occupation number, and heat current for phonon-dominated thermal transport.

Collective excitations of crystal vibrations or normal modes are customarily described using complex normal mode coordinates. While appropriate for calculating phonon dispersion, the mixed representation involving the complex conjugates does not allow the construction of equivalent phonon occupation number or modal dynamical quantities such as the energy or heat current specific to a wave-vector direction (q). Starting from a canonical solution that includes waves going to the left and right directions, we cast the Hamiltonian, normal mode population, and heat current in an exactly diagonalizable representation using real normal mode amplitudes. We show that the use of real amplitudes obviates the need for a complex modal heat current while making the passage to second quantization more apparent. Using nonequilibrium molecular dynamics simulations, we then compute the net modal energy, heat current, and equivalent phonon population in a linear lattice subjected to a thermal gradient. Our analysis paves a tractable path for probing and computing the direction-dependent thermal-phononic modal properties of dielectric lattices using atomistic simulations.

Deducing Phonon Scattering from Normal Mode Excitations.

While the quantum scattering theory has provided the theoretical underpinning for phonon interactions, the correspondence between the phonon modes and normal modes of vibrations has never been fully established; for example, the nature of energy exchange during elementary normal mode interactions remains largely unknown. In this work, by adopting a set of real asymmetric normal mode amplitudes, we first discriminate the normal and Umklapp processes directly from atomistic dynamics. We then demonstrate that the undulating harmonic and anharmonic potentials, which allow a number of interaction pathways, generate several total-energy-conserving forward and backward scattering events including those which are traditionally considered as quantum-forbidden. Although the normal mode energy is proportional to the square of the eigen-frequency, we deduce that the energy exchanged from one mode to another in each elementary interaction is proportional to the frequency - a quantum-like restriction. We anticipate that the current approach can be utilized profitably to discover unbiased scattering channels, many traditionally quantum forbidden, with complex anharmonicities. Our discovery will aid in the development of next-generation Peierls-Boltzmann transport simulations that access normal mode scattering pathways from finite temperature ab initio simulations.

Low Dimensional String-like Relaxation Underpins Superionic Conduction in Fluorites and Related Structures.

Among the superionic conductors that show a Faraday transition - the continuous increase in the ionic conductivity over a range of temperatures - the fluorite structures have enjoyed incisive examinations over the past four decades; yet the fundamental nature of superionicity has remained largely inconclusive. Departing from the traditional quasi-static defect framework, we provide weighty evidence for string-like dynamical structures that govern the fast ion conduction process in fluorites. We show that lower temperatures encourage the growth of longer but slowly relaxing strings and vice-versa - a direct manifestation of heterogeneous dynamics. Remarkably, the ionic conductivity is inversely correlated to the lifetime of the ions that participate in the strings and not explicitly to the ion population. Our analysis methodology, which resolves a long-standing disagreement on defect structures and the mechanism of ionic transport in fcc fluorite structures, is well-positioned to describe the dynamics of low dimensional conduction in a larger class of superionic conductors.

Mobility propagation and dynamic facilitation in superionic conductors.

In an earlier work [V. A. Annamareddy et al., Phys. Rev. E 89, 010301(R) (2014)], we showed the manifestation of dynamical heterogeneity (DH)-the presence of clustered mobile and immobile regions-in UO2, a model type II superionic conductor. In the current work, we demonstrate the mechanism of dynamic facilitation (DF) in two superionic conductors (CaF2 and UO2) using atomistic simulations. Using the mobility transfer function, DF is shown to vary non-monotonically with temperature with the intensity of DF peaking at temperatures close to the superionic transition temperature (Tλ). Both the metrics quantifying DH and DF show remarkable correspondence implying that DF, in the framework of kinematically constrained models, underpins the heterogeneous dynamics in type II superionic conductors.

Waxing and waning of dynamical heterogeneity in the superionic state.

Using molecular dynamics simulations of UO2-a type II superionic conductor-we identify a well-defined onset of dynamic disorder (Tα), which is remarkably correlated to a nontrivial advance of dynamical heterogeneity (DH). Quantified by the correlations in the dynamic propensity and van Hove self-correlation function, the DH is shown to grow with increasing temperature from Tα, peak at an intermediate temperature between Tα and Tλ-the superionic transition temperature-and then recede. Surprisingly, the DH attributes are not uniform across the temperatures-our investigation shows a low temperature (αT) stage DH, which is characterized by weak correlations and a plateaulike period in the correlations of the propensity, and a high temperature (λT) stage DH with strong correlations that are analogous to those in typical supercooled liquids. Our work, which has rigorously identified the onset of superionicity, gives a different direction for interpreting scattering experiments on the basis of statistical, correlated dynamics.

Absolute thermodynamic properties of molten salts using the two-phase thermodynamic (2PT) superpositioning method.

We show that the absolute thermodynamic properties of molten salts (mixtures of KCl and LiCl) can be accurately determined from the two-phase thermodynamic (2PT) method that is based on superpositioning of solid-like and gas-like (hard-sphere) vibrational density of states (DoS). The 2PT predictions are in excellent accordance with those from the thermodynamic integration method; the melting point of KCl evaluated from the free energy and the absolute entropy shows close conformity with the experimental/NIST data. The DoS partitioning shows that the Li(+) ions in the eutectic LiCl-KCl molten mixture are largely solid-like, unlike the K(+) and Cl(-) ions, which have a significant gas-like contribution, for temperatures ranging from 773 K to 1300 K. The solid-like states of the Li(+) ions may have practical implications when employed for chemical and nuclear reprocessing applications.

Multicomponent diffusion in molten LiCl-KCl: dynamical correlations and divergent Maxwell-Stefan diffusivities.

Multicomponent diffusional mechanisms in the ternary LiCl-KCl system are elucidated using the Green-Kubo formalism and equilibrium molecular dynamics simulations. The Maxwell-Stefan (MS) diffusion matrix is evaluated from the Onsager dynamical matrix that contains the diffusion flux correlation functions. From the temporal behavior of the correlation functions, we observe that the Li-Li and Li-Cl ion pairs have a pronounced cage dynamics that remains noticeably strong even at high temperatures. Even though the Onsager coefficients, which are the time integrals of the diffusion flux correlation functions, portray a relatively smooth variation across various compositions and temperatures, we observe a sign change and a divergent-like behavior for the MS diffusivity of the K-Li ion pair at a temperature of ~1100 K for the eutectic composition, and at a KCl mole fraction of ~0.49 at 1043 K. Negative MS diffusivities, while unusual, are however shown to satisfy the nonnegative entropic constraints.

Computing the viscosity of supercooled liquids: Markov Network model.

The microscopic origin of glass transition, when liquid viscosity changes continuously by more than ten orders of magnitude, is challenging to explain from first principles. Here we describe the detailed derivation and implementation of a Markovian Network model to calculate the shear viscosity of deeply supercooled liquids based on numerical sampling of an atomistic energy landscape, which sheds some light on this transition. Shear stress relaxation is calculated from a master-equation description in which the system follows a transition-state pathway trajectory of hopping among local energy minima separated by activation barriers, which is in turn sampled by a metadynamics-based algorithm. Quantitative connection is established between the temperature variation of the calculated viscosity and the underlying potential energy and inherent stress landscape, showing a different landscape topography or "terrain" is needed for low-temperature viscosity (of order 10(7) Pa·s) from that associated with high-temperature viscosity (10(-5) Pa·s). Within this range our results clearly indicate the crossover from an essentially Arrhenius scaling behavior at high temperatures to a low-temperature behavior that is clearly super-Arrhenius (fragile) for a Kob-Andersen model of binary liquid. Experimentally the manifestation of this crossover in atomic dynamics continues to raise questions concerning its fundamental origin. In this context this work explicitly demonstrates that a temperature-dependent "terrain" characterizing different parts of the same potential energy surface is sufficient to explain the signature behavior of vitrification, at the same time the notion of a temperature-dependent effective activation barrier is quantified.

Computing the viscosity of supercooled liquids. II. Silica and strong-fragile crossover behavior.

A recently developed atomistic method capable of calculating the fragile (non-Arrhenius) temperature behavior of highly viscous liquids is further tested by studying a model of SiO(2), a glass former well known for its Arrhenius temperature behavior (strong). The method predicts an Arrhenius temperature variation, in agreement with experiments, the origin of which is revealed by both quantitative and qualitative results on transition state pathways, activation barrier analysis, energy landscape connectivity, and atomistic activation mechanisms. Also predicted is a transition from fragile to strong behavior at a lower viscosity, below the range of measurements, which had been previously suggested on the basis of molecular dynamics simulations. By systematically comparing our findings with corresponding results on the binary Lennard-Jones system (fragile) we gain new insights into the topographical features of the potential energy landscape, characteristics that distinguish strong from fragile glassy systems. We interpret fragility as a universal manifestation of slowing of dynamics when the system becomes trapped in deep energy basins. As a consequence, all glass-forming systems, when cooled from their normal liquid state, should exhibit two transitions in temperature scaling of the viscosity, a strong-to-fragile crossover followed by a second transition reverting back to strong behavior.

Computing the viscosity of supercooled liquids.

We describe an atomistic method for computing the viscosity of highly viscous liquids based on activated state kinetics. A basin-filling algorithm allowing the system to climb out of deep energy minima through a series of activation and relaxation is proposed and first benchmarked on the problem of adatom diffusion on a metal surface. It is then used to generate transition state pathway trajectories in the potential energy landscape of a binary Lennard-Jones system. Analysis of a sampled trajectory shows the system moves from one deep minimum to another by a process that involves high activation energy and the crossing of many local minima and saddle points. To use the trajectory data to compute the viscosity we derive a Markov Network model within the Green-Kubo formalism and show that it is capable of producing the temperature dependence in the low-viscosity regime described by molecular dynamics simulation, and in the high-viscosity regime (10(2)-10(12) Pa s) shown by experiments on fragile glass-forming liquids. We also derive a mean-field-like description involving a coarse-grained temperature-dependent activation barrier, and show it can account qualitatively for the fragile behavior. From the standpoint of molecular studies of transport phenomena this work provides access to long relaxation time processes beyond the reach of current molecular dynamics capabilities. In a companion paper we report a similar study of silica, a representative strong liquid. A comparison of the two systems gives insight into the fundamental difference between strong and fragile temperature variations.

Mean-field versus microconvection effects in nanofluid thermal conduction.

Transient hot-wire data on thermal conductivity of suspensions of silica and perfluorinated particles show agreement with the mean-field theory of Maxwell but not with the recently postulated microconvection mechanism. The influence of interfacial thermal resistance, convective effects at microscales, and the possibility of thermal conductivity enhancements beyond the Maxwell limit are discussed.

Mechanism of thermal transport in dilute nanocolloids.

Thermal conduction modes in a nanocolloid (nanofluid) are quantitatively assessed by combining linear response theory with molecular dynamics simulations. The microscopic heat flux is decomposed into three additive fluctuation modes, namely, kinetic, potential, and collision. For low volume fractions (<1%) of nanosized platinum clusters which interact strongly with xenon host liquid, a significant thermal conductivity enhancement results from the self correlation in the potential flux. Our findings reveal a molecular-level mechanism for enhanced thermal conductivity in nanocolloids with short-ranged attraction and offer predictions that can be experimentally tested.

Beyond the Maxwell limit: thermal conduction in nanofluids with percolating fluid structures.

In a well-dispersed nanofluid with strong cluster-fluid attraction, thermal conduction paths can arise through percolating amorphouslike interfacial structures. This results in a thermal conductivity enhancement beyond the Maxwell limit of 3phi, with phi being the nanoparticle volume fraction. Our findings from nonequilibrium molecular dynamics simulations, which are amenable to experimental verification, can provide a theoretical basis for the development of future nanofluids.

Statistical field estimators for multiscale simulations.

We present a systematic approach for generating smooth and accurate fields from particle simulation data using the notions of statistical inference. As an extension to a parametric representation based on the maximum likelihood technique previously developed for velocity and temperature fields, a nonparametric estimator based on the principle of maximum entropy is proposed for particle density and stress fields. Both estimators are applied to represent molecular dynamics data on shear-driven flow in an enclosure which exhibits a high degree of nonlinear characteristics. We show that the present density estimator is a significant improvement over ad hoc bin averaging and is also free of systematic boundary artifacts that appear in the method of smoothing kernel estimates. Similarly, the velocity fields generated by the maximum likelihood estimator do not show any edge effects that can be erroneously interpreted as slip at the wall. For low Reynolds numbers, the velocity fields and streamlines generated by the present estimator are benchmarked against Newtonian continuum calculations. For shear velocities that are a significant fraction of the thermal speed, we observe a form of shear localization that is induced by the confining boundary.