Investigación mediante Dinámica Molecular del Corte Ultrapreciso de Cerio con Diamante

Resumen

El acoplamiento entre transformaciones de fase estructurales y la dinámica de dislocaciones complica la comprensión de la deformación nanométrica del cerio. En este trabajo se emplea dinámica molecular (DM) para desentrañar el mecanismo del cerio durante el corte ultrapreciso con diamante. Se construye un modelo de DM que incorpora potenciales Ce–Ce y Ce–C derivados empíricamente, y se evalúa la respuesta elástica de dos fases fcc de cerio. Las simulaciones muestran que el deslizamiento por dislocaciones domina la deformación plástica, mientras que el análisis de la función de distribución radial revela únicamente cambios menores de fase γ→δ en la superficie mecanizada y en la viruta. Los estudios paramétricos identifican las condiciones óptimas de mecanizado para lograr una superficie de alta calidad.

Antecedentes

El cerio (Ce, Z=58) es un metal lantánido valorado por sus propiedades mecánicas, físicas y químicas únicas. La morfología superficial influye críticamente en el rendimiento, la resistencia a la corrosión y el comportamiento de almacenamiento de hidrógeno de las piezas de cerio [1–3]. El corte ultrapreciso con diamante puede producir superficies ultrasuaves con daño subsuperficial mínimo [6,7], sin embargo, existe poca investigación experimental o teórica sobre el corte de cerio. Dado que el radio del filo de la herramienta en el corte ultrapreciso es comparable a la profundidad de corte, las propiedades del material dominan el proceso, lo que hace que el mecanismo de corte del cerio sea particularmente desafiante.

El Ce exhibe un diagrama de fases presión-temperatura rico, impulsado por la deslocalización de electrones 4f. A presión ambiente y <110 K, la fase α-Ce (fcc) es estable, transformándose a β-Ce (dhcp) entre 45–275 K, γ-Ce (fcc) entre 270–999 K, y δ-Ce (bcc) entre 999–1071 K [8–11]. Notablemente, la transformación γ→α a 295 K y 8 kbar colapsa el volumen en ~20% [12–14], mientras que la transición γ→δ puede desencadenarse por alta temperatura y presión. Estas transformaciones alteran el enlace y las propiedades elásticas, afectando la deformación. El Ce metálico también es dúctil, gobernado por la actividad de dislocaciones [15]. Sin embargo, la interacción entre cambios de fase y dislocaciones durante el corte con diamante permanece prácticamente inexplorada.

La simulación por DM ofrece una ruta poderosa para sondear mecanismos fundamentales de mecanizado a escala atómica. Estudios previos han examinado el rayado nanométrico y el corte ortogonal de Cu, Al, Fe y aleaciones de alta entropía, revelando conocimientos sobre nucleación de dislocaciones, acumulación superficial y comportamiento de fricción [16–27]. Dichos trabajos destacan la importancia de la orientación cristalográfica, el ángulo de ataque de la herramienta y los potenciales de interacción. Sin embargo, ningún estudio de DM ha investigado el mecanizado mecánico del cerio, lo que motiva el presente trabajo.

Métodos

Modelo de DM del Corte con Diamante

La Figura 1 ilustra la configuración de la simulación: una pieza de trabajo de Ce monocristalino (41 × 25 × 31 nm) que contiene ~1 × 10⁶ átomos de γ-Ce, y una herramienta de corte de diamante (0.1 × 10⁶ átomos de C) con un ángulo de despeje de 9°. La pieza de trabajo se divide en una capa fija inferior de 2 nm y una capa móvil que sigue la dinámica newtoniana con un paso temporal de 1 fs. Se aplican fronteras periódicas solo en la dirección longitudinal. Se consideran tres orientaciones cristalográficas: (010), (110) y (111). Se probaron siete ángulos de ataque (−30°, −20°, −10°, 0°, 10°, 20°, 30°). La herramienta de diamante se trata como cuerpo rígido; por tanto, se omiten las interacciones C–C.

Las interacciones atómicas se modelan de la siguiente manera: las fuerzas Ce–Ce usan el método del átomo empotrado (EAM) parametrizado por Sheng et al. [29], reproduciendo con precisión las constantes elásticas del Ce fcc. Las interacciones Ce–C adoptan un potencial Morse:

$$E_{\text{tot}}=\sum_{ij}D_0\left[e^{-2\alpha(r-r_0)}-2e^{-\alpha(r-r_0)}\right]$$

con D₀=0.087 eV, α=5.14, y r₀=2.93 Å; se aplica un corte de 1.0 nm [30].

El sistema se equilibra primero a 30 K, 0 bar en el ensamble NPT, y luego se somete a corte con diamante a velocidad constante de 100 m/s y profundidad de corte de 4 nm en el ensamble NVT. Las fuerzas de corte se resuelven en la dirección de corte, indicada por flechas rojas en las figuras. Aunque las dimensiones simuladas y las velocidades de corte son órdenes de magnitud menores que las condiciones experimentales ultraprecisas, el modelo de DM captura la física esencial del proceso. El análisis de vecinos comunes (CNA) distingue átomos fcc, hcp, bcc y defectos [31]. Las simulaciones se realizan con LAMMPS (paso de 1 fs) y se visualizan con OVITO [32,33].

Caracterización de las Fases del Cerio

Se modelan cinco fases de cerio —γ, α, β, ε y δ— usando parámetros de red de la literatura [8–11]. Configuraciones a granel se someten a tracción uniaxial, cizallamiento y compresión para extraer constantes elásticas y propiedades mecánicas (Tabla 2). El módulo de Young calculado del γ-Ce (24.17 GPa) se alinea con datos experimentales de nanoindentación (36.7 GPa) [10], y las constantes elásticas concuerdan con resultados de dispersión inelástica de neutrones [34]. El α-Ce, más denso, exhibe una rigidez significativamente mayor que el γ-Ce.

El análisis de la función de distribución radial (FDR) caracteriza diferencias estructurales entre fases. La Figura 2a muestra picos FDR distintos para cada fase, correspondientes a distancias de vecinos más cercanos (ej., γ-Ce: 3.64, 5.13, 6.30 Å; α-Ce: 3.41, 4.85, 5.92 Å). La compresión uniforme del γ-Ce para lograr un colapso de volumen del 20% produce una FDR que coincide con la del α-Ce (Figura 2b), confirmando la transición γ→α [12–14].

Resultados y Discusión

Mecanismos de Mecanizado del Cerio

Las simulaciones de DM del corte con diamante sobre Ce(010) con ángulo de ataque 0° revelan tres regímenes de fuerza (Figura 3). Inicialmente, fuerzas adhesivas producen fuerzas de corte y normales negativas; sigue la deformación elástica, llevando a un rápido aumento de las fuerzas. A ~2.3 nm de longitud de corte, una caída brusca indica el inicio de la plasticidad, dominada por dislocaciones parciales 1/6<112> de Shockley nucleadas cerca de la superficie libre y deslizándose sobre planos {111}. Las fuerzas de corte y normales se estabilizan entre 10–35 nm, mientras las dislocaciones se acumulan delante y debajo de la herramienta. Más allá de 35 nm, la herramienta se separa de la pieza, causando reducción de fuerzas. En todo momento, la fuerza normal permanece menor que la de corte.

La Figura 4 muestra la evolución de estructuras de defectos y morfología superficial en longitudes de corte clave. Las parciales de Shockley generan fallas de apilamiento limitadas por núcleos de dislocación. A medida que avanza el corte, una densa nube de parciales emana de la superficie frontal, formando virutas sobre la cara de ataque (Figura 4f). Las dislocaciones detrás de la herramienta ascienden a la superficie superior, promoviendo la recuperación superficial. Cuando la herramienta se acerca al límite izquierdo, la propagación de dislocaciones se ve obstaculizada, aumentando el volumen de viruta (Figura 4g). Tras la separación de la viruta, la densidad de dislocaciones disminuye por aniquilación en la superficie libre (Figura 4h).

El análisis FDR de zonas de defectos confirma una transformación γ→δ limitada. La FDR de la región de defecto bajo la superficie se alinea con picos de δ-Ce, indicando una transición γ→δ menor impulsada por calentamiento local (Figura 5a). Firmas similares de δ-Ce aparecen en la viruta (Figura 5b), pero la cantidad total permanece trivial, sugiriendo que los cambios de fase no son dominantes durante el corte.

Influencia del Ángulo de Ataque

La Figura 6 muestra que tanto la fuerza de corte como la normal disminuyen al aumentar el ángulo de ataque, de acuerdo con la teoría del plano de cizallamiento de Merchant [37]. Ángulos de ataque mayores también reducen la diferencia de fuerzas y mejoran la estabilidad del mecanizado.

El análisis de defectos (Figura 7a,b) revela que un ángulo de ataque negativo (−30°) produce mayor densidad de dislocaciones que uno positivo (30°). El número de segmentos de dislocación (Figura 7c) disminuye monótonamente con el ángulo de ataque, reflejando deformación plástica menos compleja y menor endurecimiento por trabajo.

La morfología superficial (Figura 8) corrobora estos hallazgos: ángulos de ataque negativos producen acumulación superficial asimétrica y pronunciada, mientras que 30° genera el menor volumen de acumulación y un perfil simétrico. Por tanto, un ángulo de ataque de 30° es óptimo para lograr baja fuerza de corte, densidad mínima de dislocaciones y acabado superficial superior.

Influencia de la Orientación Cristalográfica

Bajo el ángulo de ataque óptimo de 30°, se comparan las orientaciones Ce(010), Ce(110) y Ce(111) (Figura 9). Ce(010) exhibe la menor fuerza de corte y conteo de dislocaciones. Aunque Ce(111) ofrece las menores dislocaciones, su plano {111} paralelo a la superficie promueve acumulación superficial, aumentando la resistencia al corte. La morfología de la superficie mecanizada (Figura 10) muestra que Ce(010) produce el menor volumen de acumulación y el perfil más simétrico, mientras que Ce(110) sufre la mayor acumulación.

Conclusiones

Hemos desarrollado un marco de DM que incorpora potenciales EAM y Morse para modelar el corte ultrapreciso de cerio con diamante. Las constantes elásticas, propiedades mecánicas y comportamiento de transición de fase de las fases de Ce se validan, confirmando la capacidad predictiva del modelo. Las simulaciones revelan que la nucleación y deslizamiento de dislocaciones dominan la plasticidad, observándose solo cambios menores γ→δ. Las condiciones óptimas de mecanizado —ángulo de ataque 30° y orientación cristalográfica (010)— producen las menores fuerzas de corte, densidad mínima de dislocaciones y mayor calidad superficial.

Referencias

[1]–[38] (consulte el artículo original para la lista completa de citas)

Figuras

Figura 1. Modelo de DM del corte de Ce con diamante. (a) Vista frontal, (b) vista superior. Rojo = átomos de Ce inferiores, azul = átomos de Ce móviles, gris = átomos de C.

Figura 2. (a) FDR de fases de cerio. (b) FDR antes y después de la compresión del γ-Ce al 20% de colapso volumétrico.

Figura 3. Variación de fuerzas de corte y normales con la longitud de corte. Subfiguras muestran configuraciones de defectos por zona.

Figura 4. Estructuras de defectos y morfología de superficie mecanizada en varias longitudes de corte.

Figura 5. Análisis FDR de zona de defecto (a) y viruta (b) tras el corte.

Figura 6. Dependencia del ángulo de ataque de fuerzas promedio de corte y normales.

Figura 7. Estructuras de defectos para ángulos de ataque −30° (a) y 30° (b); conteo de dislocaciones (c).

Figura 8. Morfologías de superficie mecanizada para ángulos de ataque −10°, −20°, −30°, 10°, 20°, 30°.

Figura 9. Fuerza de corte y número de dislocaciones versus orientación cristalográfica.

Figura 10. Morfologías de superficie mecanizada para Ce(010), Ce(110), Ce(111).