Microstructure evolution under thermo-mechanical operating of rocksalt-structure TiN via neural network potential.

In the process of high temperature service, the mechanical properties of cutting tools decrease sharply due to the peeling of the protective coating. However, the mechanism of such coating failure remains obscure due to the complicated interaction between atomic structure, temperature, and stress. This dynamic evolution nature demands both large system sizes and accurate description on the atomic scale, raising challenges for existing atomic scale calculation methods. Here, we developed a deep neural network (DNN) potential for Ti-N binary systems based on first-principles study datasets to achieve quantum-accurate large-scale atomic simulation. Compared with empirical interatomic potential based on the embedded-atom-method, the developed DNN-potential can accurately predict lattice constants, phonon properties, and mechanical properties under various thermodynamic conditions. Moreover, for the first time, we present the atomic evolution of the fracture behavior of large-scale rocksalt-structure (B1) TiN systems coupled with temperature and stress conditions. Our study validates that interatomic brittle fractures occur when TiN stretches beyond its tensile yield point. Such simulation of coating fracture and cutting behavior based on large-scale atoms can shed new light on understanding the microstructure and mechanical properties of coating tools under extreme operating conditions.

The thermodynamic-pathway-determined microstructure evolution of copper under shock compression.

Shock-induced structural transformations in copper exhibit notable directional dependence and anisotropy, but the mechanisms that govern the responses of materials with different orientations are not yet well understood. In this study, we employ large-scale non-equilibrium molecular dynamics simulations to investigate the propagation of a shock wave through monocrystal copper and analyse the structural transformation dynamics in detail. Our results indicate that anisotropic structural evolution is determined by the thermodynamic pathway. A shock along the [Formula: see text] orientation causes a rapid and instantaneous temperature spike, resulting in a solid-solid phase transition. Conversely, a liquid metastable state is observed along the [Formula: see text] orientation due to thermodynamic supercooling. Notably, melting still occurs during the [Formula: see text]-oriented shock, even if it falls below the supercooling line in the thermodynamic pathway. These results highlight the importance of considering anisotropy, the thermodynamic pathway and solid-state disordering when interpreting phase transitions induced by shock. This article is part of the theme issue 'Dynamic and transient processes in warm dense matter'.

Strongly anisotropic ultrafast dynamic behavior of GaTe dominated by the tilted and flat bands.

The anisotropic transport properties of gallium telluride (GaTe) have been reported by several experiments, giving rise to many debates recently. The anisotropic electronic band structure of GaTe shows the extreme difference between the flat band and tilted band in two distinct directions,Γ¯-X¯andΓ¯-Y¯, and which we called as the mixed flat-tilted band (MFTB). Focusing on such two directions, the relaxation of photo-generated carriers has been studied using the non-adiabatic molecular dynamics (NAMD) method to investigate the anisotropic behavior of ultrafast dynamics. The results show that the relaxation lifetime is different in flat band direction and tilted band direction, which is evidence for the existence of anisotropic behavior of the ultrafast dynamic, and such anisotropic behavior comes from the different intensities of electron-phonon coupling of the flat band and tilted band. Furthermore, the ultrafast dynamic behavior is found to be affected strongly by spin-orbit coupling (SOC) and such anisotropic behavior of the ultrafast dynamic can be reversed by SOC. The tunable anisotropic ultrafast dynamic behavior of GaTe is expected to be detected in ultrafast spectroscopy experiments and it may provide a tunable application in nanodevice design. The results may also provide a reference for the investigation of MFTB semiconductors.

Finite-temperature density-functional-theory investigation on the nonequilibrium transient warm-dense-matter state created by laser excitation.

We present a finite-temperature density-functional-theory investigation of the nonequilibrium transient electronic structure of warm dense Li, Al, Cu, and Au created by laser excitation. Photons excite electrons either from the inner shell orbitals or from the valence bands according to the photon energy, and give rise to isochoric heating of the sample. Localized states related to the 3d orbital are observed for Cu when the hole lies in the inner shell 3s orbital. The electrical conductivity for these materials at nonequilibrium states is calculated using the Kubo-Greenwood formula. The change of the electrical conductivity, compared to the equilibrium state, is different for the case of holes in inner shell orbitals or the valence band. This is attributed to the competition of two factors: the shift of the orbital energies due to reduced screening of core electrons, and the increase of chemical potential due to the excitation of electrons. The finite-temperature effect of both the electrons and the ions on the electrical conductivity is discussed in detail. This work is helpful to better understand the physics of laser excitation experiments of warm dense matter.

Erratum: Extremely Low Electron-ion Temperature Relaxation Rates in Warm Dense Hydrogen: Interplay between Quantum Electrons and Coupled Ions [Phys. Rev. Lett. 122, 015001 (2019)].

This corrects the article DOI: 10.1103/PhysRevLett.122.015001.

Extremely Low Electron-ion Temperature Relaxation Rates in Warm Dense Hydrogen: Interplay between Quantum Electrons and Coupled Ions.

Theoretical and computational modeling of nonequilibrium processes in warm dense matter represents a significant challenge. The electron-ion relaxation process in warm dense hydrogen is investigated here by nonequilibrium molecular dynamics using the constrained electron force field (CEFF) method. CEFF evolves wave packets that incorporate dynamic quantum diffraction that obviates the Coulomb catastrophe. Predictions from this model reveal temperature relaxation times as much as three times longer than prior molecular dynamics results based on quantum statistical potentials. Through analyses of energy distributions and mean free paths, this result can be traced to delocalization. Finally, an improved GMS [Gericke, Murillo, and Schlanges, Phys. Rev. E 78, 025401 (2008)PRESCM1539-375510.1103/PhysRevE.78.025401] model is proposed, in which the Coulomb logarithms are in good agreement with CEFF results.

Dynamic electron-ion collisions and nuclear quantum effects in quantum simulation of warm dense matter.

The structural, thermodynamic and transport properties of warm dense matter (WDM) are crucial to the fields of astrophysics and planet science, as well as inertial confinement fusion. WDM refers to the states of matter in a regime of temperature and density between cold condensed matter and hot ideal plasmas, where the density is from near-solid up to ten times solid density, and the temperature between 0.1 and 100 eV. In the WDM regime, matter exhibits moderately or strongly coupled, partially degenerate properties. Therefore, the methods used to deal with condensed matter and isolated atoms need to be properly validated for WDM. It is therefore a big challenge to understand WDM within a unified theoretical description with reliable accuracy. Here, we review the progress in the theoretical study of WDM with state-of-the-art simulations, i.e. quantum Langevin molecular dynamics and first principles path integral molecular dynamics. The related applications for WDM are also included.

Stability and local magnetic moment of bilayer graphene by intercalation: first principles study.

The migration and magnetic properties of the bilayer graphene with intercalation compounds (BGICs) with magnetic elements are theoretically investigated based on first principles study. Firstly, we find that BGICs with transition metals (Sc-Zn) generate distinct magnetic properties. The intercalation with most of the transition metal atoms (TMAs) gives rise to large magnetic moments from 1.0 to 4.0 µ B, which is valuable for the spintronics. Moreover, graphene can protect the intrinsic properties of the intercalated TMAs, which can be important for applications in catalysis. These phenomena can be explained by theory of spd hybridization definitely. Secondly, weak coupling between TMAs and the surroundings indicates the possibility of implementing quantum information processing and generating controlled entanglements. For the possibility of using these materials in ultrafast electronic transistors, spintronics, catalysis, spin qubit and important applications for the extensions of graphene, we believe that BGICs can provide a significant path to synthesize novel materials.

Low toxicity and accumulation of zinc oxide nanoparticles in mice after 270-day consecutive dietary supplementation.

The toxicity and accumulation of zinc oxide nanoparticles (ZnO-NPs), ZnO microparticles (ZnO-MPs) and Zn ions were evaluated after long-term feeding with zinc-replenished food (1600 mg zinc equivalent per kg food) for 270 consecutive days. It was difficult for ZnO-NPs, ZnO-MPs and Zn ions were difficult to pass through the intestine barrier, and most of them were excreted mainly through feces. The distribution results showed that there was no noticeable difference among the distribution profiles of ZnO-NPs, ZnO-MPs and Zn ions in mice. Zn accumulated only in the digestive tract organs after the exposure to all three samples. However, the biomedical parameters and pathological investigations showed liver lesions induced by ZnO-MPs, but fewer by ZnO-NPs or Zn ions. The reason for the remarkably low in vivo toxicity of ZnO-NPs is discussed. Our findings suggest that ZnO-NPs are relatively biocompatible as the nutritional additive at the commonly used dose.

First-principles study on equation of states and electronic structures of shock compressed Ar up to warm dense regime.

The equation of states (EOS) and electronic structures of argon with temperatures from 0.02 eV to 3 eV and densities from 0.5 g/cm(3) to 5.5 g/cm(3) are calculated using the pair potential and many-body potential molecular dynamics and the density functional theory (DFT) molecular dynamics with van der Waals (vdW) corrections. First-principles molecular dynamics is implemented above 2.0 g/cm(3). For the cases of low densities below 3 g/cm(3), we performed pair potential molecular dynamics in order to obtain the ionic configurations, which are used in density functional theory to calculate the EOS and electronic structures. We checked the validity of different methods at different densities and temperatures, showing their behaviors by comparing EOS. DFT without vdW correction works well above 1 eV and 3.5 g/cm(3). Below 1 eV and 2.0 g/cm(3), it overestimates the pressure apparently and results in incorrect behaviors of the internal energy. With vdW corrections, the semi-empirical force-field correction (DFT-D2) method gives consistent results in the whole density and temperature region, and the vdW density functional (vdW-DF2) method gives good results below 2.5 g/cm(3), but it overestimates the pressure at higher densities. The interactions among the atoms are overestimated by the pair potential above 1 eV, and a temperature dependent scaled pair potential can be used to correct the ionic configurations of the pair potential up to 3 eV. The comparisons between our calculations and the experimental multi-shock compression results show that the Hugoniot line of DFT-D2 and DFT tends to give larger pressure than the results of the self-consistent fluid variational theory, and the difference increases with the density. The electronic energy gap exists for all our cases up to 5.5 g/cm(3) and 1 eV. The effect of vdW interactions on the electronic structures are also discussed.

Quantum molecular dynamics study on the structures and dc conductivity of warm dense silane.

The ionic and electronic structures of warm dense silane at the densities of 1.795, 2.260, 3.382, and 3.844 g/cm(3) have been studied with temperatures from 1000 K to 3 eV using quantum molecular dynamics simulations. At all densities, the structures are melted above 1000 K. The matter states are characterized as polymeric from 1000 to 4000 K and become dense plasma states with further increasing temperature to 1 eV. At two lower densities of 1.795 and 2.260 g/cm(3), silane first dissociates and then becomes the polymeric state via a chain state from the initial crystalline structure. At higher densities, however, no dissociation stage was found. These findings can help us understand how the warm dense matter forms. A rise is found for the direct current electric conductivity at T ∼ 1000 K, indicating the nonmetal-to-metal transition. The conductivity decreases slightly with the increase of temperature, which is due to the more disordered structures at higher temperatures.

Nuclear quantum dynamics in dense hydrogen.

Nuclear dynamics in dense hydrogen, which is determined by the key physics of large-angle scattering or many-body collisions between particles, is crucial for the dynamics of planet's evolution and hydrodynamical processes in inertial confinement confusion. Here, using improved ab initio path-integral molecular dynamics simulations, we investigated the nuclear quantum dynamics regarding transport behaviors of dense hydrogen up to the temperatures of 1 eV. With the inclusion of nuclear quantum effects (NQEs), the ionic diffusions are largely higher than the classical treatment by the magnitude from 20% to 146% as the temperature is decreased from 1 eV to 0.3 eV at 10 g/cm(3), meanwhile, electrical and thermal conductivities are significantly lowered. In particular, the ionic diffusion is found much larger than that without NQEs even when both the ionic distributions are the same at 1 eV. The significant quantum delocalization of ions introduces remarkably different scattering cross section between protons compared with classical particle treatments, which explains the large difference of transport properties induced by NQEs. The Stokes-Einstein relation, Wiedemann-Franz law, and isotope effects are re-examined, showing different behaviors in nuclear quantum dynamics.

Quantum simulation of thermally-driven phase transition and oxygen K-edge x-ray absorption of high-pressure ice.

The structure and phase transition of high-pressure ice are of long-standing interest and challenge, and there is still a huge gap between theoretical and experimental understanding. The quantum nature of protons such as delocalization, quantum tunneling and zero-point motion is crucial to the comprehension of the properties of high-pressure ice. Here we investigated the temperature-induced phase transition and oxygen K-edge x-ray absorption spectra of ice VII, VIII and X using ab initio path-integral molecular dynamics simulations. The tremendous difference between experiments and the previous theoretical predictions is closed for the phase diagram of ice below 300 K at pressures up to 110 GPa. Proton tunneling assists the proton-ordered ice VIII to transform into proton-disordered ice VII where only thermal activated proton-transfer cannot occur. The oxygen K edge with its shift is sensitive to the order-disorder transition, and therefore can be applied to diagnose the dynamics of ice structures.

Dynamic ionic clusters with flowing electron bubbles from warm to hot dense iron along the Hugoniot curve.

Complex structures of warm and hot dense matter are essential to understanding the behavior of materials in high energy density processes and provide new features of matter constitutions. Here, around a new unified first-principles determined Hugoniot curve of iron from the normal condensed condition up to 1 Gbar, the novel structures characterized by the ionic clusters with electron bubbles are found using quantum Langevin molecular dynamics. Subsistence of complex clusters can persist in the time scale of 50 fs dynamically with quantum flowing bubbles, which are produced by the interplay of Fermi electron degeneracy, the ionic coupling, and the dynamical nature. With the inclusion of those complicated features in quantum Langevin molecular dynamics, the present equation of states could serve as a first-principles based database in a wide range of temperatures and densities.

Changes of structure and dipole moment of water with temperature and pressure: a first principles study.

The changes of structure and distribution of dipole moment of water with temperatures up to 2800 K and densities up to 2.2 g/cm(3) are investigated using ab initio molecular dynamics. Along the isochore of 1.0 g/cm(3), the structure of liquid water above 800 K is dramatically different from that at ambient conditions, where the hydrogen-bonds network collapses. Along the isotherm of 1800 K, the transition from the liquid state to an amorphous superionic phase occurs at 2.0 g/cm(3) (32.9 GPa), which is not observed along the isotherm of 2800 K. With increasing temperature, the average dipole moment of water molecules is decreased arising from the weakened polarization by the collapse of the hydrogen-bonds network, while it is contrarily increased with compression due to the strengthening effect upon the polarization of water molecules. Both higher temperature and pressure broaden the distribution of dipole moment of water molecules due to the enhanced intramolecular charge fluctuations.

Structure and vibrational spectra of small water clusters from first principles simulations.

The structure and vibrational spectra of (H(2)O)(n) (n=2-5) clusters have been studied based on first-principles molecular dynamics simulations. Trends of the cluster structures with the cluster size show that water molecules in cluster are bound more tightly. The vibrational spectra as a function of cluster size and temperature are obtained using Fourier transformation of the velocity autocorrelation function. Results of the clusters in ground state show that when the cluster size increases, the librational peaks shift to blue and the bonded intramolecular OH stretching bands shift to red due to the clusterization and hydrogen-bond strengthening. Meanwhile, there are no significant shifts in the intramolecular bending and free OH stretching modes, indicating that the free hydrogen atoms are insensitive to the local bonding environment. The temperature-dependent vibrational spectra, which exhibit similar behaviors from the dimer to pentamer, show that there are significant broadenings of the spectra with temperature caused by thermal motions. Moreover, different bands shift to different directions, where librational bands shift to red while bonded OH stretching bands shift to blue, although the blueshifts are quite small for the dimer and trimer.