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

C. P. Herrero

C. P. Herrero contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.mtrl-sciphysics.chem-phcond-mat.mes-hall

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
3topics
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

9 published item(s)

preprint2020arXiv

Anharmonicity of the acoustic modes of graphene

The anharmonicity of the acoustic phonon dispersion of graphene has been studied by the harmonic linear response (HLR) approach at finite temperature. This is a non-perturbative method based on the linear response of the system to applied forces, as derived from equilibrium computer simulations. Anharmonic shifts are analyzed in the long-wavelength limit at room temperature, with emphasis in the effect of applied tensile or compressive in-plane stress. The simulation results are compared with available analytical models, based either on first-order perturbation theory or on a description by anomalous exponents. The simulations show better agreement to the expectations of the perturbational approach. The effect of temperature and zero-point vibrations on the acoustic out-of-plane anharmonic shifts of graphene are briefly reviewed.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Phonon dispersion in two-dimensional solids from atomic probability distributions

We propose a harmonic linear response (HLR) method to calculate the phonon dispersion relations of two-dimensional (2D) layers from equilibrium simulations at finite temperature. This HLR approach is based on the linear response of the system, as derived from the analysis of its centroid density in equilibrium path integral simulations. In the classical limit, this approach is closely related to those methods that study vibrational properties by the diagonalization of the covariance matrix of atomic fluctuations. The validity of the method is tested in the calculation of the phonon dispersion relations of a graphene monolayer, a graphene bilayer, and graphane. Anharmonic effects in the phonon dispersion relations of graphene are demonstrated by the calculation of the temperature dependence of the following observables: the kinetic energy of the carbon atoms, the vibrational frequency of the optical $E_{2g}$ mode, and the elastic moduli of the layer.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The phase diagram of ice: a quasi-harmonic study based on a flexible water model

The phase diagram of ice is studied by a quasi-harmonic approximation. The free energy of all experimentally known ice phases has been calculated with the flexible q-TIP4P/F model of water. The only exception is the high pressure ice X, in which the presence of symmetric O-H-O bonds prevents its modeling with this empirical interatomic potential. The simplicity of our approach allows us to study ice phases at state points of the T-P plane that have been omitted in previous simulations using free energy methods based on thermodynamic integration. The effect in the phase diagram of averaging the proton disorder that appears in several ice phases has been studied. It is found particularly relevant for ice III, at least for cell sizes typically used in phase coexistence simulations. New insight into the capability of the employed water model to describe the coexistence of ice phases is presented. We find that the H-ordered ices IX and XIV, as well as the H-disordered ice XII, are particularly stable for this water model. This fact disagrees with experimental data. The unexpected large stability of ice IX is a property related to the TIP4P-character of the water model. Only after omission of these three stable ice phases, the calculated phase diagram becomes in reasonable qualitative agreement to the experimental one in the T-P region corresponding to ice Ih, II, III, V, and VI. The calculation of the phase diagram in the quantum and classical limits shows that the most important quantum effect is the stabilization of ice II due to its lower zero-point energy when compared to that one of ices Ih, III, and V.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Anomalous Nuclear Quantum Effects in Ice

One striking anomaly of water ice has been largely neglected and never explained. Replacing hydrogen ($^1$H) by deuterium ($^2$H) causes ice to expand, whereas the "normal" isotope effect is volume contraction with increased mass. Furthermore, the anomaly increases with temperature $T$, even though a normal isotope shift should decrease with $T$ and vanish when $T$ is high enough to use classical nuclear motions. In this study, we show that these effects are very well described by {\it ab initio} density functional theory. Our theoretical modeling explains these anomalies, and allows us to predict and to experimentally confirm a counter effect, namely that replacement of $^{16}$O by $^{18}$O causes a normal lattice contraction.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Quasi-harmonic approximation of thermodynamic properties of ice Ih, II, and III

Several thermodynamic properties of ice Ih, II, and III are studied by a quasi-harmonic approximation and compared to results of quantum path integral and classical simulations. This approximation allows to obtain thermodynamic information at a fraction of the computational cost of standard simulation methods, and at the same time permits studying quantum effects related to zero point vibrations of the atoms. Specifically we have studied the crystal volume, bulk modulus, kinetic energy, enthalpy and heat capacity of the three ice phases as a function of temperature and pressure. The flexible q-TIP4P/F model of water was employed for this study, although the results concerning the capability of the quasi-harmonic approximation are expected to be valid independently of the employed water model. The quasi-harmonic approximation reproduces with reasonable accuracy the results of quantum and classical simulations showing an improved agreement at low temperatures (T < 100 K). This agreement does not deteriorate as a function of pressure as long as it is not too close to the limit of mechanical stability of the ice phases.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The phase diagram of ice Ih, II, and III: a quasi-harmonic study

The phase diagram of ice Ih, II, and III is studied by a quasi-harmonic approximation. The results of this approach are compared to phase diagrams previously derived by thermodynamic integration using path integral and classical simulations, as well as to experimental data. The studied models are based on both flexible (q-TIP4P/F) and rigid (TIP4P/2005, TIP4PQ/2005) descriptions of the water molecule. Many aspects of the simulated phase diagrams are reasonably reproduced by the quasi-harmonic approximation. Advantages of this simple approach are that it is free from the statistical errors inherent to computer simulations, both classical and quantum limits are easily accessible, and the error of the approximation is expected to decrease in the zero temperature limit. We find that the calculated phase diagram of ice Ih, II, and III depends strongly on the hydrogen disorder of ice III, at least for cell sizes typically used in phase coexistence simulations. Either ice II (in the classical limit) or ice III (in the quantum one) may become unstable depending upon the proton disorder in ice III. The comparison of quantum and classical limits shows that the stabilization of ice II is the most important quantum effect in the phase diagram. The lower vibrational zero-point energy of ice II, compared to either ice Ih or III, is the microscopic origin of this stabilization. The necessity of performing an average of the lattice energy over the proton disorder of ice III is discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Kinetic energy of protons in ice Ih and water: a path integral study

The kinetic energy of H and O nuclei has been studied by path integral molecular dynamics simulations of ice Ih and water at ambient pressure. The simulations were performed by using the q-TIP4P/F model, a point charge empirical potential that includes molecular flexibility and anharmonicity in the OH stretch of the water molecule. Ice Ih was studied in a temperature range between 210-290 K, and water between 230-320 K. Simulations of an isolated water molecule were performed in the range 210-320 K to estimate the contribution of the intramolecular vibrational modes to the kinetic energy. Our results for the proton kinetic energy, K_H, in water and ice Ih show both agreement and discrepancies with different published data based on deep inelastic neutron scattering experiments. Agreement is found for water at the experimental melting point and in the range 290-300 K. Discrepancies arise because data derived from the scattering experiments predict in water two maxima of K_H around 270 K and 277 K, and that K_H is lower in ice than in water at 269 K. As a check of the validity of the employed water potential, we show that our simulations are consistent with other experimental thermodynamic properties related to K_H, as the temperature dependence of the liquid density, the heat capacity of water and ice at constant pressure, and the isotopic shift in the melting temperature of ice upon isotopic substitution of either H or O atoms. Moreover, the temperature dependence of K_H predicted by the q-TIP4P/F model for ice Ih is found to be in good agreement to results of path integral simulations using ab initio density functional theory.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Quantum path integral simulation of isotope effects in the melting temperature of ice Ih

The isotope effect in the melting temperature of ice Ih has been studied by free energy calculations within the path integral formulation of statistical mechanics. Free energy differences between isotopes are related to the dependence of their kinetic energy on the isotope mass. The water simulations were performed by using the q-TIP4P/F model, a point charge empirical potential that includes molecular flexibility and anharmonicity in the OH stretch of the water molecule. The reported melting temperature at ambient pressure of this model (T = 251 K) increases by 6.5+-0.5 K and 8.2+-0.5 K upon isotopic substitution of hydrogen by deuterium and tritium, respectively. These temperature shifts are larger than the experimental ones (3.8 K and 4.5 K, respectively). In the classical limit, the melting temperature is nearly the same as that for tritiated ice. This unexpected behavior is rationalized by the coupling between intermolecular interactions and molecular flexibility. This coupling makes the kinetic energy of the OH stretching modes larger in the liquid than in the solid phase. However the opposite behavior is found for intramolecular modes, which display larger kinetic energy in ice than in liquid water.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Quantum path-integral study of the phase diagram and isotope effects of neon

The phase diagram of natural neon has been calculated for temperatures in the range 17-50 K and pressures between 0.01 and 2000 bar. The phase coexistence between solid, liquid, and gas phases has been determined by the calculation of the separate free energy of each phase as a function of temperature. Thus, for a given pressure, the coexistence temperature was obtained by the condition of equal free energy of coexisting phases. The free energy was calculated by using non-equilibrium techniques such as adiabatic switching and reversible scaling. The phase diagram obtained by classical Monte Carlo simulations has been compared to that obtained by quantum path-integral simulations. Quantum effects related to the finite mass of neon cause that coexistence lines are shifted towards lower temperatures when compared to the classical limit. The shift found in the triple point amounts to 1.5 K, i.e., about 6 % of the triple-point temperature. The triple-point isotope effect has been determined for 20Ne, 21Ne, 22Ne, and natural neon. The simulation data show satisfactory agreement to previous experimental results, that report a shift of about 0.15 K between triple-point temperatures of 20Ne and 22Ne. The vapor-pressure isotope effect has been calculated for both solid and liquid phases at triple-point conditions. The quantum simulations predict that this isotope effect is larger in the solid than in the liquid phase, and the calculated values show nearly quantitative agreement to available experimental data.

0 citations0 reviewsNo code linkedOpen access