Researcher profile

R. E. Cohen

R. E. Cohen contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
30works
0followers
12topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

30 published item(s)

preprint2022arXiv

First-principles Calculations of Raman and Infrared Spectroscopy For Phase Identification and Strain Calibration of Hafnia

Using density functional perturbation theory (DFPT) we computed the phonon frequencies, Raman and IR activities of hafnia polymorphs (P4$_{2}$nmc, Pca2$_{1}$, Pmn2$_{1}$, Pbca OI, brookite, and baddeleyite) for phase identification. We investigated the evolution of Raman and IR activities with respect to epitaxial strain and provide plots of frequency differences as a function of strain for experimental calibration and identification of the strain state of the sample. We found Raman signatures of different hafnia polymorphs: $ω(A_{1g})=300$ cm$^{-1}$ for P4$_{2}$nmc, $ω(A_{1})=343$ cm$^{-1}$ for Pca2$_{1}$, $ω(B_{2})=693$ cm$^{-1}$ for Pmn2$_{1}$, $ω(A_{g})=513$ cm$^{-1}$ for Pbca (OI), $ω(A_{g})=384$ cm$^{-1}$ for brookite, and $ω(A_{g}) = 496$ cm$^{-1}$ for baddeleyite. We also identified the Raman $B_{1g}$ mode, an anti-phase vibration of dipole moments, ( $ω(B_{1g}) = 758$ cm$^{-1}$ for OI, $ω(B_{1g})= 784$ cm$^{-1}$ for brookite) as the Raman signature of antipolar Pbca structures. We calculated a large splitting between longitudinal optical (LO) and transverse optical (TO) modes ($Δ{ω_{\text{LO-TO}}(A^{z}_{1})}=255$ cm$^{-1}$ in Pca2$_{1}$, and $Δ{ω_{\text{LO-TO}}(A_{1})}=263$ cm$^{-1}$ in Pmn2$_{1}$) to the same order as those observed in perovskite ferroelectrics, and related them to the anomalously large Born effective charges of Hf atoms ($Z^{*}(\text{Hf}) = 5.54$).

preprint2022arXiv

Origin of Ferroelectricity in Hafnia from Epitaxial Strain

Ferroelectric hafnia is being explored for next generation electronics due to its robust ferroelectricity in nanoscale samples and its compatibility with silicon. However, its ferroelectricity is not understood. Other ferroelectrics usually lose their ferroelectricity for nanoscopic samples and thin films, and the hafnia ground state is non-polar baddeleyite. Here we study hafnia with density functional theory (DFT) under epitaxial strain, and find that strain not only stabilizes the ferroelectric phases, but also leads to unstable modes and a downhill path in energy from the high temperature tetragonal structure. We find that for tensile epitaxial strains corresponding to a square substrate of lattice constant $a_{epi} \geqslant 5.38$ Å~the ferroelectric \oI is most stable, even more stable than baddeleyite. Furthermore, we find that under tensile epitaxial strain $η$ the tetragonal phase will distort to one of the two ferroelectric phases: for $η> 1.5$\%, the $Γ^{-}_{5}$ mode is unstable and leads to \oII , and at $η> 3.75$\% coupling between this mode and the zone boundary M1 mode leads to \oI.

preprint2020arXiv

Prediction of an Extended Ferroelectric Clathrate

Using first-principles calculations, we predict a lightweight room-temperature ferroelectric carbon-boron framework in a host/guest clathrate structure. This ferroelectric clathrate, with composition ScB$_3$C$_3$, exhibits high polarization density and low mass density compared with widely used commercial ferroelectrics. Molecular dynamics simulations show spontaneous polarization with a moderate above-room-temperature T$_c$ of $\sim$370 K, which implies large susceptibility and possibly large electrocaloric and piezoelectric constants at room temperature. Our findings open the possibility for a new class of ferroelectric materials with potential across a broad range of applications.

preprint2019arXiv

The vertical Na-O relation in the bulge globular cluster NGC 6553

In this article, we present a detailed chemical analysis of seven red giant members of NGC 6553 using high-resolution spectroscopy from VLT FLAMES. We obtained the stellar parameters (Teff, Log(g), vt, [Fe/H]) of these stars from the spectra, and we measured the chemical abundance for 20 elements, including light elements, iron-peak elements, alpha-elements and neutron-capture elements. The metallicities in our sample stars are consistent with a homogeneous distribution. We found a mean of [Fe/H]=-0.14+/-0.07 dex, in agreement with other studies. Using the alpha-elements Mg, Si, Ca and Ti we obtain the mean of [alpha/Fe]=0.11+/-0.05. We found a vertical relation between Na and O, characterized by a significant spread in Na and an almost non-existent spread in O. In fact, Na and Al are the only two light elements with a large intrinsic spread, which demonstrates the presence of Multiple Populations (MPs). An intrinsic spread in Mg is not detected in this study. The alpha, iron-peak and neutron capture elements show good agreement with the trend of the bulge field stars, indicating similar origin and evolution, in concordance with our previous studies for two other bulge GCs (NGC 6440 and NGC 6528).

preprint2016arXiv

Low-pressure phase diagram of crystalline benzene from quantum Monte Carlo

We study the low-pressure (0 to 10 GPa) phase diagram of crystalline benzene using quantum Monte Carlo (QMC) and density functional theory (DFT) methods. We consider the $Pbca$, $P4_32_12$, and $P2_1/c$ structures as the best candidates for phase I and phase II. We perform diffusion quantum Monte Carlo (DMC) calculations to obtain accurate static phase diagrams as benchmarks for modern van der Waals density functionals. We use density functional perturbation theory to compute phonon contribution in the free-energy calculations. Our DFT enthalpy-pressure phase diagram indicates that the $Pbca$ and $P2_1/c$ structures are the most stable phases within the studied pressure range. The DMC Gibbs free-energy calculations predict that the room temperature $Pbca$ to $P2_1/c$ phase transition occurs at 2.1(1) GPa. This prediction is consistent with available experimental results at room temperature. Our DMC calculations show an estimate of 50.6$\pm$0.5 kJ/mol for crystalline benzene lattice energy.

preprint2016arXiv

Structural Diversity in Lithium Carbides

The lithium-carbon binary system possesses a broad range of chemical compounds, which exhibit fascinating chemical bonding characteristics that give rise 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 \textit{ab initio} density functional theory (DFT) crystal structure searching methods. Thorough searches were performed for 29 stoichiometries ranging from Li$_{12}$C to LiC$_{12}$ at 0 GPa as well as 40 GPa. Based on formation enthalpies from optimized van der Waals density functional calculations, three thermodynamically stable phases (Li$_4$C$_3$, Li$_2$C$_2$ and LiC$_{12}$) were identified at 0 GPa, and seven thermodynamically stable phases (Li$_8$C, Li$_6$C, Li$_4$C, Li$_8$C$_3$, Li$_2$C, Li$_3$C$_4$, and Li$_2$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 was observed to increase with increasing carbon concentration. Of particular interest, we find that the well-known composition LiC$_6$ is actually a metastable one. We also find a unique coexistence of carbon monomers and dimers within the predicted thermodynamically stable phase Li$_8$C$_3$, and different widths of carbon nanoribbons coexist in a metastable phase of Li$_2$C$_2$ (\textit{Imm}2). Interesting mixed \textit{sp}$^2$-\textit{sp}$^3$ carbon frameworks are predicted in metastable phases with composition LiC$_6$.

preprint2016arXiv

Ultra-deep GEMINI near-infrared observations of the bulge globular cluster NGC 6624

We used ultra-deep $J$ and $K_s$ images secured with the near-infrared GSAOI camera assisted by the multi-conjugate adaptive optics system GeMS at the GEMINI South Telescope in Chile, to obtain a ($K_s$, $J-K_s$) color-magnitude diagram (CMD) for the bulge globular cluster NGC 6624. We obtained the deepest and most accurate near-infrared CMD from the ground for this cluster, by reaching $K_s$ $\sim$ 21.5, approximately 8 magnitudes below the horizontal branch level. The entire extension of the Main Sequence (MS) is nicely sampled and at $K_s$ $\sim$ 20 we detected the so-called MS "knee" in a purely near-infrared CMD. By taking advantage of the exquisite quality of the data, we estimated the absolute age of NGC 6624 ($t_{age}$ = 12.0 $\pm$ 0.5 Gyr), which turns out to be in good agreement with previous studies in the literature. We also analyzed the luminosity and mass functions of MS stars down to M $\sim$ 0.45 M$_{\odot}$ finding evidence of a significant increase of low-mass stars at increasing distances from the cluster center. This is a clear signature of mass segregation, confirming that NGC 6624 is in an advanced stage of dynamical evolution.

preprint2015arXiv

Chemical accuracy from quantum Monte Carlo for the Benzene Dimer

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.

preprint2014arXiv

Effects of Manganese Addition on the Electronic Structure of BaTiO$_3$

Mn is used as a dopant to improve the electromechanical properties of perovskite oxides. We elucidate the role of Mn defects and associated vacancies on the electronic, atomic and ferroelectric properties of BaTiO$_3$. Using density functional theory (DFT) and DFT+U we investigate the equilibrium geometry and electronic properties of Mn ion on A or B sites and with compensating oxygen vacancies. We study the change in the oxidation state of Mn in response to local environment changes, such as the presence of oxygen vacancies.

preprint2014arXiv

Equations of state and stability of MgSiO$_3$ perovskite and post-perovskite phases from quantum Monte Carlo simulations

We have performed quantum Monte Carlo (QMC) simulations and density functional theory (DFT) calculations to study the equations of state of MgSiO$_3$ perovskite (Pv) and post-perovskite (PPv), up to the pressure and temperature conditions of the base of Earth's lower mantle. The ground state energies were derived using QMC and the temperature dependent Helmholtz free energies were calculated within the quasi-harmonic approximation and density functional perturbation theory. The equations of state for both phases of MgSiO$_3$ agree well with experiments, and better than those from generalized gradient approximation (GGA) calculations. The Pv-PPv phase boundary calculated from our QMC equations of states is also consistent with experiments, and better than previous LDA calculations. We discuss the implications for double crossing of the Pv-PPv boundary in the Earth.

preprint2014arXiv

Prediction of a new potential high-pressure structure of FeSiO$_3$

We predict a new candidate high-temperature high-pressure structure of FeSiO$_3$ with space-group symmetry Cmmm by applying an evolutionary algorithm within DFT+U that we call post-perovskite II (PPv-II). An exhaustive search found no other competitive candidate structures with ABO$_3$ composition. We compared the X-ray diffraction (XRD) pattern of FeSiO$_3$ PPv-II with experimental results of the recently reported H-phase of (Fe,Mg)SiO$_3$. The intensities and positions of two main X-ray diffraction peaks of PPv-II FeSiO$_3$ compare well with those of the H-phase. We also calculated the static equation of state, the enthalpy and the bulk modulus of the PPv-II phase and compared it with those of perovskite (Pv) and post-perovskite (PPv) phases of FeSiO$_3$. According to the static DFT+U computations the PPv-II phase of FeSiO$_3$ is less stable than Pv and PPv phases under lower mantle pressure conditions at 0 K and has a higher volume. PPv-II may be entropically stabilized, and may be a stable phase in Earth$'$s lower mantle, coexisting with $α$-PbO$_2$ (Columbite-structured) silica and perovskite, or with magnesiowustite or ferropericlase, depending on bulk composition.

preprint2014arXiv

Pressure Suppression of Electron Correlation in the Collapsed Tetragonal Phase of $CaFe_2As_2$ : A DFT-DMFT Investigation

Recent studies reveal a pressure induced transition from a paramagnetic tetragonal phase (T) to a collapsed tetragonal phase (CT) in $CaFe_2As_2$, which was found to be superconducting with non-hydrostatic pressure at low temperature. We have investigated the effects of electron correlation and local fluctuating moment in both tetragonal and collapsed tetragonal phases of the paramagnetic $CaFe_2As_2$ using self consistent DFT+DMFT with continuous time quantum Monte Carlo as the impurity solver. From the computed optical conductivity, we find a gain in the optical kinetic energy due to the loss in Hund's rule coupling energy in the CT. We find that the transition from T to CT turns $CaFe_2As_2$ from a bad metal into a good metal. Computed mass enhancement and local moments also show significant decrease in the CT, which confirms the suppression of the electron correlation in the CT phase of $CaFe_2As_2$.

preprint2014arXiv

Strong Pressure Dependent Electron-Phonon Coupling in FeSe

We have computed the correlated electronic structure of FeSe and its dependence on the A1g mode versus compression. Using the self-consistent density functional theory - dynamical mean field theory (DFT-DMFT) with continuous time quantum Monte Carlo (CTQMC), we find that there is greatly enhanced coupling between some correlated electron states and the A1g lattice distortion. Superconductivity in FeSe shows a very strong sensitivity to pressure, with an increase in Tc of almost a factor of 5 within a few GPa, followed by a drop, despite monotonic pressure dependence of almost all electronic properties. We find that the maximum A1g deformation potential behaves similar to the experimental Tc. In contrast, the maximum deformation potential in DFT for this mode increases monotonically with increasing pressure.

preprint2013arXiv

Graphene physics and insulator-metal transition in compressed hydrogen

Compressed hydrogen passes through a series of layered structures in which the layers can be viewed as distorted graphene sheets. The electronic structures of these layered structures can be understood by studying simple model systems- an ideal single hydrogen graphene sheet and three-dimensional model lattices consisting of such sheets. The energetically stable structures result from structural distortions of model graphene-based systems due to electronic instabilities towards Peierls or other distortions associated with the opening of a band gap. Two factors play crucial roles in the metallization of compressed hydrogen: (i) crossing of conduction and valence bands in hexagonal or graphene-like layers due to topology and (ii) formation of bonding states with $2p_z$ $π$ character.

preprint2012arXiv

Giant electrocaloric effect around T$_c$

We use molecular dynamics with a first-principles-based shell model potential to study the electrocaloric effect (ECE) in lithium niobate, LiNbO$_3$, and find a giant electrocaloric effect along a line passing through the ferroelectric transition. With applied electric field, a line of maximum ECE passes through the zero field ferroelectric transition, continuing along a Widom line at high temperatures with increasing field, and along the instability that leads to homogeneous ferroelectric switching below $T_c$ with an applied field antiparallel to the spontaneous polarization. This line is defined as the minimum in the inverse capacitance under applied electric field. We investigate the effects of pressure, temperature and applied electric field on the ECE. The behavior we observe in LiNbO$_3$ should generally apply to ferroelectrics; we therefore suggest that the operating temperature for refrigeration and energy scavenging applications should be above the ferroelectric transition region to obtain large electrocaloric response. We find a relationship among $T_c$, the Widom line and homogeneous switching that should be universal among ferroelectrics, relaxors, multiferroics, and the same behavior should be found under applied magnetic fields in ferromagnets.

preprint2012arXiv

Pressure dependence of the monoclinic phase in (1-x)Pb(Mg1/3Nb2/3)O3-xPbTiO3 solid solutions

We combine high-pressure x-ray diffraction, high-pressure Raman scattering, and optical microscope to investigate a series of PMN-xPT solid solutions (x=0.2, 0.3, 0.33, 0.35, 0.37, 0.4) in diamond anvil cells up to 20 GPa at 300 K. The Raman spectra show a peak centered at 380 cm-1 above 6 GPa for all samples, in agreement with previous observations. X-ray diffraction measurements are consistent with this spectral change indicating structural phase transition; we find that the triplet at the pseudo cubic (220) Bragg peak merges into a doublet above 6 GPa. Our results indicate that the morphotropic phase boundary region (x=0.33 to 0.37) with the presence of monoclinic symmetry persists up to 7 GPa. The pressure dependence of ferroelectric domains in PMN-0.32PT single crystals was observed using a polarizing optical microscope. The domain wall density decreases with pressure and domains disappear at a modest pressure of 3 GPa. We propose a pressure-composition phase diagram for PMN-xPT solid solutions.

preprint2011arXiv

Experimental and theoretical evidence for pressure-induced metallization in FeO with the rock-salt type structure

Electrical conductivity of FeO was measured up to 141 GPa and 2480 K in a laserheated diamond-anvil cell. The results show that rock-salt (B1) type structured FeO metallizes at around 70 GPa and 1900 K without any structural phase transition. We computed fully self-consistently the electronic structure and the electrical conductivity of B1 FeO as a function of pressure and temperature, and found that although insulating as expected at ambient condition, B1 FeO metallizes at high temperatures, consistent with experiments. The observed metallization is related to spin crossover.

preprint2011arXiv

First principles based atomistic modeling of phase stability in PMN-xPT

We have performed molecular dynamics simulations using a shell model potential developed by fitting first principles results to describe the behavior of the relaxor-ferroelectric (1-x)PbMg1/3Nb2/3O3-xPbTiO3 (PMN-xPT) as function of concentration and temperature, using site occupancies within the random site model. In our simulations, PMN is cubic at all temperatures and behaves as a polar glass. As a small amount of Ti is added, a weak polar state develops, but structural disorder dominates, and the symmetry is rhombohedral. As more Ti is added the ground state is clearly polar and the system is ferroelectric, but with easy rotation of the polarization direction. In the high Ti content region, the solid solution adopts ferroelectric behavior similar to PT, with tetragonal symmetry. The ground state sequence with increasing Ti content is R-MB-O-MC-T. The high temperature phase is cubic at all compositions. Our simulations give the slope of the morphotropic phase boundaries, crucial for high temperature applications. We find that the phase diagram PMN-xPT can be understood within the random site model.

preprint2011arXiv

Orbital ordering, ferroelasticity, and the large pressure induced volume collapse in PbCrO3

We report a new tetragonal ground-state for perovskite-structured PbCrO3 from DFT+U calculations, and explain its anomalously large volume. The new structure is stabilized due to orbital ordering of Cr-d in the presence of a large tetragonal crystal field, mainly due to off-centering of the Pb atom. At higher pressures (smaller volumes) there is a first-order transition to a cubic phase where the Cr-d orbitals are orbitally liquid. This phase-transition is accompanied by a ~11.5% volume collapse, one of the largest known for transition-metal oxides. The large ferroelasticity and its strong coupling to the orbital degrees of freedom could be exploited to form potentially useful magnetostrictive materials

preprint2010arXiv

Finite-temperature magnetism in bcc Fe under compression

We investigate the contributions of finite-temperature magnetic fluctuations to the thermodynamic properties of bcc Fe as a function of pressure. First, we apply a tight-binding total-energy model parameterized to first-principles linearized augmented planewave computations to examine various ferromagnetic, anti-ferromagnetic, and noncollinear spin spiral states at zero temperature. The tight-binding data are fit to a generalized Heisenberg Hamiltonian to describe the magnetic energy functional based on local moments. We then use Monte Carlo simulations to compute the magnetic susceptibility, the Curie temperature, heat capacity, and magnetic free energy. Including the finite-temperature magnetism improves the agreement with experiment for the calculated thermal expansion coefficients.

preprint2010arXiv

First-principles studies on electrical resistivity of iron under pressure

We investigate the temperature and pressure dependences of the electrical resistivity for bcc and hcp Fe using the low-order variational approximation and theoretical transport spectral functions calculated from first principles linear response linear-muffin-tin-orbital method in the generalized-gradient approximation. The calculated values are in close agreement with available experimental data, and show strong increase with temperature and decrease with pressure. We also discuss the behavior of the electrical resistivity for the bcc -> hcp phase transition.

preprint2010arXiv

Fundamental high pressure calibration from all-electron quantum Monte Carlo calculations

We develop an all-electron quantum Monte Carlo (QMC) method for solids that does not rely on pseudopotentials, and use it to construct a primary ultra-high pressure calibration based the equation of state of cubic boron nitride(c-BN). We compute the static contribution to the free energy with QMC, and obtain the phonon contribution from density functional theory, yielding a high-accuracy calibration up to 900 GPa usable directly in experiment. Furthermore, we compute the anharmonic Raman frequency shift with QMC as a function of pressure and temperature, allowing optical pressure calibration in table-top experiments. In contrast to present experimental approaches, small systematic errors in the theoretical EOS do not increase with pressure, and no extrapolation is needed. This all-electron methodology is generally applicable to first-row solids, and can be used to provide a new reference for ab initio calculations of solids and to benchmark pseudopotential accuracy.

preprint2010arXiv

Origin of Pyroelectricity in LiNbO3

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.

preprint2010arXiv

Quantum Monte Carlo for minerals at high pressure: Phase stability, equations of state, and elasticity of silica

Silica is an abundant component of the Earth whose crystalline polymorphs play key roles in its structure and dynamics. As the simplest silicates, understanding pure silica is a prerequisite to understanding the rocky part of the Earth, its majority. First principle density functional theory (DFT) methods have often been used to accurately predict properties of silicates. Here, we study silica with quantum Monte Carlo (QMC), which until now was not computationally possible for such complex materials, and find that QMC overcomes the failures of DFT. QMC is a benchmark method that does not rely on density functionals, but rather explicitly treats the electrons and their interactions via a stochastic solution of Schrodinger's equation. Using ground state QMC plus phonons within the quasiharmonic approximation from density functional perturbation theory, we obtain the thermal pressure and equations of state of silica phases up to Earth's core-mantle boundary. Our results provide the most well-constrained equations of state and phase boundaries available for silica. QMC indicates a transition to the most dense alpha-PbO2 structure above the core-insulating D" layer, suggesting the absence of significant free silica in the bulk lower mantle, which has been assumed but never proven. We also find an accurate shear elastic constant and its geophysically important softening with pressure.

preprint2009arXiv

Elastic isotropy of hcp-Fe under Earth core conditions

Our first-principles calculations show that both the compressional and shear waves of hcp-Fe become elastically isotropic under the high temperatures of Earth inner core conditions, with the variation in sound velocities along different angles from the c axis within 1%. We computed the thermoelasticity at high pressures and temperatures from quasiharmonic linear response linear-muffin-tin-orbital calculations in the generalized-gradient approximation. The calculated anisotropic shape and magnitude in hcp-Fe at ambient temperature agree well with previous first-principles predictions, and the anisotropic effects show strong temperature dependences. This implies that other mechanisms, rather than the preferential alignment of the hcp-Fe crystal along the Earth rotation axis, account for the seismic P-wave travel time anomalies. Either the inner core is not hcp iron, and/or the seismologically observed anisotropy is caused by inhomogeneity, i.e. multiple phases.

preprint2009arXiv

First-principles thermal equation of state and thermoelasticity of hcp Fe at high pressures

We investigate the equation of state and elastic properties of hcp iron at high pressures and high temperatures using first principles linear response linear-muffin-tin-orbital method in the generalized-gradient approximation. We calculate the Helmholtz free energy as a function of volume, temperature, and volume-conserving strains, including the electronic excitation contributions from band structures and lattice vibrational contributions from quasi-harmonic lattice dynamics. We perform detailed investigations on the behavior of elastic moduli and equation of state properties as functions of temperature and pressure, including the pressure-volume equation of state, bulk modulus, the thermal expansion coefficient, the Gruneisen ratio, and the shock Hugoniot. Detailed comparison has been made with available experimental measurements and theoretical predictions.

preprint2009arXiv

Origin of Diffuse Scattering in Relaxor Ferroelectrics

High-pressure and variable temperature single crystal synchrotron X-ray measurements combined with first-principles based molecular dynamics simulations study diffuse scattering in the relaxor ferroelectric system PSN (PbSc$_{1/2}$Nb$_{1/2}$O$_3$). Constant temperature experiments show pressure induced transition to the relaxor phase at different temperatures characterized by butterfly and rod shaped diffuse scattering around the $\{$h00$\}$ and $\{$hh0$\}$ Bragg spots, respectively. The simulations reproduce the observed diffuse scattering features as well as their pressure-temperature behavior, and show that they arise from polarization correlations between chemically-ordered regions, which in previous simulations were shown to behave as polar nanoregions. Simulations also exhibit radial diffuse scattering (elongated towards and away from {\bf Q}=(000)), that persists even in the paraelectric phase, consistent with previous neutron experiments on (PbMg$_{1/3}$Nb$_{2/3}$O$_3$) (PMN)

preprint2007arXiv

First Principles Theories of Piezoelectric Materials

Piezoelectrics have long been studied using parameterized models fit to experimental data, starting with the work of Devonshire in 1954. Much has been learned using such approaches, but they can also miss major phenomena if the materials properties are not well under-stood, as is exemplified by the realization that low-symmetry monoclinic phases are common around morphotropic phase boundaries, which was missed completed by low-order Devonshire models, and can only appear in higher order models. In the last 15 years, a new approach has developed using first-principles computations, based on fundamental physics, with no essential experimental input other than the desired chemistry (nuclear charges). First-principles theory laid the framework for a basic understanding of the origins of ferroelectric behavior and piezoelectric properties. The range of properties accessible to theory continues to expand as does the accuracy of the predictions. We are moving towards the ability to design materials of desired properties computationally. Here some of the fundamental developments of our understanding of piezoelectric material behavior and ability to predict a wide range of properties using theoretical methods are reviewed.

preprint2003arXiv

A Multi-Year Photometric Study of IC 348

The extremely young cluster IC 348 has been monitored photometrically over 5 observing seasons from Dec 1998 to March 2003 in Cousins I with a 0.6 m telescope at Van Vleck Observatory. Twenty-eight periodic variables and 16 irregular variables have been identified. Among the brighter stars, 14 of the 16 known K or M-type WTTS were found to be periodic variables, while all 5 of the known CTTS were found to be irregular variables. In the full sample, which includes 150 stars with I mag as faint as 18, we find that 40% of the 63 WTTS are detected as variables, nearly all of them periodic, while 55% of the 20 CTTS are also detected as variable, with none of them periodic. Our study suggests that 80-90% of all WTTS in young clusters will be detected as periodic variables given sufficiently precise and extended monitoring, whereas CTTS will reveal themselves primarily or solely as irregular variables. This has clear consequences for PMS rotational studies based on photometric periods. We examine the stability of the periodic light curves from season to season. All periodic stars show modulations of their amplitude, mean brightness and light curve shape on time scales of less than 1 yr, presumably due to changes in spot configurations and/or physical characteristics. In no case, however, can we find definitive evidence of a change in period, indicating that differential rotation is probably much less on WTTS than it is on the Sun. Among the non-periodic stars, we report the detection of two possible UXors as well as a pre-main sequence star, HMW 15, which apparently undergoes an eclipse with a duration exceeding three years.

preprint1995arXiv

Periodic Slab LAPW Computations for Ferroelectric BaTiO$_3$

Linearized augmented plane wave (LAPW) calculations are performed for periodic (001) and (111) slabs of BaTiO$_3$ to understand the effects of surfaces on ferroelectric BaTiO$_3$. The (111) slab is found to be much less stable than the (001) slab. The average surface energies are respectively 3700 erg/cm$^2$ and 1600 erg/cm$^2$. The depolarization field is sufficiently large in the ideal unrelaxated slab to completely inhibit the ferroelectric instability. No mid-gap surface states are evident, but there are surface states in the upper gap in the unrelaxed slab and at the top of the valence band. The dangling surface Ti bonds self-heal making the Ti-O surface very reactive and an excellent epitaxial substrate. The charge density on atoms only one unit cell away from the surface are almost identical to the bulk. Keywords: ferroelectric, surface, slab, electronic structure, depolarization, BaTiO$_3$, thin film