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

Carlos P. Herrero

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

ResearcherAffiliation not importedOpen to collaborate
cond-mat.mtrl-sciphysics.chem-phcond-mat.stat-mechcond-mat.dis-nn

Trust snapshot

Quick read

Trust 21 - Emerging
18works
0followers
4topics
2close 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

18 published item(s)

preprint2020arXiv

Isotopic effects in structural properties of graphene

Isotopic effects are relevant to understand several properties of solids, and have been thoroughly analyzed along the years. These effects may depend on the dimensionality of the considered solid. Here we assess their magnitude for structural properties of graphene, a paradigmatic two-dimensional material. We use path-integral molecular dynamics simulations, a well-suited technique to quantify the influence of nuclear quantum effects on equilibrium variables, especially in cases where anharmonic effects are important. Emphasis is put on interatomic distances and mean-square displacements, as well as on the in-plane area of the graphene layer. At low temperature, the relative difference in C--C distance for $^{13}$C and $^{14}$C, with respect to $^{12}$C, is found to be $-2.5$ and $-4.6 \times 10^{-4}$, respectively, larger than in three-dimensional carbon-based materials such as diamond. For the in-plane area, the relative changes amount to $-3.9$ and $-6.9 \times 10^{-4}$. The magnitude of anharmonicity in the lattice vibrations is estimated by comparing the internal energy and atomic vibrational amplitudes with those derived from a harmonic approximation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Nuclear quantum effects in graphane

Graphane is a quasi-two-dimensional material consisting of a single layer of fully hydrogenated graphene, with a C:H ratio of 1. We study nuclear quantum effects in the so-called chair-graphane by using path-integral molecular dynamics (PIMD) simulations. The interatomic interactions are modeled by a tight-binding potential model fitted to density-functional calculations. Finite-temperature properties are studied in the range from 50 to 1500~K. To assess the magnitude of nuclear quantum effects in the properties of graphane, classical molecular dynamics simulations have been also performed. These quantum effects are significant in structural properties such as interatomic distances and layer area at finite temperatures. The in-plane compressibility of graphane is found to be about twice larger than that of graphene, and at low temperature it is 9\% higher than the classical calculation. The thermal expansion coefficient resulting from PIMD simulations vanishes in the zero-temperature limit, in agreement with the third law of Thermodynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Thermodynamic properties of graphene bilayers

Thermodynamic properties of graphene bilayers are studied by path-integral molecular dynamics (PIMD) simulations, considering quantization of vibrational modes and anharmonic effects. Bilayer graphene has been studied at temperatures between 12 and 1500~K for zero external stress, using the LCBOPII effective potential. We concentrate on the thermal expansion, in-plane and out-of-plane compressibility, and specific heat. Additional insight into the meaning of our results for bilayer graphene is obtained from a comparison with data obtained from PIMD simulations for monolayer graphene and graphite. They are also analyzed in view of experimental data for graphite. Zero-point and thermal effects on the in-plane and "real" area of bilayer graphene are studied. The thermal expansion coefficient $α_{xy}$ of the in-plane area is negative at low temperatures and positive for $T \gtrsim$ 800~K. The minimum $α_{xy}$ is $-6.6 \times 10^{-6}$ K$^{-1}$ at $T \approx 220$~K. Both in-plane ($χ_{xy}$) and out-of-plane ($χ_z$) compressibilities of graphene bilayers are found to increase for rising temperature, and turn out to be lower than that corresponding to monolayer graphene and higher than those found for graphite. At 300 K, we find for the bilayer $χ_{xy} = 9.5 \times 10^{-2}$ Å$^2$/eV and $χ_z = 2.97 \times 10^{-2}$ GPa$^{-1}$. Results for the specific heat obtained from the simulations are compared with those given by a harmonic approximation for the vibrational modes. This approach is noticeably accurate at temperatures lower than 200~K.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Quantum effects in graphene monolayers: Path-integral simulations

Path-integral molecular dynamics (PIMD) simulations have been carried out to study the influence of quantum dynamics of carbon atoms on the properties of a single graphene layer. Finite-temperature properties were analyzed in the range from 12 to 2000~K, by using the LCBOPII effective potential. To assess the magnitude of quantum effects in structural and thermodynamic properties of graphene, classical molecular dynamics simulations have been also performed. Particular emphasis has been laid on the atomic vibrations along the out-of-plane direction. Even though quantum effects are present in these vibrational modes, we show that at any finite temperature classical-like motion dominates over quantum delocalization, provided that the system size is large enough. Vibrational modes display an appreciable anharmonicity, as derived from a comparison between kinetic and potential energy of the carbon atoms. Nuclear quantum effects are found to be appreciable in the interatomic distance and layer area at finite temperatures. The thermal expansion coefficient resulting from PIMD simulations vanishes in the zero-temperature limit, in agreement with the third law of thermodynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Ising model in clustered scale-free networks

The Ising model in clustered scale-free networks has been studied by Monte Carlo simulations. These networks are characterized by a degree distribution of the form P(k) ~ k^(-gamma) for large k. Clustering is introduced in the networks by inserting triangles, i.e., triads of connected nodes. The transition from a ferromagnetic (FM) to a paramagnetic (PM) phase has been studied as a function of the exponent gamma and the triangle density. For gamma > 3 our results are in line with earlier simulations, and a phase transition appears at a temperature T_c(gamma) in the thermodynamic limit (system size N to infinity). For gamma <= 3, a FM-PM crossover appears at a size-dependent temperature T_co, so that the system remains in a FM state at any finite temperature in the limit N to infinity. Thus, for gamma = 3, T_co scales as ln N, whereas for gamma < 3, we find T_co ~ J N^z, where the exponent z decreases for increasing gamma. Adding motifs (triangles in our case) to the networks causes an increase in the transition (or crossover) temperature for exponent gamma > 3 (or <= 3). For gamma > 3, this increase is due to changes in the mean values < k > and < k^2 >, i.e., the transition is controlled by the degree distribution (nearest neighbor connectivities). For gamma <= 3, however, we find that clustered and unclustered networks with the same size and distribution P(k) have different crossover temperature, i.e., clustering favors FM correlations, thus increasing the temperature T_co. The effect of a degree cutoff k_cut on the asymptotic behavior of T_co is discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Path-integral simulation of ice VII: Pressure and temperature effects

Effects of pressure and temperature on structural and thermodynamic properties of ice VII have been studied by using path-integral molecular dynamics (PIMD) simulations. Temperatures between 25 and 450 K, as well as pressures up to 12 GPa were considered. Interatomic interactions were modeled by using the effective q-TIP4P/F potential for flexible water. We analyze the pressure dependence of the molar volume, bulk modulus, interatomic distances, kinetic energy, and atomic delocalization at various temperatures. Results of PIMD simulations are compared with those derived from a quasi-harmonic approximation (QHA) of vibrational modes, which helps to assess the importance of anharmonic effects, as well as the influence of the different modes on the properties of ice VII. The accuracy of the QHA for describing this high-pressure phase decreases for rising temperature, but this approximation becomes more reliable as pressure grows, since anharmonicity becomes less relevant. Comparisons with low-pressure cubic ice are presented.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Configurational entropy of hydrogen-disordered ice polymorphs

The configurational entropy of several H-disordered ice polymorphs is calculated by means of a thermodynamic integration along a path between a totally H-disordered state and one fulfilling the Bernal-Fowler ice rules. A Monte Carlo procedure based on a simple energy model is used, so that the employed thermodynamic path drives the system from high temperatures to the low-temperature limit. This method turns out to be precise enough to give reliable values for the configurational entropy of different ice phases in the thermodynamic limit (number of molecules N --> infinity). The precision of the method is checked for the ice model on a two-dimensional square lattice. Results for the configurational entropy are given for H-disordered arrangements on several polymorphs, including ices Ih, Ic, II, III, IV, V, VI, and XII. The highest and lowest entropy values correspond to ices VI and XII, respectively, with a difference of 3.3\% between them. The dependence of the entropy on the ice structures has been rationalized by comparing it with structural parameters of the various polymorphs, such as the mean ring size. A particularly good correlation has been found between the configurational entropy and the connective constant derived from self-avoiding walks on the ice networks.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Path-integral simulation of solids

The path-integral formulation of the statistical mechanics of quantum many-body systems is described, with the purpose of introducing practicaltechniques for the simulation of solids. Monte Carlo and molecular dynamics methods for distinguishable quantum particles are presented, with particular attention to the isothermal-isobaric ensemble. Applications of these computational techniques to different types of solids are reviewed, including noble-gas solids (helium and heavier elements), group-IV materials (diamond and elemental semiconductors), and molecular solids (with emphasis on hydrogen and ice). Structural, vibrational, and thermodynamic properties of these materials are discussed. Applications also include point defects in solids (structure and diffusion), as well as nuclear quantum effects in solid surfaces and adsorbates. Different phenomena are discussed, as solid-to-solid and orientational phase transitions, rates of quantum processes, classical-to-quantum crossover, and various finite-temperature anharmonic effects (thermal expansion, isotopic effects, electron-phonon interactions). Nuclear quantum effects are most remarkable in the presence of light atoms, so that especial emphasis is laid on solids containing hydrogen as a constituent element or as an impurity.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Structural characterization of ice polymorphs from self-avoiding walks

Topological properties of crystalline ice structures are studied by means of self-avoiding walks on their H-bond networks. The number of self-avoiding walks, C_n, for eight ice polymorphs has been obtained by direct enumeration up to walk length n = 27. This has allowed us to determine the `connective constant&#39; or effective coordination number `mu&#39; of these structures as the limit of the ratio C_n/C_{n-1} for large n. This structure-dependent parameter `mu&#39; is related with other topological characteristics of ice polymorphs, such as the mean and minimum ring size, or the topological density of network sites. A correlation between the connective constant and the configurational entropy of hydrogen-disordered ice structures is discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Configurational entropy of ice from thermodynamic integration

The configurational entropy of ice is calculated by thermodynamic integration from high to low temperatures. We use Monte Carlo simulations with a simple energy model which reproduces the Bernal-Fowler ice rules. This procedure is found to be precise enough to give reliable values for the residual entropy s_th of different ice phases in the thermodynamic limit. First, we check it for a two-dimensional ice model. Second, we calculate s_th for ice Ih, and compare our result with those previously given in the literature. Third, we obtain s_th for ice VI, for which we find a value clearly higher than for ice Ih.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

High-density amorphous ice: A path-integral simulation

Structural and thermodynamic properties of high-density amorphous (HDA) ice have been studied by path-integral molecular dynamics simulations in the isothermal-isobaric ensemble. Interatomic interactions were modeled by using the effective q-TIP4P/F potential for flexible water. Quantum nuclear motion is found to affect several observable properties of the amorphous solid. At low temperature (T = 50 K) the molar volume of HDA ice is found to increase by 6%, and the intramolecular O--H distance rises by 1.4% due to quantum motion. Peaks in the radial distribution function of HDA ice are broadened respect to their classical expectancy. The bulk modulus, B, is found to rise linearly with the pressure, with a slope dB/dP = 7.1. Our results are compared with those derived earlier from classical and path-integral simulations of HDA ice. We discuss similarities and discrepancies with those earlier simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Path-integral simulation of ice Ih: The effect of pressure

The effect of pressure on structural and thermodynamic properties of ice Ih has been studied by means of path-integral molecular dynamics simulations at temperatures between 50 and 300 K. Interatomic interactions were modeled by using the effective q-TIP4P/F potential for flexible water. Positive (compression) and negative (tension) pressures have been considered, which allowed us to approach the limits for the mechanical stability of this solid water phase. We have studied the pressure dependence of the crystal volume, bulk modulus, interatomic distances, atomic delocalization, and kinetic energy. The spinodal point at both negative and positive pressures is derived from the vanishing of the bulk modulus. For P < 0, the spinodal pressure changes from -1.38 to -0.73 GPa in the range from 50 to 300 K. At positive pressure the spinodal is associated to ice amorphization, and at low temperatures is found between 1.1 and 1.3 GPa. Quantum nuclear effects cause a reduction of the metastability region of ice Ih.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Diffusion of hydrogen in graphite: A molecular dynamics simulation

Diffusion of atomic and molecular hydrogen in the interstitial space between graphite sheets has been studied by molecular dynamics simulations. Interatomic interactions were modeled by a tight-binding potential fitted to density-functional calculations. Atomic hydrogen is found to be bounded to C atoms, and its diffusion consists in jumping from a C atom to a neighboring one, with an activation energy of about 0.4 eV. Molecular hydrogen is less attached to the host sheets and diffuses faster than isolated H. At temperatures lower than 500 K, H_2 diffuses with an activation energy of 89 meV, whereas at higher T its diffusion is enhanced by longer jumps of the molecule as well as by correlations between successive hops, yielding an effective activation energy of 190 meV.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Isotope effects in ice Ih: A path-integral simulation

Ice Ih has been studied by path-integral molecular dynamics simulations, using the effective q-TIP4P/F potential model for flexible water. This has allowed us to analyze finite-temperature quantum effects in this solid phase from 25 to 300 K at ambient pressure. Among these effects we find a negative thermal expansion of ice at low temperatures, which does not appear in classical molecular dynamics simulations. The compressibility derived from volume fluctuations gives results in line with experimental data. We have analyzed isotope effects in ice Ih by considering normal, heavy, and tritiated water. In particular, we studied the effect of changing the isotopic mass of hydrogen on the kinetic energy and atomic delocalization in the crystal, as well as on structural properties such as interatomic distances and molar volume. For D$_2$O ice Ih at 100 K we obtained a decrease in molar volume and intramolecular O--H distance of 0.6% and 0.4%, respectively, as compared to H$_2$O ice.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Molecular hydrogen in graphite: A path-integral simulation

Molecular hydrogen in the bulk of graphite has been studied by path-integral molecular dynamics simulations. Finite-temperature properties of H_2 molecules adsorbed between graphite layers were analyzed in the temperature range from 300 to 900 K. The interatomic interactions were modeled by a tight-binding potential fitted to density-functional calculations. In the lowest-energy position, an H_2 molecule is found to be disposed parallel to the sheets plane. At finite temperatures, the molecule explores other orientations, but its rotation is partially hindered by the adjacent graphite layers. Vibrational frequencies were obtained from a linear-response approach, based on correlations of atom displacements. For the stretching vibration of the molecule, we find at 300 K a frequency omega_s = 3916 cm-1, more than 100 cm-1 lower than the frequency corresponding to an isolated H_2 molecule. Isotope effects have been studied by considering also deuterium and tritium molecules. For D_2 in graphite we obtained omega_s = 2816 cm-1}, i.e., an isotopic ratio omega_s(H) / omega_s(D) = 1.39.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Isotope effects on the lattice parameter of cubic SiC

Path-integral molecular dynamics simulations in the isothermal-isobaric (NPT) ensemble have been carried out to study the dependence of the lattice parameter of 3C-SiC upon isotope mass. This computational method allows a quantitative and nonperturbative study of such anharmonic effect. Atomic nuclei were treated as quantum particles interacting via a tight-binding-type potential. At 300 K, the difference Delta a between lattice parameters of 3C-SiC crystals with 12C and 13C amounts to 2.1 x 10^{-4} A. The effect due to Si isotopes is smaller, and amounts to 3.5 x 10^{-5} A when replacing 28Si by 29Si. Results of the PIMD simulations are interpreted in terms of a quasiharmonic approximation for the lattice vibrations.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Vibrational properties and diffusion of hydrogen on graphene

Hydrogen and deuterium chemisorption on a single layer of graphene has been studied by path-integral molecular dynamics simulations. Finite-temperature properties of these point defects were analyzed in the range from 200 to 1500 K, by using a tight-binding potential fitted to density-functional calculations. On one side, vibrational properties of the adatoms are studied at their equilibrium positions, linked to C atoms. The vibrations display an appreciable anharmonicity, as derived from comparison between kinetic and potential energy, as well as between vibrational energy for hydrogen and deuterium. On the other side, adatom motion has been studied by quantum transition-state theory. At room temperature, quantum effects are found to enhance the hydrogen diffusivity on the graphene sheet by a factor of 20.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Solid helium at high pressure: A path-integral Monte Carlo simulation

Solid helium (3He and 4He) in the hcp and fcc phases has been studied by path-integral Monte Carlo. Simulations were carried out in the isothermal-isobaric (NPT) ensemble at pressures up to 52 GPa. This allows one to study the temperature and pressure dependences of isotopic effects on the crystal volume and vibrational energy in a wide parameter range. The obtained equation of state at room temperature agrees with available experimental data. The kinetic energy, E_k, of solid helium is found to be larger than the vibrational potential energy, E_p. The ratio E_k/E_p amounts to about 1.4 at low pressures, and decreases as the applied pressure is raised, converging to 1, as in a harmonic solid. Results of these simulations have been compared with those yielded by previous path integral simulations in the NVT ensemble. The validity range of earlier approximations is discussed.

0 citations0 reviewsNo code linkedOpen access