Search

We provide a fundamental insight into the microscopic mechanisms of the ageing processes. Using large scale molecular dynamics simulations of the prototypical ferroelectric material PbTiO3, we demonstrate that the experimentally observed ageing phenomena can be reproduced from intrinsic interactions of defect-dipoles related to dopant-vacancy associates, even in the absence of extrinsic effects. We show that variation of the dopant concentration modifies the material's hysteretic response. We identify a universal method to reduce loss and tune the electromechanical properties of inexpensive ceramics for efficient technologies., Comment: 6 pages, 3 figures

The role of defects in solids of mixed ionic-covalent bonds such as ferroelectric oxides is complex. Current understanding of defects on ferroelectric properties at the single-defect level remains mostly at the empirical level, and the detailed atomistic mechanisms for many defect-mediated polarization-switching processes have not been convincingly revealed quantum mechanically. We simulate the polarization-electric field (P-E) and strain-electric field ({$\epsilon$}-E) hysteresis loops for BaTiO3 in the presence of generic defect dipoles with large-scale molecular dynamics and provide a detailed atomistic picture of the defect dipole-enhanced electromechanical coupling. We develop a general first-principles-based atomistic model, enabling a quantitative understanding of the relationship between macroscopic ferroelectric properties and dipolar impurities of different orientations, concentrations, and dipole moments. We find that the collective orientation of dipolar defects relative to the external field is the key microscopic structure feature that strongly affects materials hardening/softening and electromechanical coupling. We show that a small concentration (0.1 at.%) of defect dipoles dramatically improves electromechanical response. This offers the opportunity to improve the performance of inexpensive polycrystalline ferroelectric ceramics through defect dipole engineering for a range of applications including piezoelectric sensors, actuators, and transducers.

Ferroelectric perovskite oxides possess a large electrocaloric (EC) effect, but usually at high temperatures near the ferroelectric/paraelectric phase transition temperature, which limits their potential application as next-generation solid-state cooling devices. We use classical molecular dynamics to study the electric field induced phase transitions and EC effect in PMN-PT (PbMg1/3Nb2/3O3-PbTiO3). We find that the maximum EC strength of PMN-PT occurs within the morphotropic phase boundary (MPB) region at 300 K. The large adiabatic temperature change is caused by easy rotation of polarization within the MPB region., Comment: 17 pages, 7 figures

The lithium-carbon binary system possesses a broad range of chemical compounds, which exhibit fascinating chemical bonding characteristics, which give rise to diverse and technologically important properties. While lithium carbides with various compositions have been studied or suggested previously, the crystal structures of these compounds are far from well understood. In this work, we present the first comprehensive survey of all ground state (GS) structures of lithium carbides over a broad range of thermodynamic conditions, using ab initio density functional theory (DFT) crystal structure searching methods. Thorough searches were performed for 29 stoichiometries ranging from ${\mathrm{Li}}_{12}\mathrm{C}$ to ${\mathrm{LiC}}_{12}$ at 0 and 40 GPa. Based on formation enthalpies from optimized van der Waals density functional calculations, three thermodynamically stable phases $({\mathrm{Li}}_{4}{\mathrm{C}}_{3}$, ${\mathrm{Li}}_{2}{\mathrm{C}}_{2}$, and ${\mathrm{LiC}}_{12})$ were identified at 0 GPa, and seven thermodynamically stable phases $({\mathrm{Li}}_{8}\mathrm{C}$, ${\mathrm{Li}}_{6}\mathrm{C}$, ${\mathrm{Li}}_{4}\mathrm{C}$, ${\mathrm{Li}}_{8}{\mathrm{C}}_{3}$, ${\mathrm{Li}}_{2}\mathrm{C}$, ${\mathrm{Li}}_{3}{\mathrm{C}}_{4}$, and ${\mathrm{Li}}_{2}{\mathrm{C}}_{3})$ were predicted at 40 GPa. A rich diversity of carbon bonding, including monomers, dimers, trimers, nanoribbons, sheets, and frameworks, was found within these structures, and the dimensionality of carbon connectivity existing within each phase increases with increasing carbon concentration. We find that the well-known composition ${\mathrm{LiC}}_{6}$ is actually a metastable one. We also find a unique coexistence of carbon monomers and dimers within the predicted thermodynamically stable phase ${\mathrm{Li}}_{8}{\mathrm{C}}_{3}$, and different widths of carbon nanoribbons coexist in a metastable phase of ${\mathrm{Li}}_{2}{\mathrm{C}}_{2}$ (Imm2). Interesting mixed ${\mathit{sp}}^{2}\text{\ensuremath{-}}{\mathit{sp}}^{3}$ carbon frameworks are predicted in metastable phases with composition ${\mathrm{LiC}}_{6}$.

Physical review / B 95(21), 214120 (2017). doi:10.1103/PhysRevB.95.214120, For decades, numerous attempts have been made to produce polar oxynitride perovskites, where some of the oxygen is replaced by nitrogen, but a polar ordered oxynitride has never been demonstrated. Caracas and Cohen [Appl. Phys. Lett. 91, 092902 (2007)] studied possible ordered polar oxynitrides within density-functional theory (DFT) and found a few candidates that were predicted to be insulating and at least metastable. $YSiO_{2}N$ stood out with huge predicted polarization and nonlinear optic coefficients. In this study, we demonstrate the synthesis of perovskite-structured $YSiO_{2}N$ by using a combination of a diamond-anvil cell and in situ laser-heating techniques. Subsequent in situ x-ray diffraction, second-harmonic generation, and Raman-scattering measurements confirm that it is polar and a strong nonlinear optical material, with structure and properties similar to those predicted by DFT., Published by APS, Woodbury, NY

We report an accurate study of interactions between Benzene molecules using variational quantum Monte Carlo (VMC) and diffusion quantum Monte Carlo (DMC) methods. We compare these results with density functional theory (DFT) using different van der Waals (vdW) functionals. In our QMC calculations, we use accurate correlated trial wave functions including three-body Jastrow factors, and backflow transformations. We consider two benzene molecules in the parallel displaced (PD) geometry, and find that by highly optimizing the wave function and introducing more dynamical correlation into the wave function, we compute the weak chemical binding energy between aromatic rings accurately. We find optimal VMC and DMC binding energies of -2.3(4) and -2.7(3) kcal/mol, respectively. The best estimate of the CCSD(T)/CBS limit is -2.65(2) kcal/mol [E. Miliordos et al, J. Phys. Chem. A 118, 7568 (2014)]. Our results indicate that QMC methods give chemical accuracy for weakly bound van der Waals molecular interactions, comparable to results from the best quantum chemistry methods., Comment: Accepted for publication in the Journal of Chemical Physics, Vol. 143, Issue 11, 2015

We use molecular dynamics with a first-principles based shell model potential to study pyroelec- tricity in lithium niobate. We find that the primary pyroelectric effect is dominant, and pyroelec- tricity can be understood simply from the anharmonic change in crystal structure with temperature and the Born effective charges on the ions. This opens a new experimental route to studying py- roelectricity, as candidate pyroelectric materials can be studied with X-ray diffraction as a function of temperature in conjunction with theoretical effective charges. We also predict an appreciable pressure effect on pyroelectricity, which could be used to optimize materials pyroelectricity, and the converse electrocaloric effect, peak as Tc is approached., Comment: 4 pages, 3 figures, 3 tables

We use the density functional perturbation theory to determine for the first time the pressure evolution of the Raman intensities for a mineral, the two high-pressure structures of MgSiO3 perovskite and post-perovskite. At high pressures, the Raman powder spectra reveals three main peaks for the perovskite structure and one main peak for the post-perovskite structure. Due to the large differences in the spectra of the two phases Raman spectroscopy can be used as a good experimental indication of the phase transition., Comment: 16 pages, submitted to Geophysical Research Letters

We compute the lattice-dynamical and thermal equation of state properties of ferromagnetic bcc iron using the first principles linear response linear-muffin-tin-orbital method in the generalized-gradient approximation. The calculated phonon dispersion and phonon density of states, both at ambient and high pressures, show good agreement with inelastic neutron scattering data. We find the free energy as a function of volume and temperature, including both electronic excitations and phonon contributions, and we have derived various thermodynamic properties at high pressure and temperature. The thermal equation of state at ambient temperature agrees well with diamond-anvil-cell measurements. We have performed detailed investigations on the behavior of various thermal equation of state parameters, such as the bulk modulus, the thermal expansivity, the Anderson-Gruneisen parameter, the Gruneisen ratio, and the heat capacity as function of temperature and pressure. A detailed comparison has been made with available experimental measurements, as well as results from similar theoretical studies on nonmagnetic bcc Tantalum., Comment: 34 pages, 13 figures, and 1 table

The elastic properties of the $B_1$-structured transition-metal nitrides and their carbide counterparts are studied using the {\it ab initio\} density functional perturbation theory. The linear response results of elastic constants are in excellent agreement with those obtained from numerical derivative methods, and are also consistent with measured data. We find the following trends: (1) Bulk moduli $B$ and tetragonal shear moduli $G^{\prime}=(C_{11}-C_{12})/2$, increase and lattice constants $a_{0}$ decrease rightward or downward on the Periodic Table for the metal component or if C is replaced by N; (2) The inequality $B > G^{\prime} > G > 0$ holds for $G=C_{44}$; (3) $G$ depends strongly on the number of valence electrons per unit cell ($Z_{V}$). From the fitted curve of $G$ as a function of $Z_{V}$, we can predict that MoN is unstable in $B_{1}$ structure, and transition-metal carbonitrides ($e.g.$ ZrC$_{x}$N$_{1-x}$) and di-transition-metal carbides ($e.g.$ Hf$_{x}$Ta$_{1-x}$C) have maximum $G$ at $Z_{V} \approx 8.3$., 4 pages, 2 figures, submitted to PRL. 2 typos in ref. 15 were corrected

Marcus et al. (Marcus P, Ma H and Qiu S L 2002 J. Phys.: Condens. Matter 14 L525) claim that thermodynamic properties of materials under pressure must be computed using the Gibbs free energy $G$, rather than the internal energy $E$. Marcus et al. state that ``The minima of $G$, but not of $E$, give the equilibrium structure; the second derivatives of $G$, but not of $E$, with respect to strains at the equilibrium structure give the equilibrium elastic constants.'' Both statements are incorrect., Comment: Comment