Verrouillage de vitesse de particules sur graphène sous ondes de surface propagatives

Résumé

Une diffusion rapide induite par fluctuation thermique et vibration a été détectée à l'échelle nanométrique. Dans cet article, le mouvement d'une particule sur une couche de graphène avec onde de surface propagative est étudié par simulation de dynamique moléculaire (DM) et modèle théorique. Il est prouvé que la particule maintient sa vitesse à la vitesse de l'onde sous certaines conditions préalables, à savoir l'effet de verrouillage de vitesse. En exprimant le potentiel de van der Waals (vdW) entre la particule et la surface ondulée comme fonction des courbures, le mécanisme est clarifié basé sur la flaque de potentiel dans un repère d'onde relatif. Deux conditions préalables sont proposées : la position initiale de la particule doit se situer dans la flaque de potentiel, et l'énergie cinétique initiale ne doit pas permettre à la particule de sauter hors de la flaque. L'analyse paramétrique indique que la région de verrouillage de vitesse est affectée par la longueur d'onde, l'amplitude et le potentiel de paire particule-onde. Avec une longueur d'onde plus petite, une amplitude plus grande et un potentiel vdW plus fort, la région de verrouillage est plus large. Ce travail révèle un nouveau type de mouvement cohérent pour les particules sur matériau lamellaire basé sur la théorie de la flaque de potentiel, pouvant expliquer les phénomènes de diffusion rapide à l'échelle nanométrique.

Introduction

Récemment, une série de phénomènes de transport et diffusion rapides induits par ondes de surface/phonons ont été détectés à l'échelle micro/nanométrique. D'abord, les phénomènes thermophorétiques le long d'un nanotube de carbone ou d'un ruban de graphène ont été largement étudiés. Les fluctuations thermiques permettent un écoulement d'eau continu à travers un NTC en imposant un gradient thermique axial. Des simulations DM hors équilibre ont exploré la faisabilité d'utiliser un gradient thermique sur un grand substrat de graphène pour contrôler le mouvement d'un nanoflocon de graphène. Le transport de nanogouttelettes d'eau sur graphène et nitrure de bore hexagonal (h-BN) par gradient thermique a aussi été étudié par DM. Ces phénomènes sont corrélés à certains modes de phonons. Schoen et al. ont attribué le mouvement thermophorétique dans un NTC au mode de respiration du tube. Panizon et al. ont souligné que les ondes flexurales propagatives/phonons sur graphène peuvent transmettre leur momentum aux adsorbats et causer le transport. Similairement, Angelos et al. ont montré que les ondulations propagatives induites par la température sur graphène entraînent une diffusion rapide de nanogouttelettes d'eau, 2-3 ordres de grandeur plus rapide que l'auto-diffusion de l'eau liquide.

En plus des fluctuations thermiques, des études confirment que la vibration peut aussi transporter particules et gouttelettes dans et hors d'un NTC. Par exemple, des nanogouttelettes sont transportées le long du nanotube à une vitesse ~30 nm/ns quand des ondes acoustiques transversales polarisées linéairement transmettent du momentum linéaire. Guo et al. ont démontré que les molécules d'eau dans un cantilever vibrant sont entraînées par forces centrifuges et subissent un écoulement continu de l'extrémité fixe vers l'extrémité libre. Un transport unidirectionnel nanoscopique d'eau à travers un NTC simple paroi a été conçu utilisant une charge vibrante et un NTC composite à énergie de surface asymétrique. Zhou et al. ont étudié les inversions de courant dans une pompe à eau nanoscalée basée sur un NTC alimentée par vibration mécanique. Chang et Guo ont découvert l'onde domino dans les NTC qui peut projeter la molécule interne à grande vitesse (jusqu'à 1 km/s). Un processus domino réversible est aussi prouvé dans les NTC simple paroi.

Comme divers phénomènes de diffusion et transport rapides induits par fluctuations thermiques et vibration sont détectés, il est confirmé que le mouvement haut-bas de la surface améliore la diffusion et le transport. La connexion entre l'onde et le mouvement des particules reste floue et ne peut être unifiée. Une explication principale est que le momentum de la surface est transporté vers la particule ou la gouttelette hors de la surface. Mais la relation entre amplitude, fréquence et interaction particule-surface ne peut être déduite de cette explication. De plus, Angelos et al. ont souligné qu'une préférence nette pour un signe de courbure du graphène est nécessaire pour la diffusion rapide d'adsorbats, indiquant que le potentiel d'interaction induit par la morphologie de surface ondulée est étroitement lié à la diffusion rapide. Explorer l'interaction entre surface ondulée et particule externe est donc essentiel pour comprendre le mécanisme de transport et diffusion rapides à l'échelle nanométrique.

Dans cet article, en étudiant une particule hors surface de graphène ondulée basée sur l'interaction vdW décrite par le potentiel de paire Lennard-Jones (L-J), une relation cohérente entre le mouvement ondulatoire et la vitesse de la particule est démontrée par simulations DM. La vitesse globale d'une particule déposée sur la surface ondulée est confirmée comme restant égale à la vitesse de l'onde propagative sous certaines conditions, à savoir l'effet de verrouillage de vitesse. Ensuite, une théorie de la flaque de potentiel est construite basée sur le potentiel d'interaction entre particule et surface ondulée exprimé comme fonction des courbures. Avec cette théorie, deux conditions préalables pour le verrouillage de vitesse sont proposées, et la trajectoire et vitesse prédites concordent bien avec les résultats DM. L'effet de la longueur d'onde, l'amplitude et les paramètres d'interaction vdW est analysé, en accord avec la régulation détectée pour les phénomènes de surf de gouttelettes sur graphène. Le mécanisme de verrouillage de vitesse piloté par l'onde révèle une nouvelle relation cohérente entre vitesse de particule et surface ondulée.

Méthodes

La simulation DM est implémentée dans LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator). La surface ondulée est supposée être une couche de graphène avec densité atomique ρ = 3,85×10¹⁹ m⁻². La feuille de graphène est initialement plate (z = 0 Å), longueur 6344 Å en direction x (6000 atomes par cellule unitaire). Condition aux limites périodiques en y (période 12,2 Å). Une particule sphérique de masse m = 0,83×10⁻²⁵ kg est considérée pour simplifier le modèle et se concentrer sur l'effet géométrique. Initialement, la particule est à z = 7 Å, x = 200 Å, avec vitesse initiale -50 m/s en z et ~2000 m/s en x. En décalant le temps de départ de la vitesse en z, la position initiale sur la surface ondulée est contrôlée.

Le potentiel REBO (reactive empirical bond order) modélise les atomes de graphène. Le potentiel Lennard-Jones modélise l'interaction entre particule P et chaque atome de carbone : u(R) = ε[(σ/R)¹² - (σ/R)⁶], avec ε = 5,92×10⁻²¹ J, σ = 4×10⁻¹⁰ m. La hauteur d'équilibre entre P et la surface courbe est h = 4,2×10⁻¹⁰ m (force normale nulle, validé par simulation).

La fonction d'onde propagative est sinusoïdale : y = A sin(2πx/λ - ωt), avec A = 1×10⁻⁹ m, λ = 21,75×10⁻⁹ m. Fréquence angulaire ω = 2π/10⁻¹² (période 10 ps), vitesse d'onde v_onde = 2175 m/s = λω/2π. Pour déclencher l'onde, les 10 Å gauche du graphène (x ∈ [-10, 0] Å) sont oscillés en z avec l'amplitude et fréquence ci-dessus. Les atomes avec x > 6010 Å sont fixés pour stabiliser la feuille. Pour simuler du graphène plat, les atomes non fixés sont attachés à leurs positions initiales en z par un ressort faible (0,0938 eV/Ų, A=0).

Température initiale 5 K pour les atomes de carbone libres, pour éliminer les ondulations thermiquement activées et se concentrer sur l'effet de l'onde propagative par excitation mécanique. Évolution en ensemble micro-canonique (NVE) avec pas de temps 1 fs. Température quasi inchangée pendant toute la simulation.

Résultats et Discussion

La trajectoire des particules sur graphène ondulé et plat est illustrée (Fig. 1). L'intervalle temporel est la période de l'onde. La position relative de la particule ne change pas par rapport aux crêtes ou creux, signifiant que la particule est verrouillée sur la surface avec vitesse égale à la vitesse d'onde. En comparaison, le mouvement sur surface plate est plus lent ; la vitesse diminue rapidement par frottement, alors que le frottement semble inopérant sur la surface ondulée. D'autres cas DM avec différentes températures et paramètres d'onde sont en matériel supplémentaire. Des simulations avec atome Xe et molécule C₆₀ confirment la généralisabilité.

Pour comprendre le mécanisme, un modèle est construit considérant l'interaction entre surface S et particule P (Fig. 2a, b). Longueur d'onde λ, amplitude A, hauteur minimale h, densité surfacique ρ_s. L'interaction vdW est décrite par le potentiel L-J. L'interaction P-S s'écrit comme fonction de la courbure moyenne H et courbure de Gauss K basée sur le potentiel de paire L-J [Eq. 3], validée par comparaison avec le potentiel de surface DM (Fig. 2c).

Le frottement est estimé en simulant le mouvement sur graphène plat (matériel supplémentaire). Approximation par surface plate, frottement estimé f = -5,2×10⁻¹³ N. Le potentiel relatif dans le repère d'onde (PXY, Fig. 2b, rouge) mobile à la vitesse d'onde est construit. L'onde est "gelée" dans PXY. La particule se déplaçant vers la droite subit un frottement constant vers la gauche. Le potentiel relatif est P = U_n + f·x (Eq. 4), tracé Fig. 2d.

Puisque PXY se déplace avec l'onde, la position initiale dans le potentiel détermine la trajectoire. Deux conditions préalables basées sur Fig. 2d : (1) position initiale dans la flaque de potentiel (zone rouge λ₁) ; (2) énergie cinétique relative ½m(v₀ - v_onde)² ≤ ΔU. La particule reste piégée, oscille dans la flaque, et en coordonnées absolues, oscille mais avance à la vitesse verrouillée autour de v_onde (verrouillage de vitesse). Sinon (zone bleue λ₂ ou énergie > ΔU), la particule saute vers des flaques inférieures, décroche de l'onde jusqu'à un nouvel équilibre (ex: arrêt complet, frottement nul). Panizon et al. montrent que la différence de vitesse disperse l'onde et fournit une force propulsive, suggérant une vitesse finale > 0.

L'équation du mouvement est établie par les lois de Newton. Forces motrices : force normale F_n = ∂U/∂h et force tangentielle F_t = (∂U/∂H)∇H + (∂U/∂K)∇K (Eq. 5). Pour L-J, forces attractives et répulsives existent ; P reste à hauteur h où F_n = 0. Équation du mouvement en x : mẍ = F_x - f (Eq. 6), F_x composante de F_t en x. Pour onde sinusoïdale, K=0, H=κ (courbure en Ozx). Substitution donne la trajectoire résolue numériquement.

Exemples de verrouillage et déverrouillage (Fig. 3). Position initiale dans λ₁ (Fig. 3a) : trajectoires théorie vs DM comparées (Fig. 3c). La particule s'arrête vite sur graphène plat, mais continue sur surface ondulée. Trajectoire théorique approxime DM. Vitesse (Fig. 3c) oscille autour de v_onde. Position initiale dans λ₂ (Fig. 3d) : trajectoire tend vers constante (Fig. 3e), vitesse diminue vers zéro (Fig. 3f). L'onde améliore le mouvement même en zone déverrouillée vs surface plate.

Analyse paramétrique basée sur la théorie de la flaque (Fig. 4). Paramètres : λ, A, ω, paramètres L-J (friction supposée constante). Fig. 4a : profondeur de flaque diminue avec λ, pas de verrouillage au-delà d'une λ critique. Fréquence ω plus basse → λ plus grand → région de verrouillage diminue avec ω. Fig. 4b : profondeur augmente avec A, verrouillage disparaît si A trop petit. Ratio A/λ ne doit pas être trop grand pour éviter dommages. Fig. 4c : profondeur diminue rapidement si λ et A augmentent en gardant λ/A fixe. Force motrice basée sur courbure diminue vite avec l'échelle, verrouillage disparaît pour ondes grande échelle. Fig. 4d : région de verrouillage plus large si potentiel de paire ε fort, disparaît si ε faible.

La rigidité et paramètres L-J diffèrent pour autres nanomatériaux 2D, menant à différentes fréquences et vitesses d'onde. Selon l'analyse, la flaque apparaît en choisissant λ et A appropriés. Comme la flaque est la condition d'accouplement particule-surface, le verrouillage s'établira pour de nombreux nanomatériaux 2D sous interaction courte portée.

Bien qu'une seule particule soit discutée, le cadre reste thermodynamique. La flaque est la condition essentielle pour le mouvement couplé. Pour multiples particules, si toutes dans la flaque et conditions satisfaites, elles seront piégées et se déplaceront avec l'onde. Le mouvement est contrôlable en ajustant λ et A. Région de verrouillage plus grande pour λ petit, A grand, ω élevé → diffusion rapide améliorée. L'analyse paramétrique concorde avec la régulation de diffusion rapide dans d'autres travaux. Angelos et al. : coefficient de diffusion augmente avec amplitude des ondulations graphène [22]. Amplitude accrue → préférence accrue pour les vallées, expliqué par Fig. 4b. Amplitude suffisante → région de verrouillage couvrant toute la longueur d'onde → diffusion améliorée. Potentiel vallée < potentiel crête [22] (Fig. 4), correspondant au potentiel plus bas pour région crête dans Fig. 4. Cao et al. : vitesse d'écoulement dans nano-canal augmente avec amplitude et fréquence [45], conforme à l'analyse.

La DM ne reflète que de très courts temps. Applications potentielles déductibles du mécanisme : en ajustant amplitude et fréquence, réalisation de régions quasi verrouillage ou déverrouillage pour faire bouger ou arrêter les particules. Mouvement vertical de surface transformé en mouvement transversal de particule, type mouvement cliquet, utilisable en systèmes nano-électro-mécaniques. Trajectoire différente selon potentiel de paire → séparation de phase possible.

Conclusions

Nous démontrons une relation distinctive entre particule et couche de graphène avec onde de surface propagative : phénomène de verrouillage de vitesse. Par DM, confirmé que la vitesse de la particule reste autour de la vitesse d'onde sous certaines conditions. Modèle théorique élucidant le mécanisme, où la flaque de potentiel domine l'effet de verrouillage. Conditions de verrouillage proposées : position initiale dans flaque, énergie cinétique initiale insuffisante pour en sortir. Trajectoire théorique en accord avec DM. Effet de λ, A et paramètre L-J discuté. Travail apportant nouvelle perspective pour comprendre diffusion et transport rapides sur surface ondulée et applications potentielles pour séparations de phase.