Microscopic origin of Boson peak in amorphous solids

This paper proposes a non-analytic model demonstrating that the Boson peak in amorphous solids originates solely from fluctuations in coordination numbers, while fluctuations in spring strength primarily contribute only to damping.

The Gist

Imagine you have a giant, invisible trampoline made of thousands of tiny springs and knots. In a perfect crystal (like a diamond), every knot is tied to exactly the same number of springs, and every spring is stretched to the exact same tension. If you pluck this perfect trampoline, it vibrates in a very predictable, orderly way.

Now, imagine an amorphous solid (like glass or plastic). It's still a network of springs and knots, but it's messy. The knots aren't in perfect rows, and the springs aren't all the same length. Scientists have been puzzled for decades by a strange "hiccup" in how these messy materials vibrate at low frequencies. They call this the Boson Peak. It's like an extra, unexpected drumbeat that shouldn't be there according to the standard rules of physics.

This paper by Cunyuan Jiang tries to solve the mystery of where this extra drumbeat comes from. The author breaks the problem down into two possible causes of "messiness":

1. The "Tension" Factor: Some springs are tighter or looser than others (fluctuation of spring strength).
2. The "Connection" Factor: Some knots are tied to 3 springs, some to 4, and some to 5 (fluctuation of coordination numbers).

The Experiment: Two Types of Mess

The author built a computer model of this spring network to test which factor causes the extra drumbeat.

• Scenario A (The Tension Test): Imagine a grid where every knot is tied to exactly 4 neighbors (like a perfect square grid). However, because the knots are slightly shifted from their perfect spots, the tightness of the springs varies.

  • Result: This only made the high-pitched vibrations sound a bit "muffled" or dampened. It did not create the extra low-frequency drumbeat (the Boson Peak).

• Scenario B (The Connection Test): Imagine a grid where the knots are still shifted, but now, the rules change: if two knots are close enough, they get a spring. If they are far, they don't. This means some knots end up with 3 springs, some with 4, and some with 5.

  • Result: Bingo. As soon as the number of connections varied, the extra low-frequency drumbeat (the Boson Peak) appeared.

The "Toy Model" Analogy

To explain why this happens, the author used a tiny model with just nine knots (like a 3x3 grid).

• The Perfect Grid: If every knot has exactly 4 springs, the system has two specific "notes" it can play.
• The Broken Grid: If you add one extra spring to a knot (making it have 5 connections) or remove one (making it have 3), the system suddenly gains two new notes that it couldn't play before.

These new notes are the Boson Peak. The paper shows that these new vibrations aren't just happening at the specific knot that changed; they ripple out and involve almost the entire network. It's like if one person in a choir changes their pitch slightly, and suddenly the whole choir starts humming a new, unexpected harmony.

The Big Conclusion

The paper argues that the Boson Peak isn't caused by springs being too tight or too loose. Instead, it is caused entirely by the uneven number of connections between particles.

• Spring strength (tension) just adds a little static or damping (like a blanket over a speaker).
• Coordination number (how many neighbors you have) is the sole architect of the extra vibrations.

In short: The "messiness" of amorphous solids isn't just about how far apart things are; it's about the fact that some things are holding hands with more neighbors than others. That specific type of social imbalance in the atomic network is what creates the mysterious Boson Peak.

Technical Summary

Technical Summary: Microscopic Origin of Boson Peak in Amorphous Solids

Problem Statement
The origin of the Boson peak (BP)—an anomalous excess in the vibrational density of states (DOS) of amorphous solids at low frequencies beyond the Debye prediction—remains a subject of significant debate. While various theories attribute the BP to anharmonicity, spatial fluctuations of elasticity, or phenomenological dynamical objects, the specific role of microscopic atomic-scale structural disorder has been less explored. Specifically, the dynamical matrix (Hessian) of an amorphous system contains two primary sources of disorder: the fluctuation of interaction strengths (spring constants) and the fluctuation of coordination numbers (the number of neighbors). The precise contribution of each factor to the emergence of the BP is not fully resolved.

Methodology
The authors propose a non-analytic model based on a scalar dynamical matrix representing a network of springs and nodes. To isolate the effects of different disorder types, the spatial degrees of freedom (longitudinal/transverse modes) are simplified, focusing on two distinct factors:

1. Spring Strength Fluctuation: The strength of the springs ($t_{ij}$) varies inversely with the distance between nodes ($t_{ij} = t_0/r_{ij}$).
2. Coordination Number Fluctuation: The number of springs connected to a node varies based on the spatial distribution of nodes.

The study constructs two specific scenarios on a two-dimensional square lattice with $N$ nodes, where nodes are displaced from their lattice sites by a random distance $\delta$:

• Scenario A (Fixed Coordination): Nodes are connected only to their nearest neighbors of the underlying square lattice. The coordination number is fixed at 4 for all nodes, but spring strengths fluctuate due to random displacements.
• Scenario B (Fluctuating Coordination): Nodes are connected to all neighbors within a cutoff distance ($a + \delta$). Both spring strengths and coordination numbers fluctuate.

The vibrational DOS is calculated by analyzing the eigenvalues of the dynamical matrix. The authors further validate their model by comparing the coordination number distribution of their network with simulation data from a realistic two-dimensional amorphous solid. Finally, a "toy model" consisting of a 9-node grid is used to analytically trace how local changes in coordination numbers (adding or removing a single spring) generate new dynamical modes.

Key Results

1. Origin of the Boson Peak: The simulations demonstrate that the Boson peak arises solely from the fluctuation of coordination numbers. In Scenario A, where coordination numbers are fixed at 4, the low-frequency DOS strictly follows Debye's law, regardless of the magnitude of spring strength fluctuations. The Van Hove singularity is merely weakened (damped) at high frequencies, but no low-frequency excess is generated.
2. Role of Spring Strength: Fluctuations in spring strength alone do not induce an anomalous DOS in the low-frequency region. Their primary effect is to introduce damping that weakens the Van Hove peak at higher frequencies.
3. Mechanism of Mode Generation: The toy model analysis reveals that when a node's coordination number deviates from the average (e.g., 3 or 5 instead of 4), it creates additional poles in the Green's function. These poles correspond to new dynamical modes at frequencies distinct from the crystalline spectrum. Crucially, these new modes are extended rather than localized; even nodes with the average coordination number contribute to the additional peaks, indicating that the BP involves a large fraction of particles in the system.
4. Robustness: The emergence of the BP in the low-frequency region ($\omega \in [0, 1]$) with a peak around $\omega = 0.5$ is robust against variations in the disorder level ($\delta$) and persists even when spring strength fluctuations are removed.

Significance and Claims
The paper claims to provide a direct, simplified answer to the puzzle of the microscopic origin of the Boson peak. By converting the complexity of amorphous dynamics into two fundamental factors, the work identifies the fluctuation of coordination numbers as the exclusive microscopic origin of the anomalous low-frequency DOS. The authors argue that this finding clarifies why the BP is a universal feature of amorphous solids: it is an intrinsic consequence of the topological disorder (variable connectivity) inherent in the network structure, rather than a result of elastic heterogeneity or anharmonicity. The study suggests that future analytical theories should focus on modeling the BP region using order parameters related to coordination number fluctuations.