A simulation study of the mixing, atomic flow and velocity profiles of crystalline materials during sliding

Self-contained Non-Equilibrium Molecular Dynamics (NEMD) simulations using Lennard-Jones potentials were performed to identify the origin and mechanisms of atomic scale interfacial behavior between sliding metals. The mixing sequence and velocity profiles were compared via MD simulations for three cases, viz.: sell-mated, similar and hard-softvcrystal pairs. The results showed shear instability, atomic scale mixing, and generation of eddies at the sliding interface. Vorticity at the interface suggests that atomic flow during sliding is similar to fluid flow under Kelvin-Helmholtz instability and this is supported by velocity profiles from the simulations. The initial step-function velocity profile spreads during sliding. However the velocity profile does not change much at later stages of the simulation and it eventually stops spreading. The steady state friction coefficient during simulation was monitored as a function of sliding velocity. Frictional behavior can be explained on the basis of plastic deformation and adiabatic effects. The mixing layer growth kinetics was also investigated.

Molecular dynamics simulations of sliding in an Fe-Cu tribopair system

Experimental evidence suggests that high strain rates, stresses, strains and temperatures are experienced near sliding interfaces. The associated microstructural changes are due to several dynamic an interacting phenomena. 3D non-equilibrium molecular dynamics (MD) simulations of sliding were conducted with the aim of understanding the dynamic processes taking place in crystalline tribopairs, with a focus on plastic deformation and microstructural evolution. Embedded atom potentials were employed for simulating sliding of an Fe-Cu tribopair. Sliding velocity, crystal orientation and presence of lattice defects were some of the variables in these simulations. Extensive plastic deformation involving dislocation and twin activity, dynamic recrystallization, amorphization and/or nanocrystallization, mechanical mixing and material transfer were observed. Mechanical mixing in the vicinity of the sliding interface was observed even in the Fe-Cu system, which would cluster under equilibrium conditions, hinting at the ballistic nature of the process. Flow localization was observed at high velocities implying the possible role of adiabatic heating. The presence of preexisting defects (such as dislocations and interfaces) played a pivotal role in determining friction and microstructural evolution. The study also shed light on the relationship between adhesion and plastic deformation, and friction. Comparisons with experiments suggest that such simulations can indeed provide valuable insights that are difficult to obtain from experiments.

Long waves in inviscid compressible atmosphere II

Solitary waves have been found in an adiabatic compressible atmosphere which, in ambient state, has winds and temperature gradient, generalizing our earlier results for the isothermal atmosphere. Explicit results are obtained for the special case of linear temperature and linear wind distributions in the undisturbed conditions. An important result of the study is that the number of possible critical speeds of the flow depends crucially on whether the maximum Richardson number (which is variable in the present example) is greater or less than 1/4.

Hydration numbers of some metal acetates, monochloroacetates and trichloroacetates in solution from ultrasonic velocity and conpressibility measurements

Ultrasonic velocities in aqueous solutions of some metal acetates, monochloroacelates and trichloroacetates, and the respective acids have been measured at 1 MHz frequency using the pulse technique. The ultrsonic velocity, adiabatic compressibility and apperent molal compressibility were measured as a function of concentration. The apparent molal compressibility values at infinite dilution were calculated and used to determine the hydration numbers.

Hot deformation behaviour of Mg–3Al alloy—A study using processing map

The hot deformation behaviour of Mg–3Al alloy has been studied using the processing-map technique. Compression tests were conducted in the temperature range 250–550 °C and strain rate range 3 × 10−4 to 102 s−1 and the flow stress data obtained from the tests were used to develop the processing map. The various domains in the map corresponding to different dissipative characteristics have been identified as follows: (i) grain boundary sliding (GBS) domain accommodated by slip controlled by grain boundary diffusion at slow strain-rates (<10−3 s−1) in the temperature range from 350 to 450 °C, (ii) two different dynamic recrystallization (DRX) domains with a peak efficiency of 42% at 550 °C/10−1 s−1 and 425 °C/102 s−1 governed by stress-assisted cross-slip and thermally activated climb as the respective rate controlling mechanisms and (iii) dynamic recovery (DRV) domain below 300 °C in the intermediate strain rate range from 3 × 10−2 to 3 × 10−1 s−1. The regimes of flow instability have also been delineated in the processing map using an instability criterion. Adiabatic shear banding at higher strain rates (>101 s−1) and solute drag by substitutional Al atoms at intermediate strain rates (3 × 10−2 to 3 × 10−1 s−1) in the temperature range (350–450 °C) are responsible for flow instability. The relevance of these mechanisms with reference to hot working practice of the material has been indicated. The processing maps of Mg–3Al alloy and as-cast Mg have been compared qualitatively to elucidate the effect of alloying with aluminum on the deformation behaviour of magnesium.

Primordial Isocurvature Perturbation in Light of CMB and LSS Data

The increased accuracy in the cosmological observations, especially in the measurements of the comic microwave background, allow us to study the primordial perturbations in grater detail. In this thesis, we allow the possibility for a correlated isocurvature perturbations alongside the usual adiabatic perturbations. Thus far the simplest six parameter \Lambda CDM model has been able to accommodate all the observational data rather well. However, we find that the 3-year WMAP data and the 2006 Boomerang data favour a nonzero nonadiabatic contribution to the CMB angular power sprctrum. This is primordial isocurvature perturbation that is positively correlated with the primordial curvature perturbation. Compared with the adiabatic \Lambda CMD model we have four additional parameters describing the increased complexity if the primordial perturbations. Our best-fit model has a 4% nonadiabatic contribution to the CMB temperature variance and the fit is improved by \Delta\chi^2 = 9.7. We can attribute this preference for isocurvature to a feature in the peak structure of the angular power spectrum, namely, the widths of the second and third acoustic peak. Along the way, we have improved our analysis methods by identifying some issues with the parametrisation of the primordial perturbation spectra and suggesting ways to handle these. Due to the improvements, the convergence of our Markov chains is improved. The change of parametrisation has an effect on the MCMC analysis because of the change in priors. We have checked our results against this and find only marginal differences between our parametrisation.

Non-Gaussian Cosmological Perturbations from Hybrid Inflation and Preheating

This thesis consists of four research papers and an introduction providing some background. The structure in the universe is generally considered to originate from quantum fluctuations in the very early universe. The standard lore of cosmology states that the primordial perturbations are almost scale-invariant, adiabatic, and Gaussian. A snapshot of the structure from the time when the universe became transparent can be seen in the cosmic microwave background (CMB). For a long time mainly the power spectrum of the CMB temperature fluctuations has been used to obtain observational constraints, especially on deviations from scale-invariance and pure adiabacity. Non-Gaussian perturbations provide a novel and very promising way to test theoretical predictions. They probe beyond the power spectrum, or two point correlator, since non-Gaussianity involves higher order statistics. The thesis concentrates on the non-Gaussian perturbations arising in several situations involving two scalar fields, namely, hybrid inflation and various forms of preheating. First we go through some basic concepts -- such as the cosmological inflation, reheating and preheating, and the role of scalar fields during inflation -- which are necessary for the understanding of the research papers. We also review the standard linear cosmological perturbation theory. The second order perturbation theory formalism for two scalar fields is developed. We explain what is meant by non-Gaussian perturbations, and discuss some difficulties in parametrisation and observation. In particular, we concentrate on the nonlinearity parameter. The prospects of observing non-Gaussianity are briefly discussed. We apply the formalism and calculate the evolution of the second order curvature perturbation during hybrid inflation. We estimate the amount of non-Gaussianity in the model and find that there is a possibility for an observational effect. The non-Gaussianity arising in preheating is also studied. We find that the level produced by the simplest model of instant preheating is insignificant, whereas standard preheating with parametric resonance as well as tachyonic preheating are prone to easily saturate and even exceed the observational limits. We also mention other approaches to the study of primordial non-Gaussianities, which differ from the perturbation theory method chosen in the thesis work.

Primordial Perturbations from a Self-interacting Curvaton

Inflation is a period of accelerated expansion in the very early universe, which has the appealing aspect that it can create primordial perturbations via quantum fluctuations. These primordial perturbations have been observed in the cosmic microwave background, and these perturbations also function as the seeds of all large-scale structure in the universe. Curvaton models are simple modifications of the standard inflationary paradigm, where inflation is driven by the energy density of the inflaton, but another field, the curvaton, is responsible for producing the primordial perturbations. The curvaton decays after inflation as ended, where the isocurvature perturbations of the curvaton are converted into adiabatic perturbations. Since the curvaton must decay, it must have some interactions. Additionally realistic curvaton models typically have some self-interactions. In this work we consider self-interacting curvaton models, where the self-interaction is a monomial in the potential, suppressed by the Planck scale, and thus the self-interaction is very weak. Nevertheless, since the self-interaction makes the equations of motion non-linear, it can modify the behaviour of the model very drastically. The most intriguing aspect of this behaviour is that the final properties of the perturbations become highly dependent on the initial values. Departures of Gaussian distribution are important observables of the primordial perturbations. Due to the non-linearity of the self-interacting curvaton model and its sensitivity to initial conditions, it can produce significant non-Gaussianity of the primordial perturbations. In this work we investigate the non-Gaussianity produced by the self-interacting curvaton, and demonstrate that the non-Gaussianity parameters do not obey the analytically derived approximate relations often cited in the literature. Furthermore we also consider a self-interacting curvaton with a mass in the TeV-scale. Motivated by realistic particle physics models such as the Minimally Supersymmetric Standard Model, we demonstrate that a curvaton model within the mass range can be responsible for the observed perturbations if it can decay late enough.

Propagation of a spherical shock in an inhomogeneous self-gravitating or nongravitating system

The propagation of a shock wave of finite strength due to an explosion into inhomogeneous nongravitating and self-gravitating systems has been considered, using similarity principles, supposing that the density varies as an inverse power of distance from the centre of explosion. A large number of systems, characterised by different density exponents and different adiabatic coefficients of the gas have been considered for different shock strengths. The numerical integration from the shock inward has been continued to the surface of singularity where density tends to infinity and which acts like a piston in the self-gravitating case and to the surface where the velocity gradient tends to infinity in the nongravitating case. The effect of variation of shock strength, density exponent and adiabatic coefficient on the location of these singularities and on the distribution of flow parameters behind the shock has been studied. The initial energy of the system and the manner of release of the explosion energy influence strongly the flow behind the shock. The results have been graphically depicted.

Strong shock with radiation near the surface of a star

The propagation of a shock wave, originating in a stellar interior, is considered when it approaches the surface of the star and assumes a self-similar character, "forgetting" its initial conditions. The flow behind the shock is assumed to be spatially isothermal rather than adiabatic to simulate the conditions of large radiative transfer near the stellar surface. The adiabatic and isothermal flows behind such a shock are compared. The exact shock-propagation laws, obtained by solving the equations in similarity variables, for different values of the parameter δ in the undisturbed density law, ρ0 ∝ xδ, and γ, the ratio of specific heats, are compared with the approximate values calculated by Whitham's characteristic rule and the two show a generally good agreement.

Ultrasonic propagation through aqueous solutions of acetates

Using the pulse method in the range of 2 to 26Mc's the ultrasonic absorption, velocity and the adiabatic compressibility have been studied in eleven aqueous acetate solutions up to a concentration of 1 mole/litre. The substances studied are the acetates of lithium, sodium, potassium, ammonium, magnesium, calcium, strontium, barium, zinc, cadmium and lead. Absorption in mercuric acetate has been studied only at 2 and 6 Mc/s. Two regions of relaxation are noticed, one below 10 Mc/s and the other between 10 and 26 Mc/s. The first relaxation is ascribed to the dissociation reaction of the salt and the second one to the monomerdimer reaction of the acetic acid formed by the hydrolysis of the salt in water.

A note on the Berry phase for systems having one degree of freedom

A one-dimensional arbitrary system with quantum Hamiltonian H(q, p) is shown to acquire the 'geometric' phase gamma (C)=(1/2) contour integral c(Podqo-qodpo) under adiabatic transport q to q+q+qo(t) and p to p+po(t) along a closed circuit C in the parameter space (qo(t), po(t)). The non-vanishing nature of this phase, despite only one degree of freedom (q), is due ultimately to the underlying non-Abelian Weyl group. A physical realisation in which this Berry phase results in a line spread is briefly discussed.

Annealing effect on transport properties of Nd0.67Sr0.33MnO3 thin films

Annealing dependence of the lattice parameter, resistivity, magnetoresistance and thermopower have been studied on Nd0.87Sr0.33MnO3 thin films deposited on LaAlO3 and alumina substrates by pulsed laser ablation. Upon annealing at 800 degrees C and 1000 degrees C the lattice constant of the LaAlO3 film tends toward that of the bulk target due to reduction in oxygen vacancies. This results in a metal-insulator transition at temperatures which increase with progressive annealing along with a decrease in the observed low temperature MR. Using a magnon scattering model we estimate the e(g) bandwidth of the film annealed at 1000 degrees C and show that the magnon contribution to the resistivity is suppressed in a highly oxygen deficient film and gains prominence only upon subsequent annealing. We also show that upon annealing, the polaron concentration and the spin cluster size increases in the paramagnetic phase, using an adiabatic polaron hopping model which takes into account an exchange dependent activation energy above the resistivity peak.

Time Dependent Operator-split and Unsplit Schemes for One Dimensional Premixed Flames

This paper is concerned with a study of an operator split scheme and unsplit scheme for the computation of adiabatic freely propagating one-dimensional premixed flames. The study uses unsteady method for both split and unsplit schemes employing implicit chemistry and explicit diffusion, a combination which is stable and convergent. Solution scheme is not sensitive to the initial starting estimate and provides steady state even with straight line profiles (far from steady state) in small number of time steps. Two systems H2-Air and H2-NO (involving complex nitrogen chemistry) are considered in presentinvestigation. Careful comparison shows that the operator split approach is slightly superior than the unsplit when chemistry becomes complex. Comparison of computational times with those of existing steady and unsteady methods seems to suggest that the method employing implicit-explicit algorithm is very efficient and robust.

Nonadiabatic charge pumping by oscillating potentials in one dimension: Results for infinite system and finite ring

We study charge pumping when a combination of static potentials and potentials oscillating with a time period T is applied in a one-dimensional system of noninteracting electrons. We consider both an infinite system using the Dirac equation in the continuum approximation and a periodic ring with a finite number of sites using the tight-binding model. The infinite system is taken to be coupled to reservoirs on the two sides which are at the same chemical potential and temperature. We consider a model in which oscillating potentials help the electrons to access a transmission resonance produced by the static potentials and show that nonadiabatic pumping violates the simple sin phi rule which is obeyed by adiabatic two-site pumping. For the ring, we do not introduce any reservoirs, and we present a method for calculating the current averaged over an infinite time using the time evolution operator U(T) assuming a purely Hamiltonian evolution. We analytically show that the averaged current is zero if the Hamiltonian is real and time-reversal invariant. Numerical studies indicate another interesting result, namely, that the integrated current is zero for any time dependence of the potential if it is applied to only one site. Finally we study the effects of pumping at two sites on a ring at resonant and nonresonant frequencies, and show that the pumped current has different dependences on the pumping amplitude in the two cases.