Structural phase transition study of FePt alloys using ab initio calculation

The FePt alloy undergoes the cubic to tetragonal lattice transformation in the ferromagnetic state. We calculated the electronic structure for both cubic and tetragonal structures using the FPLAPW method with APW + lo. Comparing the density of states of the cubic and tetragonal structures, it is expected that the lattice transformation is caused by the band Jahn-Teller effect. (C) 2009 Elsevier B.V. All rights reserved.

Ultrasonic material characterization using large-aperture PVDF receivers

This work describes the use of a large-aperture PVDF receiver in the measurement of liquid density and composite material elastic constants. The density measurement of several liquids is obtained with accuracy of 0.2% using a conventional NDE emitter transducer and a 70-mm-diameter, 52-mu m P(VDF-TrFE) membrane with gold electrodes. The determination of the elastic constants is based on the phase velocity measurement. Diffraction can lead to errors around 1% in velocity measurement when using alternatively the conventional pair of ultrasonic transducers (1-MHz frequency and 19-mm-diameter) operating in through-transmission mode, separated by a distance of 100 mm. This effect is negligible when using a pair of 10-MHz, 19-mm-diameter transducers. Nevertheless, the dispersion at 10 MHz can result in errors of about 0.5%, when measuring the velocity in composite materials. The use of an 80-mm diameter, 52-mu m-thick PVDF membrane receiver practically eliminates the diffraction effects in phase velocity measurement. The elastic constants of a carbon fiber reinforced polymer were determined and compared with the values obtained by a tensile test. (C) 2009 Elsevier B. V. All rights reserved.

Propriedades eletrônicas, estruturais e constantes elásticas do ZnO

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Ultrasonic material characterization using diffraction-free PVDF receivers

This work describes the use of a large aperture PVDF receiver in the measurement of density of liquids and elastic constants of composite materials. The density measurement of several liquids is obtained with the accuracy of less than 0.2% using a conventional NDT emitter transducer and a 70-mm diameter, 52-μm P(VDF-TrFE) membrane with gold electrodes. The determination of the elastic constants of composite materials is based in the measurement of phase velocity. It is shown that the diffraction can lead to errors around 1% in the velocity measurement when using a pair of ultrasonic transducers (1MHz and 19mm diameter) operating in transmission-reception mode separated by a distance of 100 mm. This effect is negligible when using a pair of 10-MHz transducers. On the other hand, the dispersion at 10 MHz can result in errors of about 0.5%, measuring the velocity in composite materials. The use of an 80-mm diameter, 52-μm thick PVDF membrane receiver allows measuring the phase velocity without the diffraction effects.

Hydrostatic and [001] uniaxial pressure on anatase TiO2 by periodic B3LYP-D* calculations

The effect of high hydrostatic and [001] uniaxial pressures on TiO 2 anatase was studied under the framework of periodic calculations with the inclusion of DFT-D2 dispersion potential adjusted for this system (B3LYP-D*). The role of dispersion in distorted unit cells was evaluated in terms of lattice parameters, elastic constants, equation of state, vibrational properties, and electronic properties (band structure and density of states). A more reliable description at high pressures was achieved because the B3LYP-D* presented an improvement in all properties for undistorted bulk over conventional B3LYP and B3LYP-D. From density of states analysis, we observed that the contribution of crystalline orbitals to the edge of valence and conduction bands changed within applied pressure. The studied distortions can give some insight into behavior of electronic and structural properties due to local stress in anatase bulk from doping, defects, and physical tensions in nanometric forms. © 2013 American Chemical Society.

Mechanical properties and fracture dynamics of silicene membranes

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Modelagem e Simulação Computacional das propriedades estruturais, eletrônicas e constantes elásticas do ZnO

The electronic, structural properties and elastic constants of the wurtzite phase of zinc oxide, ZnO, was investigated using computer simulation at Density Functional Theory level, with B3LYP hybrid functional and Hartree-Fock methodology. The electronic properties as well the band energy was investigated through the analysis of the band structures and density of states (DOS), and the mechanical properties was studied through the calculus of the elastic constants C11, C33, C44, C12 e C13. The results are in good agreement with experimental data found in the literature and in accordance with results obtained by another theoretical methodology

Inspeção de materiais compósitos utilizando ondas acústicas guiadas

Pós-graduação em Engenharia Elétrica - FEIS

Development of empirical potentials for amorphous silica

Development of empirical potentials for amorphous silica Amorphous silica (SiO2) is of great importance in geoscience and mineralogy as well as a raw material in glass industry. Its structure is characterized as a disordered continuous network of SiO4 tetrahedra. Many efforts have been undertaken to understand the microscopic properties of silica by classical molecular dynamics (MD) simulations. In this method the interatomic interactions are modeled by an effective potential that does not take explicitely into account the electronic degrees of freedom. In this work, we propose a new methodology to parameterize such a potential for silica using ab initio simulations, namely Car-Parrinello (CP) method [Phys. Rev. Lett. 55, 2471 (1985)]. The new potential proposed is compared to the BKS potential [Phys. Rev. Lett. 64, 1955 (1990)] that is considered as the benchmark potential for silica. First, CP simulations have been performed on a liquid silica sample at 3600 K. The structural features so obtained have been compared to the ones predicted by the classical BKS potential. Regarding the bond lengths the BKS tends to underestimate the Si-O bond whereas the Si-Si bond is overestimated. The inter-tetrahedral angular distribution functions are also not well described by the BKS potential. The corresponding mean value of theSiOSi angle is found to be ≃ 147◦, while the CP yields to aSiOSi angle centered around 135◦. Our aim is to fit a classical Born-Mayer/Coulomb pair potential using ab initio calculations. To this end, we use the force-matching method proposed by Ercolessi and Adams [Europhys. Lett. 26, 583 (1994)]. The CP configurations and their corresponding interatomic forces have been considered for a least square fitting procedure. The classical MD simulations with the resulting potential have lead to a structure that is very different from the CP one. Therefore, a different fitting criterion based on the CP partial pair correlation functions was applied. Using this approach the resulting potential shows a better agreement with the CP data than the BKS ones: pair correlation functions, angular distribution functions, structure factors, density of states and pressure/density were improved. At low temperature, the diffusion coefficients appear to be three times higher than those predicted by the BKS model, however showing a similar temperature dependence. Calculations have also been carried out on crystalline samples in order to check the transferability of the potential. The equilibrium geometry as well as the elastic constants of α-quartz at 0 K are well described by our new potential although the crystalline phases have not been considered for the parameterization. We have developed a new potential for silica which represents an improvement over the pair potentials class proposed so far. Furthermore, the fitting methodology that has been developed in this work can be applied to other network forming systems such as germania as well as mixtures of SiO2 with other oxides (e.g. Al2O3, K2O, Na2O).

Foliation boudinage

In this thesis foliation boudinage and related structures have been studied based on field observations and numerical modeling. Foliation boudinage occurs in foliated rocks independent of lithology contrast. The developing structures are called ‘Foliation boudinage structures (FBSs)’ and show evidence for both ductile and brittle deformation. They are recognized in rocks by perturbations in monotonous foliation adjacent to a central discontinuity, mostly filled with vein material. Foliation boudinage structures have been studied in the Çine Massif in SW-Turkey and the Furka Pass-Urseren Zone in central Switzerland. Four common types have been distinguished in the field, named after vein geometries in their boudin necks in sections normal to the boudin axis: lozenge-, crescent-, X- and double crescent- type FBSs. Lozengetype FBSs are symmetric and characterized by lozenge-shaped veins in their boudin neck with two cusps facing opposite sides. A symmetrical pair of flanking folds occurs on the two sides of the vein. Crescent-type FBSs are asymmetric with a single smoothly curved vein in the boudin neck, with vein contacts facing to one side. X- and double crescent- type FBSs are asymmetric. The geometry of the neck veins resembles that of cuspate-lobate structures. The geometry of flanking structures is related to the shape of the veins. The veins are mostly filled with massive quartz in large single crystals, commonly associated with tourmaline, feldspar and biotite and in some cases with chlorite. The dominance of large facetted single quartz crystals and spherulitic chlorite in the veins suggest that the minerals grew into open fluidfilled space. FLAC experiments show that fracture propagation during ductile deformation strongly influences the geometry of developing veins. The cusps of the veins are better developed in the case of propagating fractures. The shape of the boudin neck veins in foliation boudinage depends on the initial orientation and shape of the fracture, the propagation behaviour of the fracture, the geometry of bulk flow, and the stage at which mineral filling takes place. A two dimensional discrete element model was used to study the progressive development of foliation boudinage structures and the behavior of visco-elastic material deformed under pure shear conditions. Discrete elements are defined by particles that are connected by visco-elastic springs. Springs can break. A number of simulations was Abstract vii performed to investigate the effect of material properties (Young’s modulus, viscosity and breaking strength) and anisotropy on the developing structures. The models show the development of boudinage in single layers, multilayers and in anisotropic materials with random mica distribution. During progressive deformation different types of fractures develop from mode I, mode II to the combination of both. Voids develop along extension fractures, at intersections of conjugate shear fractures and in small pull-apart structures along shear fractures. These patterns look similar to the natural examples. Fractures are more localized in the models where the elastic constants are low and the competence contrast is high between the layers. They propagate through layers where the constants are high and the competence contrast is relatively low. Flow localize around these fractures and voids. The patterns similar to symmetric boudinage structures and extensional neck veins (e.g. lozenge type) more commonly develop in the models with lower elastic constants and anisotropy. The patterns similar to asymmetric foliation boudinage structures (e.g. X-type) develop associated with shear fractures in the models where elastic constants and anisotropy of the materials are relatively high. In these models boudin neck veins form commonly at pull-aparts along the shear fractures and at the intersection of fractures.

Hypersonic phonon propagation in mesoscopic systems by Brillouin spectroscopy

Phononic crystals, capable to block or direct the propagation of elastic/acoustic waves, have attracted increasing interdisciplinary interest across condensed matter physics and materials science. As of today, no generalized full description of elastic wave propagation in phononic structures is available, mainly due to the large number of variables determining the band diagram. Therefore, this thesis aims for a deeper understanding of the fundamental concepts governing wave propagation in mesoscopic structures by investigation of appropriate model systems. The phononic dispersion relation at hypersonic frequencies is directly investigated by the non-destructive technique of high-resolution spontaneous Brillouin light scattering (BLS) combined with computational methods. Due to the vector nature of the elastic wave propagation, we first studied the hypersonic band structure of hybrid superlattices. These 1D phononic crystals composed of alternating layers of hard and soft materials feature large Bragg gaps. BLS spectra are sensitive probes of the moduli, photo-elastic constants and structural parameters of the constituent components. Engineering of the band structure can be realized by introduction of defects. Here, cavity layers are employed to launch additional modes that modify the dispersion of the undisturbed superlattice, with extraordinary implications to the band gap region. Density of states calculations in conjunction with the associated deformation allow for unambiguous identication of surface and cavity modes, as well as their interaction with adjacent defects. Next, the role of local resonances in phononic systems is explored in 3D structures based on colloidal particles. In turbid media BLS records the particle vibration spectrum comprising resonant modes due to the spatial confinement of elastic energy. Here, the frequency and lineshapes of the particle eigenmodes are discussed as function of increased interaction and departure from spherical symmetry. The latter is realized by uniaxial stretching of polystyrene spheres, that can be aligned in an alternating electric field. The resulting spheroidal crystals clearly exhibit anisotropic phononic properties. Establishing reliable predictions of acoustic wave propagation, necessary to advance, e.g., optomechanics and phononic devices is the ultimate aim of this thesis.

Spin-strain coupling in a 3-D transition metal oxides

ab-initio Hartree Fock (HF), density functional theory (DFT) and hybrid potentials were employed to compute the optimized lattice parameters and elastic properties of perovskite 3-d transition metal oxides. The optimized lattice parameters and elastic properties are interdependent in these materials. An interaction is observed between the electronic charge, spin and lattice degrees of freedom in 3-d transition metal oxides. The coupling between the electronic charge, spin and lattice structures originates due to localization of d-atomic orbitals. The coupling between the electronic charge, spin and crystalline lattice also contributes in the ferroelectric and ferromagnetic properties in perovskites. The cubic and tetragonal crystalline structures of perovskite transition metal oxides of ABO3 are studied. The electronic structure and the physics of 3-d perovskite materials is complex and less well considered. Moreover, the novelty of the electronic structure and properties of these perovskites transition metal oxides exceeds the challenge offered by their complex crystalline structures. To achieve the objective of understanding the structure and property relationship of these materials the first-principle computational method is employed. CRYSTAL09 code is employed for computing crystalline structure, elastic, ferromagnetic and other electronic properties. Second-order elastic constants (SOEC) and bulk moduli (B) are computed in an automated process by employing ELASTCON (elastic constants) and EOS (equation of state) programs in CRYSTAL09 code. ELASTCON, EOS and other computational algorithms are utilized to determine the elastic properties of tetragonal BaTiO3, rutile TiO2, cubic and tetragonal BaFeO3 and the ferromagentic properties of 3-d transition metal oxides. Multiple methods are employed to crosscheck the consistency of our computational results. Computational results have motivated us to explore the ferromagnetic properties of 3-d transition metal oxides. Billyscript and CRYSTAL09 code are employed to compute the optimized geometry of the cubic and tetragonal crystalline structure of transition metal oxides of Sc to Cu. Cubic crystalline structure is initially chosen to determine the effect of lattice strains on ferromagnetism due to the spin angular momentum of an electron. The 3-d transition metals and their oxides are challenging as the basis functions and potentials are not fully developed to address the complex physics of the transition metals. Moreover, perovskite crystalline structures are extremely challenging with respect to the quality of computations as the latter requires the well established methods. Ferroelectric and ferromagnetic properties of bulk, surfaces and interfaces are explored by employing CRYSTAL09 code. In our computations done on cubic TMOs of Sc-Fe it is observed that there is a coupling between the crystalline structure and FM/AFM spin polarization. Strained crystalline structures of 3-d transition metal oxides are subjected to changes in the electromagnetic and electronic properties. The electronic structure and properties of bulk, composites, surfaces of 3-d transition metal oxides are computed successfully.

Comparison of Mixed and Kinematic Uniform Boundary Conditions in Homogenized Elasticity of Femoral Trabecular Bone Using Microfinite Element Analyses

Mechanical properties of human trabecular bone play an important role in age-related bone fragility and implant stability. Micro-finite element (microFE) analysis allows computing the apparent elastic properties of trabecular bone biopsies, but the results depend on the type of applied boundary conditions (BCs). In this study, 167 femoral trabecular cubic biopsies with a side length of 5.3 mm were analyzed using microFE analysis to compare their stiffness systematically with kinematic uniform boundary conditions (KUBCs) and periodicity-compatible mixed uniform boundary conditions (PMUBCs). The obtained elastic constants were then used in the volume fraction and fabric-based orthotropic Zysset-Curnier model to identify their respective model parameters. As expected, PMUBCs lead to more compliant apparent elastic properties than KUBCs, especially in shear. The differences in stiffness decreased with bone volume fraction and mean intercept length. Unlike KUBCs, PMUBCs were sensitive to heterogeneity of the biopsies. The Zysset-Curnier model predicted apparent elastic constants successfully in both cases with adjusted coefficients of determination of 0.986 for KUBCs and 0.975 for PMUBCs. The role of these boundary conditions in finite element analyses of whole bones and bone-implant systems will need to be investigated in future work.

Seismic investigations of the Greenland ice sheet at EGIG points

The compressional and the shear wave velocities in the Greenland ice sheet are derived from seismic records of the EGIG 1959. Further the variation of velocities in the firn and the dependance of Poisson's ratio from depth are determined. At Station Centrale, two P-waves are recorded from underground layers. Their velocities show that the ice basement consists of crystalline rocks. The P-wave velocities derived from reflections agree well with those obtained by refraction shooting. From this agreement results that the ice is ± homogenous and ± isotropic for Pwaves. The elastic constants for isotropic ice are calculated. Finally the temperature dependence of the velocities is discussed.