[en] The development of new materials impacts on all branches of engineering and, in particular, in aerospace engineering. Metallic glasses (MG) are relatively newcomers to the world of materials science and have excellent mechanical properties; its study is mandatory to allow its technological implantation. The macroscopic mechanical properties of a material are linked to its atomic structure. In particular, the fracture behaviour of brittle materials is initiated by the generation of vibrational modes. In metallic glasses, with amorphous structure, the vibrational spectrum has specific features. In this work, the vibrational properties of metallic glasses are examined by Molecular Dynamics simulations. The study was focused in binary systems, which were simulated using different interatomic potentials: Lennard-Jones (LJ), Morse and the semiempirical Embedded atom (EAM). As in Pd-based metallic glasses, the ratio of masses of both species was high, namely 2 in Lennard-Jones potentials and 1.67 in Morse and EAM potentials. The large scale simulations allowed us to simulate systems with nm-scale heterogeneities frozen-in during the quenching process. Different relaxation states were obtained by changing the quenching rates of the simulated MGs. The collective vibrational atomic dynamics of metallic glasses is a longstanding subject of debate. The origin of the excess of vibrational modes known as Boson Peak (BP) is not clear. In the systems analysed the dependence of the BP position and intensity on the system size is found to be weak. On the contrary, the BP intensity increases with the quenching rate, while its position shifts slightly to smaller frequencies. The results obtained by using realistic, semiempirical EAM potentials compare well with the experimental data available in glasses of similar compositions. The dynamic structure factor, S(q, ¿) is also computed in large systems to get information on the behaviour of acoustic excitations at low wavenumbers. The dominant frequencies O L,T (q) are determined for each considered wavevector, in order to compute the relation of dispersion of longitudinal and transverse phonons. In all studied cases the width of the peak, GL,T(q), increases as the frequency increases. A linear region at low wavenumbers and a bending when approaching the limit of the first pseudo-Brillouin zone are found. This behaviour is the same than that observed experimentally by Inelastic X-Ray scattering. The macroscopic sound speed is obtained for wavenumbers tending to zero. The values obtained with EAM and Morse potentials are in qualitative agreement with those obtained experimentally in systems of similar composition. The Ioffer-Regel limit (IR), where the coherence length of the phonon is similar to the phonon wavelenght, was computed. It is found that the Ioffe-Regel frequency decreases slightly when applying faster quenching rates. The longitudinal Ioffe-Regel limit was found at frequencies higher than the Boson peak frequency for all the cases, although the diference in EAM systems is much reduced. Contrary to the typical results obtained in the LJ systems, using EAM potentials in both Cu20Pd80 and Cu50Pd50 the longitudinal IR limit is very close to the position of the BP while the transversal IR limit is found well below. This behaviour is coherent with that found by IXS measurements. It is inferred that the EAM potential increases the interaction between the longitudinal modes and the BP excess states. Finally the fragility of the studied systems was obtained by calculating the viscosity at different temperatures. Lennard-Jones systems showed a much larger fragility than EAM and Morse systems. However, even systems simulated with more realistic potentials showed fragility values much higher than those obtained experimentally. This is attributed to the extremely high quenching rates used in simulations, as it is known experimentally that fragility increases with the quenching rate.

[es] El desarrollo de nuevos materiales tiene impacto en todas las áreas de la ingeniería, y en particular en la ingeniería aeronáutica. Los vidrios metálicos son materiales relativamente nuevos, con excelentes propiedades mecánicas; su estudio es imprescindible para su implantación tecnológica. Las propiedades mecánicas macroscópicas de un material están determinadas por su estructura atómica. En particular, la fractura de materiales frágiles se inicia mediante la generación de modos de vibración. En los vidrios metálicos, con estructura amorfa, el espectro vibracional tiene características específicas. En este trabajo es estudian las propiedades vibracionales de los vidrios metálicos usando simulaciones de Dinámica Molecular. Se estudian sistemas binarios utilizando potenciales atómicos de Lennard-Jones (LJ), Morse y Embedded atom method (EAM), de tipo semiempírico. Al igual que en los vidrios metálicos a base de Pd, la relación de las masas de ambas especies en las simulaciones es alta, siendo 2 en potenciales LJ y 1,67 en potenciales Morse y EAM. Las simulaciones a gran escala permiten simular sistemas con heterogeneidades a escala nm, y la simulación a diferente velocidad de enfriamiento permite obtener configuraciones en diferentes estados de relajación. La dinámica de las vibraciones atómicas de los vidrios metálicos es un tema de debate. El origen del exceso de los modos de vibración conocidos como Boson Peak (BP) no es claro. En los sistemas analizados se observa que la dependencia de la posición y la intensitat del BP con el tamaño del sistema es débil. Por el contrario, la intensidad del BP aumenta con la velocidad de enfriamiento, mientras que su posición se desplaza a frecuencias menores. Los resultados obtenidos mediante el uso de potenciales realistas tipo EAM coinciden con los datos experimentales disponibles en vidres de composiciones similares. El factor de estructura dinámica, S (q, ω) se calculó también en sistemas grandes para obtener información sobre el comportamiento de las excitaciones acústicas en vectores de onda bajos. También se obtuvo la relación de dispersión de fonones longitudinales y transversales. En todos los casos estudiados, la anchura del pico, ΓL,T(q), aumenta con la frecuencia. Se observa una región lineal a bajos números de onda y una flexión cuando se aproxima al límite de la pseudo-zona de Brillouin. Este comportamiento es el mismo que el observado experimentalmente por Inelastic X-Ray scattering (IXS). La velocidad macroscópica del sonido se obtiene cuando el número de onda tiende a cero. Los valores obtenidos con potenciales EAM y Morse están en concordancia con los valores medidos experimentalmente en sistemas de composición similar. El límite Ioffer-Regel (IR) se define como la frecuencia a la que la longitud de coherencia del fonón es similar a la longitud de onda del fonón. Se observa que la frecuencia IR aumenta ligeramente en sistemas más relajados. El límite IR longitudinal se encuentra en todos los casos en frecuencias superiores a la frecuencia del BP. Contrariamente a los resultados obtenidos en los sistemas LJ, el límite longitudinal IR en Cu20Pd80 y Cu50Pd50 está muy cerca de la posición del BP, mientras que el límite transversal de IR se encuentra muy por debajo; estos resultados coinciden con los obtenidos experimentalmente mediante IXS. Se infiere que el potencial EAM aumenta la interacción entre los modos longitudinales y el exceso de estados de BP. Finalmente se obtuvo la fragilidad en los sistemas estudiados mediante el cálculo de la viscosidad a diferentes temperaturas. Es conocido experimentalmente que la fragilidad aumenta con la velocidad de enfriamiento. Los sistemas LJ mostraron una fragilidad mucho mayor que los sistemas EAM y Morse. Sin embargo, incluso los sistemas simulados con potenciales realistas mostraron valores de fragilidad muy superiores a los obtenidos experimentalmente debido a las extremas velocidades de enfriamiento de las simulaciones.