Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userThe Stochastic Self-Consistent Harmonic Approximation: Calculating Vibrational Properties of Materials with Full Quantum and Anharmonic Effects
94
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
The efficient and accurate calculation of how ionic quantum and thermal fluctuations impact the free energy of a crystal, its atomic structure, and phonon spectrum is one of the main challenges of solid state physics, especially when strong anharmonicy invalidates any perturbative approach. To tackle this problem, we present the implementation on a modular Python code of the stochastic self-consistent harmonic approximation method. This technique rigorously describes the full thermodyamics of crystals accounting for nuclear quantum and thermal anharmonic fluctuations. The approach requires the evaluation of the Born-Oppenheimer energy, as well as its derivatives with respect to ionic positions (forces) and cell parameters (stress tensor) in supercells, which can be provided, for instance, by first principles density-functional-theory codes. The method performs crystal geometry relaxation on the quantum free energy landscape, optimizing the free energy with respect to all degrees of freedom of the crystal structure. It can be used to determine the phase diagram of any crystal at finite temperature. It enables the calculation of phase boundaries for both first-order and second-order phase transitions from the Hessian of the free energy. Finally, the code can also compute the anharmonic phonon spectra, including the phonon linewidths, as well as phonon spectral functions. We review the theoretical framework of the stochastic self-consistent harmonic approximation and its dynamical extension, making particular emphasis on the physical interpretation of the variables present in the theory that can enlighten the comparison with any other anharmonic theory. A modular and flexible Python environment is used for the implementation, which allows for a clean interaction with other packages. We briefly present a toy-model calculation to illustrate the potential of the code.
rate research
Read More
At room temperature, PbTe and SnTe are efficient thermoelectrics with a cubic structure. At low temperature, SnTe undergoes a ferroelectric transition with a critical temperature strongly dependent on the hole concentration, while PbTe is an incipient ferroelectric. By using the stochastic self-consistent harmonic approximation, we investigate the anharmonic phonon spectra and the occurrence of a ferroelectric transition in both systems. We find that vibrational spectra strongly depends on the approximation used for the exchange-correlation kernel in density functional theory. If gradient corrections and the theoretical volume are employed, then the calculation of the free energy Hessian leads to phonon spectra in good agreement with experimental data for both systems. In PbTe, we reproduce the transverse optical mode phonon satellite detected in inelastic neutron scattering and the crossing between the transverse optical and the longitudinal acoustic modes along the $Gamma$X direction. In the case of SnTe, we describe the occurrence of a ferroelectric transition from the high temperature Fm$overline{3}$m structure to the low temperature R3m one.
The self-consistent harmonic approximation is extended in order to account for the existence of Klein factors in bosonized Hamiltonians. This is important for the study of finite systems where Klein factors cannot be ignored a priori. As a test we apply the method to interacting spinless fermions with modulated hopping. We calculate the finite-size corrections to the energy gap and the Drude weight and compare our results with the exact solution for special values of the model parameters.
The vibrational density of states (VDOS) of nanoclusters and nanocrystalline materials are derived from molecular-dynamics simulations using empirical tight-binding potentials. The results show that the VDOS inside nanoclusters can be understood as that of the corresponding bulk system compressed by the capillary pressure. At the surface of the nanoparticles the VDOS exhibits a strong enhancement at low energies and shows structures similar to that found near flat crystalline surfaces. For the nanocrystalline materials an increased VDOS is found at high and low phonon energies, in agreement with experimental findings. The individual VDOS contributions from the grain centers, grain boundaries, and internal surfaces show that, in the nanocrystalline materials, the VDOS enhancements are mainly caused by the grain-boundary contributions and that surface atoms play only a minor role. Although capillary pressures are also present inside the grains of nanocrystalline materials, their effect on the VDOS is different than in the cluster case which is probably due to the inter-grain coupling of the modes via the grain-boundaries.
We present quasiparticle (QP) energies from fully self-consistent $GW$ (sc$GW$) calculations for a set of prototypical semiconductors and insulators within the framework of the projector-augmented wave methodology. To obtain converged results, both finite basis-set corrections and $k$-point corrections are included, and a simple procedure is suggested to deal with the singularity of the Coulomb kernel in the long-wavelength limit, the so called head correction. It is shown that the inclusion of the head corrections in the sc$GW$ calculations is critical to obtain accurate QP energies with a reasonable $k$-point set. We first validate our implementation by presenting detailed results for the selected case of diamond, and then we discuss the converged QP energies, in particular the band gaps, for a set of gapped compounds and compare them to single-shot $G_0W_0$, QP self-consistent $GW$, and previously available sc$GW$ results as well as experimental results.
A theoretical study is reported of the molecular-to-atomic transition in solid hydrogen at high pressure. We use the diffusion quantum Monte Carlo method to calculate the static lattice energies of the competing phases and a density-functional-theory-based vibrational self-consistent field method to calculate anharmonic vibrational properties. We find a small but significant contribution to the vibrational energy from anharmonicity. A transition from the molecular Cmca-12 direct to the atomic I4_1/amd phase is found at 374 GPa. The vibrational contribution lowers the transition pressure by 91 GPa. The dissociation pressure is not very sensitive to the isotopic composition. Our results suggest that quantum melting occurs at finite temperature.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
