Modelling the nano-scale phenomena in condensed matter physics via computer-based numerical simulations

A review of the atomistic modelling of the behaviour of nano-scale structures and processes via molecular dynamics (MD) simulation method of a canonical ensemble is presented. Three areas of application in condensed matter physics are considered. We focus on the adhesive and indentation properties of the solid surfaces in nano-contacts, the nucleation and growth of nano-phase metallic and semi-conducting atomic and molecular films on supporting substrates, and the nano- and multi-scale crack propagation properties of metallic lattices. A set of simulations selected from these fields are discussed, together with a brief introduction to the methodology of the MD simulation. The pertinent inter-atomic potentials that model the energetics of the metallic and semi-conducting systems are also given.

One Subject, Two Lands: My Journey in Condensed Matter Physics

This is an account of a professional life in the field that was generally known as solid-state physics when I started working in it; India and the United States of America are the countries in which this life was largely played out. My attempts to understand various things in condensed matter physics, and efforts to put together people and activities in India in this field, are mainly the story.

Interacting fermionic atoms in optical lattices - A quantum simulator for condensed matter physics

This thesis reports on the creation and analysis of many-body states of interacting fermionic atoms in optical lattices. The realized system can be described by the Fermi-Hubbard hamiltonian, which is an important model for correlated electrons in modern condensed matter physics. In this way, ultra-cold atoms can be utilized as a quantum simulator to study solid state phenomena. The use of a Feshbach resonance in combination with a blue-detuned optical lattice and a red-detuned dipole trap enables an independent control over all relevant parameters in the many-body hamiltonian. By measuring the in-situ density distribution and doublon fraction it has been possible to identify both metallic and insulating phases in the repulsive Hubbard model, including the experimental observation of the fermionic Mott insulator. In the attractive case, the appearance of strong correlations has been detected via an anomalous expansion of the cloud that is caused by the formation of non-condensed pairs. By monitoring the in-situ density distribution of initially localized atoms during the free expansion in a homogeneous optical lattice, a strong influence of interactions on the out-of-equilibrium dynamics within the Hubbard model has been found. The reported experiments pave the way for future studies on magnetic order and fermionic superfluidity in a clean and well-controlled experimental system.

Relative energetics and structural properties of zirconia using a self-consistent tight-binding model

We describe an empirical, self-consistent, orthogonal tight-binding model for zirconia, which allows for the polarizability of the anions at dipole and quadrupole levels and for crystal field splitting of the cation d orbitals, This is achieved by mixing the orbitals of different symmetry on a site with coupling coefficients driven by the Coulomb potentials up to octapole level. The additional forces on atoms due to the self-consistency and polarizabilities are exactly obtained by straightforward electrostatics, by analogy with the Hellmann-Feynman theorem as applied in first-principles calculations. The model correctly orders the zero temperature energies of all zirconia polymorphs. The Zr-O matrix elements of the Hamiltonian, which measure covalency, make a greater contribution than the polarizability to the energy differences between phases. Results for elastic constants of the cubic and tetragonal phases and phonon frequencies of the cubic phase are also presented and compared with some experimental data and first-principles calculations. We suggest that the model will be useful for studying finite temperature effects by means of molecular dynamics.

Effect of quantization of vibrations on the structural properties of crystals

We study the structural effects produced by the quantization of vibrational degrees of freedom in periodic crystals at zero temperature. To this end we introduce a methodology based on mapping a suitable subspace of the vibrational manifold and solving the Schrödinger equation in it. A number of increasingly accurate approximations ranging from the quasiharmonic approximation (QHA) to the vibrational self-consistent field (VSCF) method and the exact solution are described. A thorough analysis of the approximations is presented for model monatomic and hydrogen-bonded chains, and results are presented for a linear H-F chain where the potential-energy surface is obtained via first-principles electronic structure calculations. We focus on quantum nuclear effects on the lattice constant and show that the VSCF is an excellent approximation, meaning that correlation between modes is not extremely important. The QHA is excellent for covalently bonded mildly anharmonic systems, but it fails for hydrogen-bonded ones. In the latter, the zero-point energy exhibits a nonanalytic behavior at the lattice constant where the H atoms center, which leads to a spurious secondary minimum in the quantum-corrected energy curve. An inexpensive anharmonic approximation of noninteracting modes appears to produce rather good results for hydrogen-bonded chains for small system sizes. However, it converges to the incorrect QHA results for increasing size. Isotope effects are studied for the first-principles H-F chain. We show how the lattice constant and the H-F distance increase with decreasing mass and how the QHA proves to be insufficient to reproduce this behavior.

The Effect of Antimony Substitution on the Magnetic and Structural Properties of Fe0.75–xSi0.25Sbx Alloys

The results of the investigation of the magnetic and structural properties of the alloy system Fe0.75–xSi0.25Sbx, where x = 0, 0.05, 0.1, 0.15, 0.2, and 0.25 synthesized by mechanical alloying followed by heat treatment are described. The x-ray diffraction reveals that all samples crystallize in the DO3-type cubic phase structure. Substituting Fe by Sb led to a de-crease in the lattice constant and the unit cell volume. The magnetic properties are investigated by vibrating sample magnetometer and show that all the samples are ferromagnetically ordered at room temperature. The Curie temperature is found to decrease linearly from (850 ± 5) K for the parent alloy to (620 ± 5) K for the alloyith x = 0.25. The satura-tion magnetizations at room temperature and at 100 K are found to decrease with increasing the antimony concentration. The above results indicate that Sb dissolves in the cubic structure of this alloy system.

Thermodynamic properties and atomistic structure of the dry amorphous silica surface from a reactive force field model

A force field model of the Keating type supplemented by rules to break, form, and interchange bonds is applied to investigate thermodynamic and structural properties of the amorphous SiO2 surface. A simulated quench from the liquid phase has been carried out for a silica sample made of 3888 silicon and 7776 oxygen atoms arranged on a slab similar to 40 angstrom thick, periodically repeated along two directions. The quench results into an amorphous sample, exposing two parallel square surfaces of similar to 42 nm(2) area each. Thermal averages computed during the quench allow us to determine the surface thermodynamic properties as a function of temperature. The surface tension turns out to be gamma=310 +/- 20 erg/cm(2) at room temperature and gamma=270 +/- 30 at T=2000 K, in fair agreement with available experimental estimates. The entropy contribution Ts-s to the surface tension is relatively low at all temperatures, representing at most similar to 20% of the surface energy. Almost without exceptions, Si atoms are fourfold coordinated and oxygen atoms are twofold coordinated. Twofold and threefold rings appear only at low concentration and are preferentially found in proximity of the surface. Above the glass temperature T-g=1660 +/- 50 K, the mobility of surface atoms is, as expected, slightly higher than that of bulk atoms. The computation of the height-height correlation function shows that the silica surface is rough in the equilibrium and undercooled liquid phase, becoming smooth below the glass temperature T-g.

VESUVIO : a novel instrument for performing spectroscopic studies in condensed matter with eV neutrons at the ISIS facility

The VESUVIO project aims to provide unique prototype instrumentation at the ISIS-pulsed neutron source and to establish a routine experimental and theoretical program in neutron scattering spectroscopy at eV energies. This instrumentation will be specifically designed for high momentum, , and energy transfer inelastic neutron scattering studies of microscopic dynamical processes in materials and will represent a unique facility for EU researchers. It will allow to derive single-particle kinetic energies and single-particle momentum distributions, n(p), providing additional and/or complementary information to other neutron inelastic spectroscopic techniques.

Structural and chemical embrittlement of grain boundaries by impurities: A general theory and first-principles calculations for copper

First-principles calculations of the Sigma 5(310)[001] symmetric tilt grain boundary in Cu with Bi, Na, and Ag substitutional impurities provide evidence that in the phenomenon of Bi embrittlement of Cu grain boundaries electronic effects do not play a major role; on the contrary, the embrittlement is mostly a structural or "size" effect. Na is predicted to be nearly as good an embrittler as Bi, whereas Ag does not embrittle the boundary in agreement with experiment. While we reject the prevailing view that "electronic" effects (i.e., charge transfer) are responsible for embrittlement, we do not exclude the role of chemistry. However, numerical results show a striking equivalence between the alkali metal Na and the semimetal Bi, small differences being accounted for by their contrasting "size" and "softness" (defined here). In order to separate structural and chemical effects unambiguously if not uniquely, we model the embrittlement process by taking the system of grain boundary and free surfaces through a sequence of precisely defined gedanken processes; each of these representing a putative mechanism. We thereby identify three mechanisms of embrittlement by substitutional impurities, two of which survive in the case of embrittlement or cohesion enhancement by interstitials. Two of the three are purely structural and the third contains both structural and chemical elements that by their very nature cannot be further unraveled. We are able to take the systems we study through each of these stages by explicit computer simulations and assess the contribution of each to the net reduction in intergranular cohesion. The conclusion we reach is that embrittlement by both Bi and Na is almost exclusively structural in origin; that is, the embrittlement is a size effect.

Structural phase transitions in low-dimensional ion crystals

A chain of singly charged particles, confined by a harmonic potential, exhibits a sudden transition to a zigzag configuration when the radial potential reaches a critical value, depending on the particle number. This structural change is a phase transition of second order, whose order parameter is the crystal displacement from the chain axis. We study analytically the transition using Landau theory and find full agreement with numerical predictions by Schiffer [Phys. Rev. Lett. 70, 818 (1993)] and Piacente [Phys. Rev. B 69, 045324 (2004)]. Our theory allows us to determine analytically the system's behavior at the transition point.

Electrical and structural properties of zirconia thin films prepared by reactive magnetron sputtering

Thin films of ZrO2 were prepared by reactive magnetron sputtering. Annealing of the films exhibited a drastic change in the properties due to improved crystallinity and packing density. The root mean square roughness of the sample observed from atomic force microscope is about 5.75 nm which is comparable to the average grain size of the thin film which is about 6 nm obtained from X-ray diffraction. The film annealed at 873 K exhibits an optical band gap of around 4.83 eV and shows +4 oxidation state of zirconium indicating fully oxidized zirconium, whereas higher annealing temperatures lead to oxygen deficiency in the films and this is reflected in their properties. A discontinuity in the imaginary part of the AC conductivity was observed in the frequency range of tens of thousands of Hz, where as, the real part does not show such behavior.