BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Renata M. Wentzcovitch

Renata M. Wentzcovitch contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.mtrl-sciphysics.geo-phcond-mat.str-elcond-mat.otherphysics.comp-ph

Trust snapshot

Quick read

Trust 21 - Emerging
21works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

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

21 published item(s)

preprint2025arXiv

Ab initio superionic-liquid phase diagram of Fe1-xOx under Earth's inner core conditions

The superionic state is a phase of matter in which liquid-like ionic mobility coexists with a solid crystalline lattice. Recently identified in Earth's inner core (IC), this state has attracted considerable attention for its unique kinetic behavior and geophysical implications. However, the ab initio phase diagram describing the equilibrium between the superionic phase and the liquid solution under core conditions remains largely unexplored. Here, we present a thermodynamic approach to compute the Gibbs free energy and construct the ab initio superionic-liquid phase diagram for the Fe1-xOx system under IC conditions. We find that oxygen forms superionic states in both hcp and bcc Fe phases, with a pronounced influence on cooperative diffusion of iron in the bcc lattice. The stability fields of these superionic phases are sensitive to oxygen stoichiometry. The presence of superionic states leads to a higher oxygen concentration in the IC than previously estimated. Our work establishes a framework for investigating superionic-liquid equilibria under extreme conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A deep machine learning potential for atomistic simulation of Fe-Si-O systems under Earth's outer core conditions

Using artificial neural-network machine learning (ANN-ML) to generate interatomic potentials has been demonstrated to be a promising approach to address the long-standing challenge of accuracy versus efficiency in molecular dynamics (MD) simulations. Here, taking the Fe-Si-O system as a prototype, we show that accurate and transferable ANN-ML potentials can be developed for reliable MD simulations of materials at high-pressure and high-temperature conditions of the Earth's outer core. The ANN-ML potential for Fe-Si-O system is trained by fitting to the energies and forces of related binaries and ternary liquid structures at high pressures and temperatures obtained by first-principles calculations based on density functional theory (DFT). We show that the generated ANN-ML potential describes well the structure and dynamics of liquid phases of this complex system. The efficient ANN-ML potential with DFT accuracy provides a promising scheme for accurate atomistic simulations of structures and dynamics of complex Fe-Si-O system in the Earth's outer core.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Ab initio investigation of H-bond disordering in $δ$-AlOOH

$δ$-AlOOH ($δ$) is a high-pressure hydrous phase that participates in the deep geological water cycle. At 0 GPa, $δ$ has asymmetric hydrogen bonds (H-bonds). Under pressure, it exhibits H-bond disordering, tunneling, and finally, H-bond symmetrization at ~18 GPa. This study investigates these 300 K pressure-induced state changes in $δ$ with ab initio calculations. H-bond disordering in $δ$ was modeled using supercell multi-configuration quasiharmonic calculations. We examine: (a) energy barriers for proton jumps, (b) the pressure dependence of phonon frequencies, (c) 300 K compressibility, (d) neutron diffraction pattern anomalies, and (e) compare ab initio bond lengths with measured ones. Such thorough and systematic comparisons indicate that: (a) proton "disorder" has a restricted meaning when applied to $δ$. Nevertheless, H-bonds are disordered between 0 and 8 GPa, and a gradual change in H-bond configuration results in enhanced compressibility. (b) several structural and vibrational anomalies at ~8 GPa are consistent with the disappearance of a particular (HOC-12) H-bond configuration and its change into another one (HOC-11*). (c) between 8-11 GPa, H-bond configuration (HOC-11*) is generally ordered, at least in short- to mid-range scale. (d) between 11.5-18 GPa, H-bond lengths approach a critical value that impedes compression, resulting in decreased compressibility. In this pressure range, especially approaching H-bond symmetrization at ~18 GPa, anharmonicity and tunneling should play an essential role in the proton dynamics. Further simulations accounting for these effects are desirable to clarify the protons' state in this pressure range.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Structure and motifs of iron oxides from 1 to 3 TPa

Iron oxides are fundamental components of planet-forming materials. Understanding the Fe-O system's behavior and properties under high pressure can help us identify many new phases and states possible in exoplanetary interiors, especially terrestrial ones. Using the adaptive genetic algorithm (AGA), we investigate the structure of iron oxides for a wide range of stoichiometries ($0.25\leq x_O \leq 0.8$) at 1, 2, and 3 TPa. Five unreported ground-state structures with Fe$_2$O, FeO, Fe$_3$O$_5$, FeO$_2$, and FeO$_4$ compositions are identified. The calculated density of states (DOS) suggests that, except for FeO$_4$, all phases are metallic, but their carrier densities decrease with increasing pressure and oxygen content. The cluster alignment analysis of stable and metastable phases shows that several motifs may co-exist in a structure of iron oxides with low O content. In contrast, most iron oxides with high O content adopt a simple BCC motif at TPa pressures. Our results provide a crystal structure database of iron oxides for modeling and understanding the interiors of exoplanets.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Thermodynamics of spin crossover in ferropericlase: an improved LDA+$U_{sc}$ calculation

We present LDA+$U_{sc}$ calculations of high-spin (HS) and low-spin (LS) states in ferropericlase (fp) with an iron concentration of 18.75$\%$. The Hubbard parameter $U$ is determined self-consistently with structures optimized at arbitrary pressures. We confirm a strong dependence of $U$ on the pressure and spin state. Static calculations confirm that the antiferromagnetic configuration is more stable than the ferromagnetic one in the HS state, consistent with low-temperature measurements. Phonon calculations guarantee the dynamical stability of HS and LS states throughout the pressure range of the Earth mantle. Compression curves for HS and LS states agree well with experiments. Using a non-ideal mixing model for the HS to LS states solid solution, we obtain a crossover starting at $\sim$45 GPa at room temperature and considerably broader than previous results. The spin-crossover phase diagram is calculated, including vibrational, magnetic, electronic, and non-ideal HS-LS entropic contributions. Our results suggest the mixed-spin state predominates in fp in most of the lower mantle.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Two-step nucleation of the Earth's inner core

It has long been assumed the Earth's solid inner core started to grow when molten iron cooled to its melting point. However, the nucleation mechanism, which is a necessary step of crystallization, has not been well understood. Recent studies found it requires an unrealistic degree of undercooling to nucleate the stable hexagonal close-packed (hcp) phase of iron, which can never be reached under the actual Earth's core conditions. This contradiction leads to the inner core nucleation paradox [1]. Here, using a persistent-embryo method and molecular dynamics simulations, we demonstrate that the metastable body-centered cubic (bcc) phase of iron has a much higher nucleation rate than the hcp phase under inner-core conditions. Thus, the bcc nucleation is likely to be the first step of inner core formation instead of direct nucleation of the hcp phase. This mechanism reduces the required undercooling of iron nucleation, which provides a key factor to solve the inner-core nucleation paradox. The two-step nucleation scenario of the inner core also opens a new avenue for understanding the structure and anisotropy of the present inner core.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The Phonon Quasiparticle Approach for Anharmonic Properties of Solids

Knowledge of lattice anharmonicity is essential to elucidate distinctive thermal properties in crystalline solids. Yet, accurate \textit{ab initio} investigations of lattice anharmonicity encounter difficulties owing to the cumbersome computations. Here we introduce the phonon quasiparticle approach and review its application to various materials. This method efficiently and reliably addresses lattice anharmonicity by combining \textit{ab initio} molecular dynamics and lattice dynamics calculations. Thus, in principle, it accounts for full anharmonic effects and overcomes finite-size effects typical of \textit{ab initio} molecular dynamics. The validity and effectiveness of the current approach are demonstrated in the computation of thermodynamic and heat transport properties of weakly and strongly anharmonic systems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

LDA+Usc calculations of phase relations in FeO

Using the LDA+$\textit{U}_\text{sc}$ method, we present calculations phase relations of iron monoxides involving five polytypes in multiple spin-state configurations. The Hubbard parameter $U$ is determined self-consistently simultaneously with the occupation matrix and structures at arbitrary pressures. The Hubbard parameter strongly depends on pressure, structure, and spin state. Comparison with experimental structural data indicates the LDA+$\textit{U}_\text{sc}$ can predict structure, compression curves, phase relations, and transition pressures very well for the insulating B1 and iB8 states. However, it requires additional calculations using the Mermin functional that includes the electronic entropic contribution to the free energy to obtain an nB8 metallic state and a consistent iB8 to nB8 insulator-to-metal transition pressure.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A New Line Defect in NdTiO3 Perovskite

Perovskite oxides form an eclectic class of materials owing to their structural flexibility in accommodating cations of different sizes and valences. They host well known point and planar defects, but so far no line defect has been identified other than dislocations. Using analytical scanning transmission electron microscopy (STEM) and ab initio calculations we have detected and characterized the atomic and electronic structures of a novel line defect in NdTiO3 perovskite. It appears in STEM images as a perovskite cell rotated by 45 degrees. It consists of self-organized Ti-O vacancy lines replaced by Nd columns surrounding a central Ti-O octahedral chain containing Ti4+ ions, as opposed to Ti3+ in the host. The distinct Ti valence in this line defect introduces the possibility of engineering exotic conducting properties in a single preferred direction and tailoring novel desirable functionalities in this Mott insulator.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Spin crossover in (Mg,Fe$^{3+}$)(Si,Fe$^{3+}$)O$_3$ bridgmanite: effects of disorder, iron concentration, and temperature

The spin crossover of iron in Fe$^{3+}$-bearing bridgmanite, the most abundant mineral of the Earth's lower mantle, is by now a well-established phenomenon, though several aspects of this crossover remain unclear. Here we investigate effects of disorder, iron concentration, and temperature on this crossover using ab initio LDA + U$_{sc}$ calculations. The effect of concentration and disorder are addressed using complete statistical samplings of coupled substituted configurations in super-cells containing up to 80 atoms. Vibrational/thermal effects on the crossover are addressed within the quasiharmonic approximation. The effect of disorder seems quite small, while increasing iron concentration results in considerable increase in crossover pressure. Our calculated compression curves for iron-free, Fe$^{2+}$-, and Fe$^{3+}$-bearing bridgmanite compare well with the latest experimental measurements. The comparison also suggests that in a close system, Fe$^{2+}$ present in the sample may transform into Fe$^{3+}$ by introduction of Mg and O vacancies with increasing pressure. As in the spin crossover in ferropericlase, this crossover in bridgmanite is accompanied by a clear volume reduction and an anomalous softening of the bulk modulus throughout the crossover pressure range. These effects reduce significantly with increasing temperature. Though the concentration of [Fe$^{3+}$]$_{Si}$ in bridgmanite may be small, related elastic anomalies may impact the interpretation of radial and lateral velocity structures of the Earth's lower mantle.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Thermoelasticity of Fe$^{3+}$- and Al-bearing bridgmanite

We report \textit{ab initio} (LDA + U$_{sc}$) calculations of thermoelastic properties of ferric iron (Fe$^{3+}$)- and aluminum (Al)-bearing bridgmanite (MgSiO$_3$ perovskite), the main Earth forming phase, at relevant pressure and temperature conditions and compositions. Three coupled substitutions, namely, [Al]$_{Mg}$-[Al]$_{Si}$, [Fe$^{3+}$]$_{Mg}$-[Fe$^{3+}$]$_{Si}$, and [Fe$^{3+}$]$_{Mg}$-[Al]$_{Si}$ have been investigated. Aggregate elastic moduli and sound velocities are successfully compared with limited experimental data available. In the case of the [Fe$^{3+}$]$_{Mg}$-[Fe$^{3+}$]$_{Si}$ substitution, the high-spin (S=5/2) to low-spin (S=1/2) crossover in [Fe$^{3+}$]$_{Si}$ induces a volume collapse and elastic anomalies across the transition region. However, the associated anomalies should disappear in the presence of aluminum in the most favorable substitution, i.e., [Fe$^{3+}$]$_{Mg}$-[Al]$_{Si}$. Calculated elastic properties along a lower mantle model geotherm suggest that the elastic behavior of bridgmanite with simultaneous substitution of Fe$_{2}$O$_3$ and Al$_{2}$O$_3$ in equal proportions or with Al$_{2}$O$_3$ in excess should be similar to that of (Mg,Fe$^{2+}$)SiO$_3$ bridgmanite. Excess Fe$_{2}$O$_3$ should produce elastic anomalies though.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Iron spin crossover and its influence on post-perovskite transitions in MgSiO$_3$ and MgGeO$_3$

MgGeO$_3$-perovskite is known to be a low-pressure analog of MgSiO$_3$-perovskite in many respects, but especially in regard to the post-perovskite transition. As such, investigation of spin state changes in Fe-bearing MgGeO$_3$ might help to clarify some aspects of this type of state change in Fe-bearing MgSiO$_3$. Using DFT+U calculations, we have investigated pressure induced spin state changes in Fe$^{2+}$ and Fe$^{3+}$ in MgGeO$_3$ perovskite and post-perovskite. Owing to the relatively larger atomic size of germanium compared to silicon, germanate phases have larger unit cell volume and inter-atomic distances than equivalent silicate phases at same pressures. As a result, all pressure induced state changes in iron occur at higher pressures in germanate phases than in the silicate ones, be it a spin state change or position change of (ferrous) iron in the perovskite cage. We showed that iron state transitions occur at particular average Fe-O bond-length irrespective of mineral composition (silicate or germanate) or functionals (LDA+U$_{sc}$ or GGA+U$_{sc}$). Ferrous iron substitution decreases the perovskite to post-perovskite (PPv) transition pressure while coupled ferric iron substitution increases it noticeably.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Post-stishovite transition in hydrous aluminous SiO2

Lakshtanov et al. (2007) showed that incorporation of aluminum and some water into SiO2 significantly reduces the post-stishovite transition pressure in SiO2. This discovery suggested that the ferroelastic post-stishovite transition in subducted MORB crust could be the source of reflectors/scatterers with low shear velocities observed in the mid to upper lower mantle. A few years later, a similar effect was observed in anhydrous Al-bearing silica. In this paper, we show using first principles static calculations and molecular dynamics using inter-atomic potentials that hydrogen bonds and hydrogen mobility play a crucial role in lowering the post-stishovite transition pressure. A cooperative redistribution of hydrogen atoms is the main mechanism responsible for the transition pressure reduction in hydrous aluminous stishovite. The effect is enhanced by increasing hydrogen concentration. This perspective suggests a potential relationship between the depth of seismic scatterers and the water content in stishovite.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Thermoelasticity of Fe2+-bearing bridgmanite

We present LDA+U calculations of high temperature elastic properties of bridgmanite with composition (Mg$_{(1-x)}$Fe$_{x}^{2+}$)SiO$_3$ for $0\le{x}\le0.125$. Results of elastic moduli and acoustic velocities for the Mg-end member (x=0) agree very well with the latest high pressure and high temperature experimental measurements. In the iron-bearing system, we focus particularly on the change in thermoelastic parameters across the state change that occurs in ferrous iron above $\sim$30 GPa, often attributed to a high-spin (HS) to intermediate spin (IS) crossover but explained by first principles calculations as a lateral displacement of substitutional iron in the perovskite cage. We show that the measured effect of this change on the equation of state of this system can be explained by the lateral displacement of substitutional iron, not by the HS to IS crossover. The calculated elastic properties of (Mg$_{0.875}$Fe$_{0.125}^{2+}$)SiO$_3$ along an adiabatic mantle geotherm, somewhat overestimate longitudinal velocities but produce densities and shear velocities quite consistent with Preliminary Reference Earth Model data throughout most of the lower mantle.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

First-principles study of intermediate-spin ferrous iron in the Earth's lower mantle

Spin crossover of iron is of central importance in solid Earth geophysics. It impacts all physical properties of minerals that altogether constitute $\sim 95$ vol\% of the Earth's lower mantle: ferropericlase [(Mg,Fe)O] and Fe-bearing magnesium silicate (MgSiO$_3$) perovskite. Despite great strides made in the past decade, the existence of intermediate-spin (IS) state in ferrous iron (Fe$^{2+}$) (with total electron spin $S=1$) and its possible role in the pressure-induced spin crossover in these lower-mantle minerals still remain controversial. Using density functional theory $+$ self-consistent Hubbard $U$ (DFT$+U_{sc}$) calculations, we investigate all possible types of IS states of Fe$^{2+}$ in (Mg,Fe)O and (Mg,Fe)SiO$_3$ perovskite. Among the possible IS states in these minerals, the most probable IS state has an electronic configuration that significantly reduces the electron overlap and the iron nuclear quadrupole splitting (QS). These most probable IS states, however, are still energetically disfavored, and their QSs are inconsistent with Mössbauer spectra. We therefore conclude that IS Fe$^{2+}$ is highly unlikely in the Earth's lower mantle.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Searching for high magnetization density in bulk Fe: the new metastable Fe$_6$ phase

We report the discovery of a new allotrope of iron by first principles calculations. This phase has $Pmn2_1$ symmetry, a six-atom unit cell (hence the name Fe$_6$), and the highest magnetization density (M$_s$) among all known crystalline phases of iron. Obtained from the structural optimizations of the Fe$_3$C-cementite crystal upon carbon removal, $Pmn2_1$ Fe$_6$ is shown to result from the stabilization of a ferromagnetic FCC phase, further strained along the Bain path. Although metastable from 0 to 50 GPa, the new phase is more stable, at low pressures, than the other well-known HCP and FCC allotropes and smoothly transforms into the FCC phase under compression. If stabilized to room temperature, e.g., by interstitial impurities, Fe$_{6}$ could become the basis material for high M$_s$ rare-earth-free permanent magnets and high-impact applications such as, light-weight electric engine rotors or high-density recording media. The new phase could also be key to explain the enigmatic high M$_s$ of Fe$_{16}$N$_2$, which is currently attracting an intense research activity.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Phonon Quasi-Particles and Anharmonic Free Energy in Complex Systems

We use a hybrid strategy to obtain anharmonic frequency shifts and lifetimes of phonon quasi-particles from first principles molecular dynamics simulations in modest size supercells. This approach is effective irrespective of crystal structure complexity and facilitates calculation of full anharmonic phonon dispersions, as long as phonon quasi-particles are well defined. We validate this approach to obtain anharmonic effects with calculations in MgSiO3-perovskite, the major Earth forming mineral phase. First, we reproduce irregular temperature induced frequency shifts of well characterized Raman modes. Second, we combine the phonon gas model (PGM) with quasi-particle frequencies and reproduce free energies obtained using thermodynamic integration. Combining thoroughly sampled quasi-particle dispersions with the PGM we then obtain first-principles anharmonic free energy in the thermodynamic limit (N right arrow infinity).

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Ferromagnetic insulating state in tensile-strained LaCoO$_3$ thin films

With local density approximation + Hubbard $U$ (LDA+$U$) calculations, we show that the ferromagnetic (FM) insulating state observed in tensile-strained LaCoO$_3$ epitaxial thin films is most likely a mixture of low-spin (LS) and high-spin (HS) Co, namely, a HS/LS mixture state. Compared with other FM states, including the intermediate-spin (IS) state (\textit{metallic} within LDA+$U$), which consists of IS Co only, and the insulating IS/LS mixture state, the HS/LS state is the most favorable one. The FM order in HS/LS state is stabilized via the superexchange interactions between adjacent LS and HS Co. We also show that Co spin state can be identified by measuring the electric field gradient (EFG) at Co nucleus via nuclear magnetic resonance (NMR) spectroscopy.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

First principles study of electronic and structural properties of CuO

We investigate the electronic and structural properties of CuO, which shows significant deviations from the trends obeyed by other transition-metal monoxides. Using an extended Hubbard corrective functional, we uncover an orbitally ordered insulating ground state for the cubic phase of this material, which was expected but never found before. This insulating state results from a fine balance between the tendency of Cu to complete its d-shell and Hund's rule magnetism. Starting from the ground state for the cubic phase, we also study tetragonal distortions of the unit cell (recently reported in experiments), the consequent electronic reorganizations and identify the equilibrium structure. Our calculations reveal an unexpected richness of possible magnetic and orbital orders, relatively close in energy to the ground state, whose stability depends on the sign and entity of distortion.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

New ultrahigh pressure phases of H2O ice predicted using an adaptive genetic algorithm

We propose three new phases of H2O under ultrahigh pressure. Our structural search was performed using an adaptive genetic algorithm which allows an extensive exploration of crystal structure. The new sequence of pressure-induced transitions beyond ice X at 0 K should be ice X - Pbcm - Pbca - Pmc21 - P21 - P21/c phases. Across the Pmc21 - P21 transition, the coordination number of oxygen increases from 4 to 5 with a significant increase of density. All stable crystalline phases have nonmetallic band structures up to 7 TPa.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Spin-state crossover and hyperfine interactions of ferric iron in MgSiO$_3$ perovskite

Using density functional theory plus Hubbard $U$ calculations, we show that the ground state of (Mg,Fe)(Si,Fe)O$_3$ perovskite, a major mineral phase in the Earth's lower mantle, has high-spin ferric iron ($S=5/2$) at both the dodecahedral (A) and octahedral (B) site. As the pressure increases, the B-site iron undergoes a spin-state crossover to the low-spin state ($S=1/2$), while the A-site iron remains in the high-spin state. Our calculation shows that the B-site spin-state crossover in the pressure range of 40-70 GPa is accompanied by a noticeable volume reduction and an increase in quadrupole splitting, consistent with recent X-ray diffraction and Mössbauer spectroscopy measurements. The volume reduction leads to a significant softening in the bulk modulus, which suggests a possible source of seismic velocity anomalies in the lower mantle.

0 citations0 reviewsNo code linkedOpen access