Path Integral Molecular Dynamics of Liquid Water in a Mean-Field Particle Reservoir

We present a simulation scheme for path integral simulation of molecular liquids where a small open region is embedded in a large reservoir of non interacting point-particles. The scheme is based on the latest development of the adaptive resolution technique AdResS and allows for the space-dependent change of molecular resolution from a path integral representation with 120 degrees of freedom to a point particle that does not interact with other molecules and vice versa. The method is applied to liquid water and implies a sizable gain regarding the request of computational resources compared to full path integral simulations. Given the role of water as universal solvent with a specific hydrogen bonding network, the path integral treatment of water molecules is important to describe the quantum effects of hydrogen atoms’ delocalization in space on the hydrogen bonding network. The method presented here implies feasible computational efforts compared to full path integral simulations of liquid water which, on large scales, are often prohibitive

Artificial Neural Network Methods in Quantum Mechanics

In a previous article we have shown how one can employ Artificial Neural Networks (ANNs) in order to solve non-homogeneous ordinary and partial differential equations. In the present work we consider the solution of eigenvalue problems for differential and integrodifferential operators, using ANNs. We start by considering the Schr\"odinger equation for the Morse potential that has an analytically known solution, to test the accuracy of the method. We then proceed with the Schr\"odinger and the Dirac equations for a muonic atom, as well as with a non-local Schr\"odinger integrodifferential equation that models the $n+\alpha$ system in the framework of the resonating group method. In two dimensions we consider the well studied Henon-Heiles Hamiltonian and in three dimensions the model problem of three coupled anharmonic oscillators. The method in all of the treated cases proved to be highly accurate, robust and efficient. Hence it is a promising tool for tackling problems of higher complexity and dimensionality.Comment: Latex file, 29pages, 11 psfigs, submitted in CP

Electronic Origin of α″ to β Phase Transformation in Ti-Nb-Based Thin Films upon Hf Microalloying.

We present results on thin Ti-Nb-based films containing Hf at various concentrations grown by magnetron sputtering. The films exhibit α" patterns at Hf concentrations up to 11 at.%, while at 16 at.% Hf, the β-phase emerges as a stable structure. These findings were consolidated by ab initio calculations, according to which the α"-β transformation is manifested in the calculation of the electronic band energies for Hf contents between 11 and 18 at.%. It turns out that the β-phase transition originates from the Hf 5d contributions at the Fermi level and the Hf 6s hybridizations at low energies in the electronic density of states. Bonding-anti-bonding first neighbor features existing in the shifted plane destabilize the α″-phase, especially at high Hf concentrations, while the covalent-like features in the first neighborhood stabilize the corresponding plane of the β-phase. Thin films measurements and bulk total energy calculations agree that the lattice constants of both α″ and β phases increase upon Hf substitution. These results are important for the understanding of β-Ti-based alloys formation mechanisms and can be used for the design of suitable biocompatible materials

Immersed nano-sized Al dispersoids in an Al matrix; effects on the structural and mechanical properties by Molecular Dynamics simulations

We used molecular dynamics simulations based on a potential model in analogy to the Tight Binding scheme in the Second Moment Approximation to simulate the effects of aluminum icosahedral grains (dispersoids) on the structure and the mechanical properties of an aluminum matrix. First we validated our model by calculating several thermodynamic properties referring to the bulk Al case and we found good agreement with available experimental and theoretical data. Afterwards, we simulated Al systems containing Al clusters of various sizes. We found that the structure of the Al matrix is affected by the presence of the dispersoids resulting in well ordered domains of different symmetries that were identified using suitable Voronoi analysis. In addition, we found that the increase of the grain size has negative effect on the mechanical properties of the nanocomposite as manifested by the lowering of the calculated bulk moduli. The obtained results are in line with available experimental data.Comment: 15 pages, 8 figures. Submitted to J. Phys: Condens. Matte

An ab initio study of the structural and mechanical alterations of Ti-Nb alloys

This article describes a systematic theoretical investigation of the role of Nb substitution on the structural and mechanical properties of Ti-Nb alloys. The aim is to understand the origin of the low-rigidity of some of these materials. This quality makes these materials suitable for metallic implants. The mechanical stability conditions in conjunction with the calculated elastic constants of Ti-Nb predict the destabilization of α′ and ω structures, while the β-phase can be stabilized for Nb content above 10 at. %. The evaluated Young's moduli (E) follow the sequence of Eω > Eα′ > Εα″ > Εβ, revealing high Eω and Eα′ values (greater than 120 GPa), while the Eβ value converges to approximately 87 GPa. The averaged E, estimated from a weighted average of Eω, Eα′, Εα″, and Εβ ab initio values, reproduces the experimental Ti-Nb Young's modulus w-shaped curve. Young's modulus surface reveals highly anisotropic E values for all Ti-Nb phases. Eβ exhibits values under 30 GPa along the [100] direction for Nb compositions larger than 12 at. %, suggesting that the orientational growth of a Ti-Nb alloy is important for the design of low-rigidity alloys, especially at small Nb concentrations. These results can be used as a guide for the design of novel low-rigidity alloys for biomedical applications

Local atomic order, electronic structure and electron transport properties of Cu-Zr metallic glasses

We studied atomic and electronic structures of binary Cu-Zr metallic glasses (MGs) using combined experimental and computational methods including X-ray absorption fine structure spectroscopy, electrical resistivity, thermoelectric power (TEP) measurements, molecular dynamics (MD) simulations, and ab-initio calculations. The results of MD simulations and extended X-ray absorption fine structure analysis indicate that atomic order of Cu-Zr MGs and can be described in terms of interpenetrating icosahedral-like clusters involving five-fold symmetry. MD configurations were used as an input for calculations of theoretical electronic density of states (DOS) functions which exhibits good agreement with the experimental X-ray absorption near-edge spectra. We found no indication of minimum of DOS at Fermi energy predicted by Mott's nearly free electron (NFE) model for glass-forming alloys. The theoretical DOS was subsequently used to test Mott's model describing the temperature variation of electrical resistivity and thermoelectric power of transition metal-based MGs. We demonstrate that the measured temperature variations of electrical resistivity and TEP remain in a contradiction with this model. On the other hand, the experimental temperature dependence of electrical resistivity can be explained by incipient localization of conduction electrons. It is shown that weak localization model works up to relatively high temperatures when localization is destroyed by phonons. Our results indicate that electron transport properties of Cu-Zr MGs are dominated by localization effects rather than by electronic structure. We suggest that NFE model fails to explain a relatively high glass-forming ability of binary Cu-Zr alloy

Self-diffusion of adatoms, dimers, and vacancies on Cu(100)

We use ab initio static relaxation methods and semi-empirical molecular-dynamics simulations to investigate the energetics and dynamics of the diffusion of adatoms, dimers, and vacancies on Cu(100). It is found that the dynamical energy barriers for diffusion are well approximated by the static, 0 K barriers and that prefactors do not depend sensitively on the species undergoing diffusion. The ab initio barriers are observed to be significantly lower when calculated within the generalized-gradient approximation (GGA) rather than in the local-density approximation (LDA). Our calculations predict that surface diffusion should proceed primarily via the diffusion of vacancies. Adatoms are found to migrate most easily via a jump mechanism. This is the case, also, of dimers, even though the corresponding barrier is slightly larger than it is for adatoms. We observe, further, that dimers diffuse more readily than they can dissociate. Our results are discussed in the context of recent submonolayer growth experiments of Cu(100).Comment: Submitted to the Physical Review B; 15 pages including postscript figures; see also http://www.centrcn.umontreal.ca/~lewi

Ab-initio and experimental study of phase stability of Ti-Nb alloys

A systematic theoretical and experimental study concerning the crystallographic structure and electronic properties of Ti-xNb (x 18.75 at% the β-phase is favoured against all other crystallographic structures in line with the experimental results. Interestingly, at high Nb content the α′ and ω hexagonal phases become unstable due to the electronic band filling close to the Fermi level EF, which is mainly due to Nb-p and Ti-d antibonding hybridizations. On the contrary, in the cubic β-Ti-25Nb (at%) the depletion of the occupied electronic states at the EF occurs mainly due to Nb-d and Ti-d bonding interactions, resulting in a stable β-TiNb structure. These data could enlighten the electronic origin of the Ti-Nb phase stability, thus, may contribute to the design of β stabilized low moduli Ti-based alloys suitable for load-bearing biomedical applications