Source author record

Claudio Cazorla

Claudio Cazorla appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

32works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

32 published item(s)

preprint2023arXiv

Absence of off-diagonal long-range order in hcp $^{\bf 4}$He dislocation cores

The mass transport properties along dislocation cores in hcp $^4$He are revisited by considering two types of edge dislocations as well as a screw dislocation, using a fully correlated quantum simulation approach. Specifically, we employ the zero-temperature path-integral ground state (PIGS) method together with ergodic sampling of the permutation space to investigate the fundamental dislocation core structures and their off-diagonal long-range order properties. It is found that the Bose-Einstein condensate fraction of such defective $^4$He systems is practically null ($\le 10^{-6}$), just as in the bulk defect-free crystal. These results provide compelling evidence for the absence of intrinsic superfluidity in dislocation cores in hcp $^4$He and challenge the superfluid dislocation-network interpretation of the mass-flux-experiment observations, calling for further experimental investigation.

preprint2022arXiv

Atomistic simulation of barocaloric effects

Due to critical environmental issues there is a pressing need to switch from current refrigeration methods based on compression of greenhouse gases to novel solid-state cooling technologies. Solid-state cooling capitalizes on the thermal response of materials to external fields named "caloric effect". The barocaloric (BC) effect driven by hydrostatic pressure is particularly promising from a technological point of view since typically presents larger cooling potential than other caloric variants (e.g., magnetocaloric and electrocaloric effects driven by magnetic and electric fields, respectively). Atomistic simulation of BC effects represents an efficient and physically insightful strategy for advancing solid-state cooling by complementing, and in some cases even guiding, experiments. Atomistic simulation of BC effects involves approaches ranging from computationally inexpensive force fields to computationally very demanding, but quantitatively accurate, first-principles methods. Here, we survey several methods and strategies involved in atomistic simulation of BC effects like the quasi-harmonic approximation and direct/quasi-direct estimation approaches. The Review finalizes with a collection of case studies in which some of these methods were employed to simulate and predict original BC effects.

preprint2022arXiv

Direct Hydrogen Production from Water/Seawater by Irradiation/Vibration-Activated Using Defective Ferroelectric BaTiO3-x Nanoparticles

Hydrogen is a promising fossil-fuel alternative fuel owing to its environmentally neutral emissions and high energy density. However, the need for purified water and external power are critical hindrances to implementation of hydrogen production. The present work reveals the potential to overcome these shortcomings through piezo-photocatalysis of seawater using BaTiO3-x (BTO) nanoparticles. This material was made piezoelectrically active by annealing under different atmospheres, including O2, N2, Ar, and H2, the latter of which caused Ti4+ to Ti(4-x)+ multiple reductions and structural expansions that stabilized piezoelectric tetragonal BTO domains. The resultant defect equilibria combine ionic and electron effects, including Ti redox reactions, charge-compensating surface oxygen vacancy formation, and color centre alterations. Further, variety of experimental techniques revealed the effects of reduction on the energy band structure. A strong piezoelectric effect and the presence of self-polarization were confirmed by piezoresponse force microscopy, while simulation work clarified the role of vibration on band bending deriving from the former. The performance data contrasted H2 evolution using deionized (DI) water, simulated seawater, and natural seawater subjected to photocatalysis, piezocatalysis, and piezo-photocatalysis. An efficient H2 evolution rate of 132.4 micromol/g/h was achieved from DI water using piezo-photocatalysis for 5 h. In contrast, piezocatalysis for 2 h followed by piezo-photocatalysis for 3 h resulted in H2 evolution rates of 100.7 micromol/g/h for DI water, 63.4 micromol/g/h for simulated seawater, and 48.7 micromol/g/h for natural seawater. This work provides potential new strategies for large-scale green H2 production using abundant natural resources with conventional piezoelectric material while leveraging the effects of ions dissolved in seawater.

preprint2022arXiv

High-throughput screening of piezo-photocatalytic materials for hydrogen production

Finding cost-effective and efficient photocatalytic materials able to catalyse the water splitting reaction under visible light is one of the greatest challenges in current environmental material science. Despite that many photocatalysts are already known in the context of green hydrogen production, strategies to systematically and rationally modify their optoelectronic properties to achieve desired photocatalytic performance are yet to be established. Piezoelectric materials react to mechanical stimuli by adjusting their band gaps and band alignments, thus offering a possible route to precise photocatalyst design. However, piezo-photocatalysts are relatively scarce and have been seldom investigated to date. Here, we present a high-throughput screening of piezo-photocatalytic materials performed over $\sim 1,000$ bulk piezoelectrics that relies on a simple electrostatic model and first-principles calculations. A total of $\sim 10$ previously overlooked binary and tertiary bulk compounds are theoretically identified as highly promising piezo-photocatalysts due to their appropriate optoelectronic properties and superb band alignment tunability driven by uniaxial strain.

preprint2020arXiv

Colossal barocaloric effects in the complex hydride Li$_{2}$B$_{12}$H$_{12}$

Traditional refrigeration technologies based on compression cycles of greenhouse gases pose serious threats to the environment and cannot be downscaled to electronic device dimensions. Solid-state cooling exploits the thermal response of caloric materials to external fields and represents a promising alternative to current refrigeration methods. However, most of the caloric materials known to date present relatively small adiabatic temperature changes ($|ΔT| \sim 1$ K) and/or limiting irreversibility issues resulting from significant phase-transition hysteresis. Here, we predict the existence of colossal barocaloric effects (isothermal entropy changes of $|ΔS| \sim 100$ JK$^{-1}$kg$^{-1}$) in the energy material Li$_{2}$B$_{12}$H$_{12}$ by means of molecular dynamics simulations. Specifically, we estimate $|ΔS| = 387$ JK$^{-1}$kg$^{-1}$ and $|ΔT| = 26$ K for an applied pressure of $P = 0.4$ GPa at $T = 475$ K. The disclosed colossal barocaloric effects are originated by an order-disorder phase transformation that exhibits a fair degree of reversibility and involves coexisting Li$^{+}$ diffusion and (BH)$_{12}^{-2}$ reorientational motion at high temperatures.

preprint2020arXiv

Combining density functional theories to correctly describe the energy, lattice structure and electronic density of functional oxide perovskites

Functional oxide perovskites are the pillar of cutting-edge technological applications. Density functional theory (DFT) simulations are the theoretical methods of choice to understand and design perovskite materials. However, tests on the reliability of DFT to describe fundamental properties of oxide perovskites are scarce and mostly ill-defined due to a lack of rigorous theoretical benchmarks for solids. Here, we present a quantum Monte Carlo benchmark study of DFT on the archetypal perovskite BaTiO$_{3}$ (BTO). It shows that no DFT approximation can simultaneously reproduce the energy, structure, and electronic density of BTO. Traditional protocols to select DFT approximations are empirical and fail to detect this shortcoming. An approach combining two different non-empirical DFT schemes, "SCAN" and "HSE06", is able to holistically describe BTO with accuracy. Combined DFT approaches should thus be considered as a promising alternative to standard methods for simulating oxide perovskites.

preprint2020arXiv

Flexible complementary logic circuit built from two identical organic electrochemical transistors

The organic electrochemical transistor (OECT) with a conjugated polymer as the active material is the elementary unit of organic bioelectronic devices. Increased functionalities, such as low power consumption, can be achieved by building complementary circuits featuring two or more OECTs. Complementary circuits commonly combine both p- and n-type transistors to reduce power draw. While p-type OECTs are readily available, n-type OECTs are less common mainly due to poor stability of the n-type active channel material in aqueous electrolyte. Additionally, an OECT based complementary circuit requires well matched transport properties in the p- and n-type materials. Here, a complementary circuit is made using a pair of OECTs having polyaniline (PANI) as the channel material in both transistors. PANI is chosen due to its unique behaviour exhibiting a peak in current versus gate voltage when used as an active channel in an OECT. The PANI based circuit is shown to have excellent performance with gain of ~ 7 and could be transferred on a flexible biocompatible chitosan substrate with demonstrated operation in aqueous electrolyte. This study extends the capabilities of conjugated polymer based OECTs.

preprint2020arXiv

Giant Thermal Enhancement of the Electric Polarization in Ferrimagnetic BiFe$_{1-x}$Co$_{x}$O$_{3}$ Solid Solutions Near Room Temperature

Thermal excitations typically reduce the electric polarization in ferroelectric materials. Here, we show by means of first-principles calculations that multiferroic BiFe$_{1-x}$Co$_{x}$O$_{3}$ solid solutions with $0.25 \le x \le 0.50$ (BFCO) represent a noteworthy exception to this behaviour. In particular, we find that at room temperature and for moderate pressures of $0.1$-$1.0$ GPa, depending on the composition, the electric polarization of bulk BFCO increases by $\sim 200$%. The origin of such an exceptional behavior is a phase transformation involving a low-$T$ rhombohedral (${\cal R}$) phase and a high-$T$ super-tetragonal (${\cal T}$) phase. Both ${\cal R}$ and ${\cal T}$ phases are ferrimagnetic near room temperature with an approximate net magnetization of $0.13$$μ_{B}$ per formula unit. Contrarily to what occurs in either bulk BiFeO$_{3}$ or BiCoO$_{3}$, the ${\cal T}$ phase is stabilized over the ${\cal R}$ by increasing temperature due to its higher vibrational entropy. This extraordinary $T$-induced ${\cal R} \to {\cal T}$ phase transition is originated by polar phonon modes involving concerted displacements of transition-metal and oxygen ions.

preprint2020arXiv

Is the normal to superionic transformation occurring in type-II fast-ion conductors a real thermodynamic phase transition?

Fedorov has written a Comment on our recent paper J. Phys. Chem. C 112, 1267 (2018) in which he suggests that the transformation from a normal ($α$) to a superionic ($β$) state occurring in type-II fast-ion conductors (e.g., CaF$_{2}$ and Li-based superionic materials) cannot be regarded as a phase transition in a rigorous thermodynamic sense. Fedorov's arguments for such a criticism can be summarized as follows: (1) the available heat capacity data for the $α\to β$ transformation do not reproduce the expected behaviour for a thermodynamic phase transition, (2) there is not evidence for molar volume discontinuity during the superionic $α\to β$ transformation, and (3) the $α\to β$ transformation is symmetry preserving. Here, we demonstrate either by reproducing experimental data published by other authors or by explaining some simple solid-state arguments that the statements put forward by Fedorov are incorrect and do not correspond to reality. In fact, based on mounting experimental evidence, the normal to superionic phase transformation occurring in type-II fast-ion conductors can be rigorously considered, and should be referred to, as a phase transition. Notwithstanding, it remains less obvious to establish with generality which is the order of the $α\to β$ phase transition.

preprint2020arXiv

Simulation of the electrocaloric effect based on first-principles methods

Due to critical environmental and technological issues, there is a pressing need to switch from current refrigeration methods based on compression of gases to novel solid-state cooling technologies. Solid-state cooling is based on the thermal response of materials to external magnetic, electric, or mechanic fields, the so-called caloric effect. The electrocaloric (EC) effect, which is caused by electric fields and typically occurs in polar materials, is particularly promising from a technological point of view owing to its good scalability and natural implementation in circuitry. Simulation of EC effects represents an efficient and physically insightful strategy for advancing the field of solid-state cooling by complementing, and in some cases guiding, experiments. Theoretical estimation of EC effects can be achieved with different approaches ranging from computationally inexpensive but physically insightful phenomenological free-energy models to computationally very demanding and quantitatively accurate first-principles methods. In this Chapter, we review EC simulation approaches that rely on first-principles methods. In this category, we include ab initio quasi-harmonic methods, bond-valence and classical interatomic potentials and effective Hamiltonians. In analogy to the experiments, these simulation approaches can be used to estimate EC effects either directly or indirectly and we review here well-established protocols that can be followed for each case. The Chapter finalises with a collection of representative examples in which first-principles based approaches have been used to predict and understand original EC effects.

preprint2019arXiv

First-principles prediction of extraordinary thermoelectric efficiency in superionic Li2SnX3(X=S,Se)

Thermoelectric materials create an electric potential when subject to a temperature gradient and vice versa hence they can be used to harvest waste heat into electricity and in thermal management applications. However, finding highly efficient thermoelectrics with high figures of merit, zT$\geq$1, is very challenging because the combination of high power factor and low thermal conductivity is rare in materials. Here, we use first-principles methods to analyze the thermoelectric properties of Li$_2$Sn$X_3$ ($X$=S,Se), a recently synthesized class of lithium fast-ion conductors presenting high thermal stability. In p-type Li$_2$Sn$X_3$, we estimate highly flat electronic valence bands that render high Seebeck coefficients exceeding 400 $μ$VK$^{-1}$ at 700K. In n-type Li$_2$Sn$X_3$, the electronic conduction bands are slightly dispersive however the accompanying weak electron-acoustic phonon scattering induces high electrical conductivity. The combination of high Seebeck coefficient and electrical conductivity gives rise to high power factors, reaching a maximum of 4 mWm$^{-1}$K$^{-2}$ in p-type Li$_2$SnS$_3$ and 8 mWm$^{-1}$K$^{-2}$ in n-type Li$_2$SnSe$_3$ at 300 K. Likewise, the thermal conductivity in Li$_2$Sn$X_3$ is low as compared to conventional thermoelectric materials, 2-5 Wm$^{-1}$K$^{-1}$ at room temperature. As a result, we estimate a maximum zT = 1.05 in p-type Li$_2$SnS$_3$ at 700 K and an extraordinary 3.07 (1.5) in n-type Li$_2$SnSe$_3$ at the same temperature (300 K). Our findings of huge zT in Li$_2$Sn$X_3$ suggest that lithium fast-ion conductors, typically employed as electrolytes in solid-state batteries, hold exceptional promise as thermoelectric materials.

preprint2019arXiv

Prediction of large barocaloric effects in thermoelectric superionic materials

We predict the existence of large barocaloric effects above room temperature in the thermoelectric fast-ion conductor Cu$_{2}$Se by using classical molecular dynamics simulations and first-principles computational methods. A hydrostatic pressure of $1$ GPa induces large isothermal entropy changes of $|ΔS| \sim 15$-$45$ Jkg$^{-1}$K$^{-1}$ and adiabatic temperature shifts of $|ΔT| \sim 10$ K in the temperature interval $400 \le T \le 700$ K. Structural phase transitions are absent in the analysed thermodynamic range. The causes of such large barocaloric effects are significant $P$-induced variations on the ionic conductivity of Cu$_{2}$Se and the inherently high anharmonicity of the material. Uniaxial stresses of the same magnitude, either compressive or tensile, produce comparatively much smaller caloric effects, namely, $|ΔS| \sim 1$ Jkg$^{-1}$K$^{-1}$ and $|ΔT| \sim 0.1$ K, due to practically null influence on the ionic diffusivity of Cu$_{2}$Se. Our simulation work shows that thermoelectric compounds presenting high ionic disorder, like copper and silver-based chalcogenides, may render large mechanocaloric effects and thus are promising materials for engineering solid-state cooling applications that do not require the application of electric fields.

preprint2016arXiv

Multiple structural transitions driven by spin-phonon couplings in a perovskite oxide

Spin-phonon interactions are central to many interesting phenomena, ranging from superconductivity to magnetoelectric effects. Yet, they are believed to have a negligible influence on the structural behavior of most materials. For example, magnetic perovskite oxides often undergo structural transitions accompanied by magnetic signatures whose minuteness suggests that the underlying spin-phonon couplings are largely irrelevant. Here we present an exception to this rule, showing that novel effects can occur as a consequence. Our first-principles calculations reveal that spin-phonon interactions are essential to reproduce the experimental observations on the phase diagram of magnetoelectric multiferroic BiCoO$_{3}$. Moreover, we predict that, under compression, these couplings lead to an unprecedented temperature-driven double-reentrant sequence of ferroelectric transitions. We propose how to modify BiCoO$_{3}$ via chemical doping to reproduce such striking effects at ambient conditions, thereby yielding useful multifunctionality.

preprint2016arXiv

Prediction of giant mechanocaloric effects in fluorite-structured superionic materials

Mechanocaloric materials experience a change in temperature when a mechanical stress is adiabatically applied on them. Thus far, only ferroelectrics and superelastic metallic alloys have been considered as potential mechanocaloric compounds to be exploited in solid-state cooling applications. Here we show that giant mechanocaloric effects occur in hitherto overlooked fast ion conductors (FIC), a class of multicomponent materials in which above a critical temperature, Ts, a constituent ionic species undergoes a sudden increase in mobility. Using first-principles and molecular dynamics simulations, we found that the superionic transition in fluorite-structured FIC, which is characterised by a large entropy increase of the order of 100 J/K*Kg, can be externally tuned with hydrostatic, biaxial or uniaxial stresses. In particular, Ts can be reduced several hundreds of degrees through the application of moderate tensile stresses due to the concomitant drop in the formation energy of Frenkel pair defects. We predict that the adiabatic temperature change in CaF2 and PbF2, two archetypal fluorite-structured FIC, close to their critical points are of the order of 100 and 10 K, respectively. This work advocates that FIC constitute a new family of mechanocaloric materials showing great promise for prospective solid-state refrigeration applications.

preprint2016arXiv

Revisiting the zero-temperature phase diagram of stoichiometric SrCoO3 with first-principles methods

By using first-principles methods based on density functional theory we revisited the zero-temperature phase diagram of stoichiometric SrCoO3, a ferromagnetic metallic perovskite that undergoes significant structural, electronic, and magnetic changes as its content of oxygen is decreased. We considered both bulk and epitaxial thin film geometries. In the bulk case, we found that a tetragonal P4/mbm phase with moderate Jahn-Teller distortions and c/a ratio of ~1/sqrt{2} is consistently predicted to have a lower energy than the thus far assumed ground-state cubic Pm-3m phase. In thin films, we found two phase transitions occurring at compressive and tensile epitaxial strains. However, in contrast to previous theoretical predictions, our results show that: (i) the phase transition induced by tensile strain is isostructural and involves only a change in magnetic spin order (that is, not a metallic to insulator transformation), and (ii) the phase transition induced by compressive strain comprises simultaneous structural, electronic and magnetic spin order changes, but the required epitaxial stress is so large (<-6%) that is unlikely to be observed in practice. Our findings call for a revision of the crystallographic analysis performed in fully oxidised SrCoO3 samples at low temperatures, as well as of previous first-principles studies.

preprint2015arXiv

Electrostatic engineering of strained ferroelectric perovskites from first-principles

Design of novel artificial materials based on ferroelectric perovskites relies on the basic principles of electrostatic coupling and in-plane lattice matching. These rules state that the out-of-plane component of the electric displacement field and the in-plane components of the strain are preserved across a layered superlattice, provided that certain growth conditions are respected. Intense research is currently directed at optimizing materials functionalities based on these guidelines, often with remarkable success. Such principles, however, are of limited practical use unless one disposes of reliable data on how a given material behaves under arbitrary electrical and mechanical boundary conditions. Here we demonstrate, by focusing on the prototypical ferroelectrics PbTiO3 and BiFeO3 as testcases, how such information can be calculated from first principles in a systematic and efficient way. In particular, we construct a series of two-dimensional maps that describe the behavior of either compound (e.g. concerning the ferroelectric polarization and antiferrodistortive instabilities) at any conceivable choice of the in-plane lattice parameter, a, and out-of-plane electric displacement, D. In addition to being of immediate practical applicability to superlattice design, our results bring new insight into the complex interplay of competing degrees of freedom in perovskite materials, and reveal some notable instances where the behavior of these materials depart from what naively is expected.

preprint2015arXiv

First-principles modeling of three-body interactions in highly compressed solid helium

We present a new set of three-body interaction models based on the Bruch-McGee (BM) potential that are suitable for the study of the energy, structural and elastic properties of solid 4He at high pressure. Our ab initio three-body potentials are obtained from the fit to total energies and atomic forces computed with the van der Waals density functional theory method due to Grimme, and represent an improvement with respect to previously reported three-body interaction models. In particular, we show that some of the introduced BM parametrizations reproduce closely the experimental equation of state and bulk modulus of solid helium up to a pressure of ~ 60 GPa, when used in combination with standard pairwise interaction models in diffusion Monte Carlo simulations. Importantly, we find that recent predictions reporting a surprisingly small variation of the kinetic energy and Lindeman ratio on quantum crystals under increasing pressure are likely to be artifacts produced by the use of incomplete interaction models. Also, we show that the experimental variation of the shear modulus, C44, at P < 25 GPa can be quantitatively described with the new set of three-body BM potentials. At higher pressures, however, the agreement between our C44 results and experiments deteriorates and thus we argue that higher order many-body terms in the expansion of the atomic interactions probably are necessary in order to better describe elasticity in very dense solid 4He.

preprint2015arXiv

In the search of new electrocaloric materials: Fast ion conductors

We analyse the effects of applying an electric field on the critical temperature, Ts, at which superionicity appears in archetypal fast ion conductor CaF2 by means of molecular dynamics simulations. We find that the onset of superionicity can be reduced by about 100K when relatively small electric fields of ~50KV/cm are employed. Under large enough electric fields, however, ionic conductivity is depleted. The normal to superionic phase transition is characterised by a large increase of entropy, thereby sizeable electrocaloric effects can be realised in fast ion conductors that are promising for solid-state cooling applications.

preprint2015arXiv

The Role of Density Functional Theory Methods in the Prediction of Nanostructured Gas-Adsorbent Materials

With the advent of new synthesis and large-scale production technologies, nanostructured gas-adsorbent materials (GAM) like carbon nanocomposites and metal-organic frameworks are becoming increasingly more influential in our everyday lives. First-principles methods based on density functional theory (DFT) have been pivotal in establishing the rational design of GAM, a factor which has tremendously boosted their development. However, DFT methods are not perfect and due to the stringent accuracy thresholds demanded in modelling of GAM (i.e., exact binding energies to within ~0.01 eV) these techniques may provide erroneous conclusions in some challenging situations. Examples of problematic circumstances include gas-adsorption processes in which both electronic long-range exchange and nonlocal correlations are important, and systems where many-body energy and Coulomb screening effects cannot be disregarded. In this critical review, we analyse recent efforts done in the assessment of the performance of DFT methods in the prediction and understanding of GAM. Our inquiry is constrained to the areas of hydrogen storage and carbon capture and sequestration, for which we expose a number of unresolved modelling controversies and define a set of best practice simulation principles. Also, we identify the subtle problems found in the generalization of DFT benchmark studies performed in model cluster systems to real materials, and discuss effective approaches to circumvent them. The increasing awareness of the strengths and imperfections of DFT methods in the simulation of gas-adsorption phenomena should lead in the medium term to more precise, and hence even more fruitful, ab initio engineering of GAM.

preprint2014arXiv

Ab initio design of charge-mismatched ferroelectric superlattices

We present a systematic approach to modeling the electrical and structural properties of charge-mismatched superlattices from first principles. Our strategy is based on bulk calculations of the parent compounds, which we perform as a function of in-plane strain and out-of-plane electric displacement field. The resulting two-dimensional phase diagrams allow us to accurately predict, without performing further calculations, the behavior of a layered heterostructure where the aforementioned building blocks are electrostatically and elastically coupled, with an arbitrary choice of the interface charge (originated from the polar discontinuity) and volume ratio. By using the [PbTiO3]_m/[BiFeO3]_n system as test case, we demonstrate that interface polarity has a dramatic impact on the ferroelectric behavior of the superlattice, leading to the stabilization of otherwise inaccessible bulk phases.

preprint2014arXiv

First-Principles Modeling of Quantum Nuclear Effects and Atomic Interactions in Solid 4He at High Pressure

We present a first-principles computational study of solid 4He at T=0K and pressures up to 160GPa. Our computational strategy consists in using van der Waals density functional theory (DFT-vdW) to describe the electronic degrees of freedom in this material, and the diffusion Monte Carlo (DMC) method to solve the Schrodinger equation describing the behavior of the quantum nuclei. For this, we construct an analytical interaction function based on the pairwise Aziz potential that closely matches the volume variation of the cohesive energy calculated with DFT-vdW in dense helium. Interestingly, we find that the kinetic energy of solid 4He does not increase appreciably with compression for P > 85GPa. Also, we show that the Lindemann ratio in dense solid 4He amounts to 0.10 almost independently of pressure. The reliability of customary quasi-harmonic DFT (QH DFT) approaches in the description of quantum nuclear effects in solids is also studied. We find that QH DFT simulations, although provide a reasonable equation of state in agreement with experiments, are not able to reproduce correctly these critical effects in compressed 4He. In particular, we disclose huge discrepancies of at least ~50% in the calculated 4He kinetic energies using both the QH DFT and present DFT-DMC methods.

preprint2014arXiv

Superionicity and Polymorphism in Calcium Fluoride at High Pressure

We present a combined experimental and computational first-principles study of the superionic and structural properties of CaF2 at high P-T conditions. We observe an anomalous superionic behavior in the low-P fluorite phase that consists in a decrease of the normal-> superionic critical temperature with compression. This unexpected effect can be explained in terms of a P-induced softening of a zone-boundary $X$ phonon which involves exclusively fluorine displacements. Also we find that superionic conductivity is absent in the high-P cotunnite phase. Instead, superionicity develops in a new low-symmetry high-T phase that we identify as monoclinic (space group P2_1/c). We discuss the possibility of observing these intriguing phenomena in related isomorphic materials.

preprint2014arXiv

Supporting Material for: Superionicity and Polymorphism in Calcium Fluoride at High Pressure

In this supplemental document we explain the details of our ab initio density functional theory calculations and crystal structure searches performed in CaF2 under pressure. Also we report the computed P-dependence of the formation energy of Frenkel pair defects in the cubic fluorite phase, and the vibrational phonon spectra and structural data of the predicted high-T monoclinic P2_1/c phase. Calculated enthalpy energies in compressed SrF2 and BaF2 are also presented.

preprint2014arXiv

The Limit of Mechanical Stability in Quantum Crystals: A Diffusion Monte Carlo Study of Solid 4He

We present a first-principles study of the energy and elastic properties of solid helium at pressures below the range in which is energetically stable. We find that the limit of mechanical stability in hcp 4He is $P_{s}$ = -33.82 bar, which lies significantly below the spinodal pressure found in the liquid phase (i.e., -9.6 bar). Furthermore, we show that the pressure variation of the transverse and longitudinal sound velocities close to $P_{s}$ do not follow a power law of the form $\propto \left( P - P_{s} \right)^γ$, in contrast to what is observed on the fluid.

preprint2013arXiv

Accuracy of Density Functional Theory in Prediction of Carbon Dioxide Adsorbent Materials

We have performed a thorough computational study to assess the accuracy of density functional theory (DFT) methods in describing the interactions of CO2 with model alkali-earth-metal (AEM, Ca and Li) decorated carbon structures, namely anthracene (C14H10) molecules. We find that gas-adsorption energy and equilibrium structure results obtained with both standard (i.e. LDA and GGA) and hybrid (i.e. PBE0 and B3LYP) exchange-correlation functionals of DFT differ significantly from results obtained with second-order Moller-Plesset perturbation theory (MP2), an accurate computational quantum chemistry method. The major disagreements found can be mostly rationalized in terms of electron correlation errors that lead to inaccurate charge transfers and electrostatic Coulomb interactions between the molecules. Interestingly, we show that when the concentration of AEM atoms in anthracene is tuned to resemble as closely as possible to the electronic structure of AEM-decorated graphene, hybrid exchange-correlation DFT and MP2 methods provide quantitatively similar results. We discuss the implications of our work in modeling and characterization of currently sought carbon capture and sequestration materials.

preprint2013arXiv

Elastic Constants of Incommensurate Solid 4He

We study the elastic properties of incommensurate solid 4He in the limit of zero temperature. Specifically, we calculate the pressure dependence of the five elastic constants (C11, C12, C13, C33, and C44), longitudinal and transversal speeds of sound, and the T = 0 Debye temperature of incommensurate and commensurate hcp 4He using the diffusion Monte Carlo method. Our results show that under compression the commensurate crystal is globally stiffer than the incommensurate, however at pressures close to melting (i.e. P ~ 25 bars) some of the elastic constants accounting for strain deformations of the hcp basal plane (C12 and C13) are slightly larger in the incommensurate solid. Also, we find that upon the introduction of tiny concentrations of point defects the shear modulus of 4He (C44) undergoes a small reduction.

preprint2013arXiv

Insights into the Phase Diagram of Bismuth Ferrite from Quasi-Harmonic Free Energy Calculations

We have used first-principles methods to investigate the phase diagram of multiferroic bismuth ferrite (BiFeO3 or BFO), revealing the energetic and vibrational features that control the occurrence of various relevant structures. More precisely, we have studied the relative stability of four low-energy BFO polymorphs by computing their free energies within the quasi-harmonic approximation, introducing a practical scheme that allows us to account for the main effects of spin disorder. As expected, we find that the ferroelectric ground state of the material (with R3c space group) transforms into an orthorhombic paraelectric phase (Pnma) upon heating. We show that this transition is not significantly affected by magnetic disorder, and that the occurrence of the Pnma structure relies on its being vibrationally (although not elastically) softer than the R3c phase. We also investigate a representative member of the family of nano-twinned polymorphs recently predicted for BFO [Prosandeev et al., Adv. Funct. Mater. 23, 234 (2013)] and discuss their possible stabilization at the boundaries separating the R3c and Pnma regions in the corresponding pressure-temperature phase diagram. Finally, we elucidate the intriguing case of the so-called super-tetragonal phases of BFO: Our results explain why such structures have never been observed in the bulk material, despite their being stable polymorphs of very low energy. Quantitative comparison with experiment is provided whenever possible, and the relative importance of various physical effects (zero-point motion, spin fluctuations, thermal expansion) and technical features (employed exchange-correlation energy density functional) is discussed. Our work attests the validity and usefulness of the quasi-harmonic scheme to investigate the phase diagram of this complex oxide, and prospective applications are discussed.

preprint2013arXiv

Possible superfluidity of molecular hydrogen in a two-dimensional crystal phase of sodium

We theoretically investigate the ground-state properties of a molecular para-hydrogen (p-H2) film in which crystallization is energetically frustrated by embedding sodium (Na) atoms periodically distributed in a triangular lattice. In order to fully deal with the quantum nature of p-H2 molecules, we employ the diffusion Monte Carlo method and realistic semi-empirical pairwise potentials describing the interactions between H2-H2 and Na-H2 species. In particular, we calculate the energetic, structural and superfluid properties of two-dimensional Na-H2 systems within a narrow density interval around equilibrium at zero temperature. In contrast to previous computational studies considering other alkali metal species such as rubidium and potassium, we find that the p-H2 ground-state is a liquid with a significantly large superfluid fraction of ~30%. The appearance of p-H2 superfluid response is due to the fact that the interactions between Na atoms and H2 molecules are less attractive than between H2 molecules. This induces a considerable reduction of the hydrogen density which favours the stabilization of the liquid phase.

preprint2013arXiv

The high-pressure high-temperature phase diagram of calcium fluoride from classical atomistic simulations

We study the phase diagram of calcium fluoride (CaF2) under pressure using classical molecular dynamic simulations performed with a simple pairwise interatomic potential of the Born-Mayer-Huggings form. Our results obtained under conditions 0 < P < 20 GPa and 0 < T < 4000 K reveal a rich variety of multi-phase boundaries involving different crystal, superionic and liquid phases, for all which we provide an accurate parametrization. Interestingly, we predict the existence of three special triple points (i.e. solid-solid-superionic, solid-superionic-superionic and superionic-superionic-liquid coexisting states) within a narrow and experimentally accessible thermodynamic range of 6 < P < 8 GPa and 1500 < T < 2750 K. Also, we examine the role of short-ranged repulsive (SR) and long-ranged attractive (LA) atomic interactions in the prediction of melting lines with the finding that SR Ca-F and LA F-F contributions are most decisive.

preprint2011arXiv

Ab initio study of compressed Ar(H2)2: structural stability and anomalous melting

We study the structural stability and dynamical properties of Ar(H2)2 under pressure using first-principles and ab initio molecular dynamics techniques. At low temperatures, Ar(H2)2 is found to stabilize in the cubic C15 Laves structure (MgCu2) and not in the hexagonal C14 Laves structure (MgZn2) as it has been assumed previously. Based on enthalpy energy and phonon calculations, we propose a temperature-induced MgCu2 -> MgZn2 phase transition that may rationalize the existing discrepancies between the sets of Raman and infrared vibron measurements. Our AIMD simulations suggest that the melting line of Ar(H2)2 presents negative slope in the interval 60 < P < 110 GPa. We explain the origin of this intriguing physical phenomenon in terms of decoupling of the Ar and H2 degrees of freedom and effective thermal-like excitations arising from coexisting liquid H2 and solid Ar phases.

preprint2011arXiv

Constraints on the phase diagram of molybdenum from first-principles free-energy calculations

We use first-principles techniques to re-examine the suggestion that transitions seen in high-P experiments on Mo are solid-solid transitions from the bcc structure to either the fcc or hcp structures. We confirm that in the harmonic approximation the free energies of fcc and hcp structures become lower than that of bcc at P > 325 GPa and T below the melting curve, as reported recently. However, we show that if anharmonic effects are fully included this is no longer true. We calculate fully anharmonic free energies of high-T crystal phases by integration of the thermal average stress with respect to strain as structures are deformed into each other, and also by thermodynamic integration from harmonic reference systems to the fully anharmonic system. Our finding that fcc is thermodynamically less stable than bcc in the relevant high-P/high-T region is supported by comparing the melting curves of the two structures calculated using the first-principles reference-coexistence technique. We present first-principles simulations based on the recently proposed Z method which also support the stability of bcc over fcc.

preprint2011arXiv

First-Principles Modeling of Pt/LaAlO3/SrTiO3 Capacitors Under an External Bias Potential

We study the electrical properties of Pt/LaAlO3/SrTiO3 capacitors under the action of an external bias potential, using first-principles simulations performed at constrained electric displacement field. A complete set of band diagrams, together with the relevant electrical characteristics (capacitance and built-in fields), are determined as a function of LaAlO3 thickness and the applied potential.We find that the internal field in LaAlO3 monotonically decreases with increasing thickness; hence, the occurrence of spontaneous Zener tunneling is ruled out in this system.We discuss the implications of our results in the light of recent experimental observations on biased LaAlO3/SrTiO3 junctions involving metallic top electrodes.