The reconstruction of Rh(001) upon oxygen adsorption

We report on a first-principles study of the structure of O/Rh(001) at half a monolayer of oxygen coverage, performed using density-functional theory. We find that oxygen atoms sit at the center of the black squares of a chess-board, $c(2\times 2)$, pattern. This structure is unstable against a rhomboid distortion of the black squares, which shortens the distance between an O atom and two of the four neighboring Rh atoms, while lengthening the distance with respect to the other two. We actually find that the surface energy is further lowered by allowing the O atom to get off the short diagonal of the rhombus so formed. We predict that the latter distortion is associated with an order-disorder transition, occurring below room temperature. The above rhomboid distortion of the square lattice may be seen as a rotation of the empty, white, squares. Our findings are at variance with recent claims based on STM images, according to which it is instead the black squares which would rotate. We argue that these images are indeed compatible with our predicted reconstruction pattern.Comment: 14 pages (inclusive of 5 figures). To appear on Surface Scienc

Rotational dynamics of CO solvated in small He clusters: a quantum Monte Carlo study

The rotational dynamics of CO single molecules solvated in small He clusters (CO@He_N) has been studied using Reptation Quantum Monte Carlo for cluster sizes up to N=30. Our results are in good agreement with the roto-vibrational features of the infrared spectrum recently determined for this system, and provide a deep insight into the relation between the structure of the cluster and its dynamics. Simulations for large N also provide a prediction of the effective moment of inertia of CO in the He nano-droplet regime, which has not been measured so far

Reptation quantum Monte Carlo for lattice Hamiltonians with a directed-update scheme

We provide an extension to lattice systems of the reptation quantum Monte Carlo algorithm, originally devised for continuous Hamiltonians. For systems affected by the sign problem, a method to systematically improve upon the so-called fixed-node approximation is also proposed. The generality of the method, which also takes advantage of a canonical worm algorithm scheme to measure off-diagonal observables, makes it applicable to a vast variety of quantum systems and eases the study of their ground-state and excited-states properties. As a case study, we investigate the quantum dynamics of the one-dimensional Heisenberg model and we provide accurate estimates of the ground-state energy of the two-dimensional fermionic Hubbard model

The itinerant ferromagnetic phase of the Hubbard model

Using a newly developed quantum Monte Carlo technique, we provide strong evidence for the stability of a saturated ferromagnetic phase in the high-density regime of the two-dimensional infinite-U Hubbard model. By decreasing the electron density, a discontinuous transition to a paramagnetic phase is observed, accompanied by a divergence of the susceptibility on the paramagnetic side. This behavior, resulting from a high degeneracy among different spin sectors, is consistent with an infinite-order phase transition. The remarkable stability of itinerant ferromagnetism renews the hope to describe this phenomenon within a purely kinetic mechanism and will facilitate the validation of experimental quantum simulators with cold atoms loaded in optical lattices

Cross-sectional imaging of sharp Si interlayers embedded in gallium arsenide

We investigate the electronic properties of the (110) cross-sectional surface of Si-doped GaAs using first-principles techniques. We focus on doping configurations with an equal concentration of Si impurities in cationic and anionic sites, such as occurring in a self-compensating doping regime. In particular we study a bilayer of Si atoms uniformly distributed over two consecutive (001) atomic layers. The simulated cross-sectional scanning tunneling microscopy images show a bright signal at negative bias, which is strongly attenuated when the bias is reversed. This scenario is consistent with experimental results which had been attributed to hitherto unidentified Si complexes.Comment: 10 pages, 3 figure

Structure, rotational dynamics, and superfluidity of small OCS-doped He clusters

The structural and dynamical properties of OCS molecules solvated in Helium clusters are studied using reptation quantum Monte Carlo, for cluster sizes n=3-20 He atoms. Computer simulations allow us to establish a relation between the rotational spectrum of the solvated molecule and the structure of the He solvent, and of both with the onset of superfluidity. Our results agree with a recent spectroscopic study of this system, and provide a more complex and detailed microscopic picture of this system than inferred from experiments.Comment: 4 pages. TeX (requires revtex4) + 3 ps figures (1 color

Seebeck coefficient of liquid water from equilibrium molecular dynamics

The application of a temperature gradient to an extended system generates an electromotive force that induces an electric current in conductors and a macroscopic polarization in insulators. The ratio of the electromotive force to the temperature difference, usually referred to as the Seebeck coefficient, is often computed using non-equilibrium techniques, such as non-equilibrium molecular dynamics (NEMD). In this paper we argue that thermo-polarization effects in insulating fluids can be conveniently treated by standard equilibrium thermodynamics and devise a protocol, based on a combination of equilibrium molecular dynamics and Bayesian inference methods, that allows one to compute the Seebeck coefficient in these systems along with a rigorous estimate of the resulting statistical accuracy. The application of our methodology to liquid SPC/E water results in good agreement with previous studies, based on more elaborate NEMD simulations, and in a more reliable estimate of the statistical accuracy of the results.Comment: 6 pages, 5 figure