Thermal Nanoquakes: Terahertz Frequency Surface Rayleigh Waves in Diamond Nanocrystals

Mechanical THz vibrations in nanocrystals have recently been harnessed for quantum sensing and thermal management, and this study elucidates these features by using atomistic molecular dynamics simulations and experimental neutron spectroscopy to reveal thermally induced Rayleigh surface phonons as the cause of linear scaling in the low-energy vibrational density of states.

The gist: Thermally induced Rayleigh surface phonons, which have a low group velocity and an amplitude that decays exponentially away from the surface, are responsible for the previously observed but unexplained linear scaling of the low-energy vibrational density of states in nanocrystals.

Unique Nanodiamond Lattice Vibrations from Molecular Dynamics

Molecular dynamics (MD) simulations reveal characteristic surface features in nanoscale diamond distinct from infinite bulk diamond, including several unique energy regions:

1. Region i) Acoustic gap: This is determined by the length scale 2π/d, representing a forbidden region of k-space for wavelengths exceeding the particle diameter.

2. Region ii) Spheroidal Lamb mode resonances: Resonant modes appear in spherical NDs, with the strongest coinciding with band minima above the acoustic gap, matching the eigenfrequency of the fundamental Lamb mode (radial particle breathing). These modes are sensitive to structural disorder.

3. Region iii) Rayleigh surface acoustic wave: This is a large energy range (0–40 meV) where a linear scaling of the DOS, g(E) E, appears, distinct from the quadratic scaling found in bulk diamond. The mode is described as a primarily transverse SAW analogous to the slow surface component of an earthquake.

4. Region iv) Boson-like peak: A flat mode at 50 meV causes a bump in the DOS, with a view that it may be related to the resonant coupling between Rayleigh SAWs and a high-order Lamb wave or disorder-induced surface mode.

5. Region v) Modified bulk acoustic modes: Broadening and softening of bulk diamond modes are attributed to phonon scattering induced by disorder and confinement, leading to reduced phonon lifetimes, with a scattering rate d/c resulting in linewidth broadening 1/d at the meV scale.

6. Region vi) Modified bulk optical modes: Similar modifications occur due to significant lifetime broadening and mode softening observed via Raman spectroscopy.

7. Region vii) High-energy surface rattling modes: These arise from low-coordination (sp2 or sp) surface carbons, with the DOS increased and broadened in disordered core-shell NDs.

Isolating the Contribution of the Surface Acoustic Waves

To explicitly demonstrate that linearity arises from Rayleigh SAWs, the authors dissected the vibrational dynamics by separating core and surface contributions. The core atoms strongly resemble infinite diamond with quadratic Debye scaling of low-energy DOS, whereas outer shells clearly transition to a linear energy dependence. The corresponding dynamic structure factor maps confirm that the additional linearly dispersive mode is localized on the surface and provides the extra contribution to the DOS. Furthermore, imposing rigid (Dirichlet) boundary conditions causes this excess low-energy DOS associated with SAWs to vanish, returning the quadratic scaling of the DOS.

Experimental Confirmation of the Linear DOS from Rayleigh Waves

Inelastic neutron spectroscopy (INS) experiments were used to confirm these findings on nanodiamond powder and microdiamond powders. The INS-derived DOS for NDs showed an unambiguous linear scaling between 0-10 meV, with the linear component dominating up to 40 meV, consistent with the prediction of Rayleigh SAWs. This linearity was highlighted when the DOS was scaled by 1/E. Clear differences were observed between the microdiamond (quadratic dependence) and ND samples in this low-energy region, suggesting that the additional signal in the ND is primarily from the surface modes that give rise to the linear DOS.

Theoretical Formalism for 2D-DOS

The paper proposes a heuristic argument linking the linear DOS to the dimensionality of the NC surface. Since Debye DOS scales as g(E) E(n-1), where n is dimensionality, a linear DOS likely originates from the "2D" NC surface. The authors derive this by considering an elastic medium with a surface extending along x and y directions, leading to real wave-vectors in those directions and imaginary wave-vectors perpendicular to the surface. The resulting DOS is given by g(E) E(cR) squared, where cR is the Rayleigh wave group velocity. This linear, 2D-DOS is proposed as a unique feature of NCs hosting specific phononic surface states with radial exponential decay.

Conclusions and Implications

The study concludes that the surfaces of NDs undergo persistent violent thermal nanoquakes, which are ultra-fast and ultra-large relative to object diameter. Rayleigh phonons, demonstrated as the origin of the linearly scaling DOS, are expected to significantly alter thermal conductivity through altered phonon scattering and surface localization.

As a fastidious and diligent researcher, I have analyzed the provided scientific paper, Thermal Nanoquakes: Terahertz Frequency Surface Rayleigh Waves in Diamond Nanocrystals. The core findings revolve around identifying and characterizing novel vibrational modes (Rayleigh surface waves) in diamond nanocrystals (NDs) using atomistic molecular dynamics (MD) simulations and experimental neutron spectroscopy.

Here are the specific improvements that can be made to AI systems, derived from the scientific insights of this paper:

The scientific findings strongly suggest that AI systems can be significantly improved by integrating physics-informed modeling and high-fidelity surface dynamics into their architectures.

1. The discovery of a distinct, linearly scaling vibrational density of states (DOS) in nanocrystals is attributed to thermally induced Rayleigh surface phonons whose amplitude decays exponentially away from the surface.

2. This linear scaling is proposed to originate from an effective 2D DOS arising from confined, THz-scale Rayleigh SAWs propagating along the free surface of the nanocrystal.

3. The behavior of these modes is characterized by a low group velocity and exponential amplitude decay, distinct from bulk lattice phonons (Debye scaling).

Here are the specific improvements to AI systems:

1. A specialized physics-informed neural network (PINN) module should be developed for materials science simulations, specifically trained on classical molecular dynamics trajectories (like those used in this paper) to predict vibrational properties of nanoscale structures.

2. This PINN should be capable of distinguishing between bulk lattice modes and surface-localized, Rayleigh wave contributions based on the energy-dependent DOS scaling.

Here is what the improved AI system can do:

The improved AI system, leveraging the physics derived from this research, can perform the following specific tasks:

1. Analyse experimental data (e.g., Inelastic Neutron Scattering (INS) or Raman spectroscopy) obtained on nanomaterials to automatically identify and isolate features that correspond to surface modes versus bulk modes.

2. Predict the presence and nature of surface-localized phonon states in novel nanostructures by analyzing their predicted vibrational DOS scaling, specifically identifying linear scaling indicative of 2D Rayleigh wave excitations.

3. Design or optimize nanoscale materials (e.g., diamond coatings or composites) to enhance or suppress specific surface vibrational modes (like the Rayleigh SAW), thereby tailoring thermal conductivity, sensing capabilities, and phonon scattering properties for applications in phononic devices and thermometry.

4. Model the non-equilibrium dynamics of surface modes under different environmental conditions (e.g., temperature changes, chemical environments) to predict how these nanoacoustic features will influence energy transfer and catalytic reactions at the nanoscale.

Mechanical THz vibrations in nanocrystals have recently been harnessed for quantum sensing and thermal management. The free boundaries of nanocrystals introduce new surface wave solutions, analogous to the seismic waves on Earth, yet the implications of these surface waves on nanocrystals have remained largely unexplored. Here, we use atomistic molecular dynamics simulations and experimental neutron spectroscopy to elucidate these THz-scale features in nanodiamond. Our key insight is that thermally induced Rayleigh surface phonons, which have a low group velocity and an amplitude that decays exponentially away from the surface, are responsible for the previously observed but unexplained linear scaling of the low-energy vibrational density of states in nanocrystals. Large thermal atomic displacements, relative to the nanoparticle radius, induce perpetual surface quakes, even at ambient conditions. Normalised to the radius, the surface displacement ratio in diamond nanocrystals exceeds that of the largest recorded earthquakes by a factor of 10 5 plus or minus1. We explicate how these dramatic Rayleigh waves coexist with other distinctive features including confined lattice phonons, soft surface modes, the acoustic gap, Love waves, and Lamb modes, thereby offering a complete framework for the vibrational dynamics of nanocrystals.