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

Rafael Ramirez

Rafael Ramirez contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.mtrl-sciphysics.chem-phcond-mat.dis-nncond-mat.stat-mechmath-phmath.DSmath.MP

Trust snapshot

Quick read

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

16 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

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
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

Cartesian approach for constrained mechanical system with three degree of freedom

In the history of mechanics, there have been two points of view for studying mechanical systems: The Newtonian and the Cartesian. According the Descartes point of view, the motion of mechanical systems is described by the first-order differential equations in the $N$ dimensional configuration space $\textsc{Q}.$ In this paper we develop the Cartesian approach for mechanical systems with three degrees of freedom and with constraint which are linear with respect to velocity. The obtained results we apply to discuss the integrability of the geodesic flows on the surface in the three dimensional Euclidian space and to analyze the integrability of a heavy rigid body in the Suslov and the Veselov cases .

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