Carlos P. Herrero

Ising model in small-world networks.

The Ising model in small-world networks generated from two- and three-dimensional regular lattices has been studied. Monte Carlo simulations were carried out to characterize the ferromagnetic transition appearing in these systems. In the thermodynamic limit, the phase transition has a mean-field character for any finite value of the rewiring probability p, which measures the disorder strength of a given network. For small values of p, both the transition temperature and critical energy change with p as a power law. In the limit p-->0, the heat capacity at the transition temperature diverges logarithmically in two-dimensional (2D) networks and as a power law in 3D.

Anomalous nuclear quantum effects in ice

One striking anomaly of water ice has been largely neglected and never explained. Replacing hydrogen (1H) by deuterium (2H) 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 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 16O by 18O causes a normal lattice contraction.

Self-avoiding walks on scale-free networks.

Several kinds of walks on complex networks are currently used to analyze search and navigation in different systems. Many analytical and computational results are known for random walks on such networks. Self-avoiding walks (SAWs) are expected to be more suitable than unrestricted random walks to explore various kinds of real-life networks. Here we study long-range properties of random SAWs on scale-free networks, characterized by a degree distribution P(k) approximately k(-gamma). In the limit of large networks (system size N-->infinity), the average number sn of SAWs starting from a generic site increases as mu(n) , with mu = k2/k-1 . For finite N, sn is reduced due to the presence of loops in the network, which causes the emergence of attrition of the paths. For kinetic growth walks, the average maximum length L increases as a power of the system size: L approximately Nalpha, with an exponent alpha increasing as the parameter gamma is raised. We discuss the dependence of alpha on the minimum allowed degree in the network. A similar power-law dependence is found for the mean self-intersection length of nonreversal random walks. Simulation results support our approximate analytical calculations.

Ising model in scale-free networks: a Monte Carlo simulation.

The Ising model in uncorrelated scale-free networks has been studied by means of Monte Carlo simulations. These networks are characterized by a degree (or connectivity) distribution P (k) approximately k(-gamma). The ferromagnetic-paramagnetic transition temperature has been studied as a function of the parameter gamma. For gamma>3 our results agree with earlier analytical calculations, which found a phase transition at a temperature T(c) (gamma) in the thermodynamic limit. For gamma< or =3, a ferromagnetic-paramagnetic crossover occurs at a size-dependent temperature T(co), and the system is in the ordered ferromagnetic state at any temperature for a system size N--> infinity. For gamma=3 and large enough N, the crossover temperature is found to be T(co) approximately A ln N, with a prefactor A proportional to the mean degree. For 2N(z), with an exponent z that decreases as gamma increases. This exponent is found to be lower than predicted by earlier calculations.

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 practical techniques 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.

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

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.

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.

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.

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 of 17-50 K and pressures between 10(-2) and 2 x 10(3) 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 nonequilibrium 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 toward 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 (20)Ne, (21)Ne, (22)Ne, and natural neon. The simulation data show satisfactory agreement to previous experimental results, which report a shift of about 0.15 K between triple-point temperatures of (20)Ne and (22)Ne. 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.