First-Principles Methods in the Investigation of the Chemical and Transport Properties of Materials under Extreme Conditions

Earth is a dynamic system. The thermodynamics conditions of Earth vary drastically depending on the depth, ranging from ambient temperature and pressure at the surface to 360 GPa and 6600 K at the core. Consequently, the physical and chemical properties of Earth’s constituents (e.g., silicate and carbonate minerals) are strongly affected by their immediate environment. In the past 30 years, there has been a tremendous amount of progress in both experimental techniques and theoretical modeling methods for material characterization under extreme conditions. These advancements have elevated our understanding of the properties of minerals, which is essential in order to achieve full comprehension of the formation of this planet and the origin of life on it. This article reviews recent computational techniques for predicting the behavior of materials under extreme conditions. This survey is limited to the application of the first-principles molecular dynamics (FPMD) method to the investigation of chemical and thermodynamic transport processes relevant to Earth Science.

First-principles study on stability and superconductivity of ternary hydride LaYH_(x)(x=2,3,6 and 8)

Recent studies have shown that the La-and Y-hydrides can exhibit significant superconducting properties under high pressures.In this paper,we investigate the stability,electronic and superconducting properties of LaYH_(x)(x=2,3,6 and 8)under 0-200 GPa.It is found that LaYH_(2) stabilizes in the C2/m phase at ambient pressure,and transforms to the Pmmn phase at 67 GPa.LaYH_(3) stabilizes in the C2/m phase at ambient pressure,and undergoes phase transitions of C2/m→P2_(1)/m→R3m at 12 GPa and 87 GPa,respectively.LaYH_(6) stabilizes in the P4_32_12 phase at ambient pressure,and undergoes phase transitions of P4_(3)2_(1)2→P4/mmm→Cmcm at 28 GPa and 79 GPa,respectively.LaYH_(8) stabilizes in the Imma phase at 60 GPa and transforms to the P4/mmm phase at 117 GPa.Calculations of the electronic band structures show that the P4/mmm-LaYH_(8) and all phases of LaYH_(2) and LaYH_(3) exhibit metallic character.For the metallic phases,we then study their superconducting properties.The calculated superconducting transition temperatures(T_c)are 0.47 K for C2/m-LaYH_(2) at 0 GPa,0 K for C2/m-LaYH_(3) at 0 GPa,and 55.51 K for P4/mmm-LaYH_(8) at 50 GPa.

The growth ofβphase in Mg-Gd-Y-Ni alloy by experimental and first-principles study

The paper reports on the atomic investigation aboutβphase in Mg_(96)Gd_(2)Y_(1)Ni_(1) alloy by using the first-principles study and the high-angle annular dark-field scanning transmission electron microscope(HAADF-STEM)corrected by atomic Cs.By using HAADF-STEM,the rectangularβphases were observed in the underage and peak aging stages in Mg_(96)Gd_(2)Y_(1)Ni_(1) alloy.Theβphase could be precipitated from the previously precipitatedβphase,and theβphase grew in steps when it was precipitated.A special transition structure of three atomic layer thicknesses was first observed at the edge of theβphase and the structure of this interface is probably as theβ/Mg_(1) interface for the analysis of thermodynamic characterization and electronic characterization.Theβ'phase and theβ_(H) structure were precipitated only at the edge of the length directions of theβphase.Theβ'phase continues to grow into aβphase directly without the formation ofβ_(1) phase,resulting in an increase in the length of theβphase,which is discovered for the first time.

Transformation of long-period stacking ordered structures in Mg-Gd-Y-Zn alloys upon synergistic characterization of first-principles calculation and experiment and its effects on mechanical properties

Based on experiments and first-principles calculations,the microstructures and mechanical properties of as-cast and solution treated Mg-10Gd-4Y-xZn-0.6Zr(x=0,1,2,wt.%)alloys are investigated.The transformation process of long-period stacking ordered(LPSO)structure during solidification and heat treatment and its effect on the mechanical properties of experimental alloys are discussed.Results reveal that the stacking faults and 18R LPSO phases appear in the as-cast Mg-10Gd-4Y-1Zn-0.6Zr and Mg-10Gd-4Y-2Zn-0.6Zr alloys,respectively.After solution treatment,the stacking faults and 18R LPSO phase transform into 14H LPSO phase.The Enthalpies of formation and reaction energy of 14H and 18R LPSO are calculated based on first-principles.Results show that the alloying ability of 18R is stronger than that of 14H.The reaction energies show that the 14H LPSO phase is more stable than the 18R LPSO.The elastic properties of the 14H and 18R LPSO phases are also evaluated by first-principles calculations,and the results are in good agreement with the experimental results.The precipitation of LPSO phase improves the tensile strength,yield strength and elongation of the alloy.After solution treatment,the Mg-10Gd-4Y-2Zn-0.6Zr alloy has the best mechanical properties,and its ultimate tensile strength and yield strength are 278.7 MPa and 196.4 MPa,respectively.The elongation of Mg-10Gd-4Y-2Zn-0.6Zr reaches 15.1,which is higher than that of Mg-10Gd-4Y0.6Zr alloy.The improving mechanism of elastic modulus by the LPSO phases and the influence on the alloy mechanical properties are also analyzed.

First-principles study on the diffusion behavior of Cs and I in Cr coating

Cs and I can migrate through fuel-cladding interfaces and accelerate the cladding corrosion process induced by the fuel-cladding chemical interaction.Cr coating has emerged as an important candidate for mitigating this chemical interaction.In this study,first-principles calculations were employed to investigate the diffusion behavior of Cs and I in the Cr bulk and grain boundaries to reveal the microscopic interaction mitigation mechanisms at the fuel-cladding interface.The interaction between these two fission products and the Cr coating were studied systematically,and the Cs and I temperature-dependent diffusion coefficients in Cr were obtained using Bocquet’s oversized solute-atom model and Le Claire’s nine-frequency model,respectively.The results showed that the Cs and I migration barriers were significantly lower than that of Cr,and the Cs and I diffusion coefficients were more than three orders of magnitude larger than the Cr self-diffusion coefficient within the temperature range of Generation-IV fast reactors(below 1000 K),demonstrating the strong penetration ability of Cs and I.Furthermore,Cs and I are more likely to diffuse along the grain boundary because of the generally low migration barriers,indicating that the grain boundary serves as a fast diffusion channel for Cs and I.

First-Principles Investigation of Charge Transfer Mechanism of B-Doped 3C-SiC Semiconductor Material

This study delves into the charge transfer mechanism of boron (B)-doped 3C-SiC through first-principles investigations. We explore the effects of B doping on the electronic properties of 3C-SiC, focusing on a 12.5% impurity concentration. Our comprehensive analysis encompasses structural properties, electronic band structures, and charge density distributions. The optimized lattice constant and band gap energy of 3C-SiC were found to be 4.373 Å and 1.36 eV respectively, which is in agreement with previous research (Bui, 2012;Muchiri et al., 2018). Our results show that B doping narrows the band gap, enhances electrical conductivity, and influences charge transfer interactions. The charge density analysis reveals substantial interactions between B dopants and surrounding carbon atoms. This work not only enhances our understanding of the material’s electronic properties, but also highlights the importance of charge density analysis for characterizing charge transfer mechanisms and their implications in the 3C-SiC semiconductors.

Ideal tensile strength of chromium by first-principles method

The ideal tensile strengths of Cr along [001],[110] and [111] directions were calculated based on the first-principles method.The results show that the ideal tensile strengths are 30.83,37.2 and 35.49 GPa for antiferromagnetic Cr,while they are 33.09,47.15 and38.11 GPa for non-magnetic Cr along [001],[110] and[111] directions,respectively.It is obvious that [001] is the weakest direction.When the loading is applied on the direction [001],the ideal tensile strength is reached before the shear instability for both the anti-ferromagnetic and non-magnetic Cr;thus,Cr fails by cleavage and it is deemed to be intrinsically brittle.Meanwhile,for the antiferromagnetic Cr,the correlation between the magnetic moment and volume was analyzed,and the result shows that the magnetic moment increases with the increase in volume and eventually disappears with the increase in strain.In addition,the density of states in the process of loading was also discussed.

Latent Heat of TB18 Titanium Alloy during β to α Phase Transition by DSC and First-Principles Methods

The phase transition of titanium alloys is sensitive to the heat-treatment procedure,accompanied with the latent heat induced by phase transition.However,the latent heat during phase transition of titanium alloy has not been systematically studied,which may result in the gap between designed and actual temperature of the sample and affect the final mechanical properties.In this work,DSC(differential scanning calorimetry)and first-principles simulate methods were used to study theβ→αphase transition process of TB18(Ti–Al-Mo-V-Cr-Nb-Fe system)metastableβtitanium alloy,especially to reveal the influence of the heating rate on latent heat.The ratio of latent heat to endothermic heat of the sample during temperature rising was introduced to interpret the effect of latent heat to actual temperature.The ratio of latent heat to endothermic heat at 1℃/min is about 15 to 20 times higher than that at 10℃/min.The higher ratio indicates that the latent heat of phase transition has a more significant effect on the temperature,which is related to the temperature range of phase transition and theαvolume fraction.Compared with the heating rate of 1℃/min,theβ→αphase transition takes place at higher temperature and the volume fraction ofαis smaller at 10℃/min.Meanwhile,there is a precipitation free zone between grain boundaryαand intragranularαand the distribution ofαlamellae is heterogeneous when the heating rate is 10℃/min.Both of the experimental and theoretical results suggest that the latent heat of phase transition is the main cause of the temperature fluctuation during heat-treatment process.This work has guiding significance for microstructure optimization affected by temperature,to achieve the desired mechanical properties.

First-principles calculations of Ni–(Co)–Mn–Cu–Ti all-d-metal Heusler alloy on martensitic transformation,mechanical and magnetic properties

The martensitic transformation,mechanical,and magnetic properties of the Ni_(2)Mn_(1.5-x)Cu_(x)Ti_(0.5) (x=0.125,0.25,0.375,0.5) and Ni_(2-y)Co_(y)Mn_(1.5-x)Cu_(x)Ti_(0.5)[(x=0.125,y=0.125,0.25,0.375,0.5) and (x=0.125,0.25,0.375,y=0.625)]alloys were systematically studied by the first-principles calculations.For the formation energy,the martensite is smaller than the austenite,the Ni–(Co)–Mn–Cu–Ti alloys studied in this work can undergo martensitic transformation.The austenite and non-modulated (NM) martensite always present antiferromagnetic state in the Ni_(2)Mn_(1.5-x)Cu_(x)Ti_(0.5) and Ni_(2-y)Co_(y)Mn_(1.5-x)Cu_(x)Ti_(0.5) (y<0.625) alloys.When y=0.625 in the Ni_(2-y)Co_(y)Mn_(1.5-x)Cu_(x)Ti_(0.5) series,the austenite presents ferromagnetic state while the NM martensite shows antiferromagnetic state.Cu doping can decrease the thermal hysteresis and anisotropy of the Ni–(Co)–Mn–Ti alloy.Increasing Mn and decreasing Ti content can improve the shear resistance and normal stress resistance,but reduce the toughness in the Ni–Mn–Cu–Ti alloy.And the ductility of the Co–Cu co-doping alloy is inferior to that of the Ni–Mn–Cu–Ti and Ni–Co–Mn–Ti alloys.The electronic density of states was studied to reveal the essence of the mechanical and magnetic properties.

A first-principles study on remote van der Waals epitaxy through a graphene monolayer on semiconductor substrates

To investigate the mechanism of remote epitaxy, where the overlayer can follow the same crystalline structure as the underlying semiconductor substrate through a thin two-dimensional interlayer, we systematically study the potential fluctuations of graphene covered Si, Ga As, and Ga N substrates from first-principles. We find that the uneven semiconductor surface, the distorted graphene, and the non-uniform interface charge transfer make significant contributions to the potential fluctuation. The semiconductor substrate with different surface reconstructions and orientations will generate different potential fluctuations through the graphene interlayer. We also calculate and compare the adsorption of adatoms on graphene covered substrates. The adsorption energies of adatoms not only depend on their distances to the underlying semiconductor surface, but are also sensitive to the direction of the charge transfer at the graphene/substrate interface. Changing the semiconductor reconstruction or orientation could even reverse the order of the adsorption energies of cation and anion adatoms by reversing the interface charge transfer direction, leading to a change in the growth orientation of the overlayer.Our study improves the understanding of the mechanism of remote epitaxy, and reveals that it is possible to control the initial nucleation and orientation of overlayers by changing the semiconductor reconstructions and/or orientations in remote epitaxy.

First-principles Study of Electronic Structural and Mechanical Properties of Mg_(x)La(x=1,2,3)Compounds under Pressure

The effects of pressure on structural,elastic and electronic properties of Mg_(x)La(x=1,2,3)compounds are investigated by using CASTEP program based on the density functional theory.The calculated equilibrium lattice parameters at zero pressure agree well with available experimental and theoretical values.The calculated DOS show that the structure of these compounds remains mechanically stable and structural phase transformation is not induced with increasing pressure from 0 to 30 GPa,and their structural stability increases with pressure.The ductility of MgLa can be improved by increasing pressure,which is the same as Mg_(2)La in 0-20 GPa,while brittle behavior turns into ductile behavior in 0-5 GPa for Mg_(3)La.The resistance to volume deformation of Mg_(x)La(x=1,2,3)compounds can be improved as the pressure increases.The shear deformation resistance and elastic stiffness of Mg_(3)La can be enhanced by rising pressure,but MgLa and Mg_(2)La increase first and then decrease when pressure is up to 25 GPa.In addition,the three compounds exhibit the elastic anisotropy with pressure.

Pressure-dependent electronic,optical,and mechanical properties of antiperovskite X_(3)NP(X=Ca,Mg):A first-principles study

Hydrostatic pressure provides an efficient way to tune and optimize the properties of solid materials without chang-ing their composition.In this work,we investigate the electronic,optical,and mechanical properties of antiperovskite X_(3)NP(X^(2+)=Ca,Mg)upon compression by first-principles calculations.Our results reveal that the system is anisotropic,and the lat-tice constant a of X_(3)NP exhibits the fastest rate of decrease upon compression among the three directions,which is different from the typical Pnma phase of halide and chalcogenide perovskites.Meanwhile,Ca_(3)NP has higher compressibility than Mg_(3)NP due to its small bulk modulus.The electronic and optical properties of Mg_(3)NP show small fluctuations upon compression,but those of Ca_(3)NP are more sensitive to pressure due to its higher compressibility and lower unoccupied 3d orbital energy.For example,the band gap,lattice dielectric constant,and exciton binding energy of Ca_(3)NP decrease rapidly as the pressure increases.In addition,the increase in pressure significantly improves the optical absorption and theoretical conversion effi-ciency of Ca_(3)NP.Finally,the mechanical properties of X_(3)NP are also increased upon compression due to the reduction in bond length,while inducing a brittle-to-ductile transition.Our research provides theoretical guidance and insights for future experi-mental tuning of the physical properties of antiperovskite semiconductors by pressure.

Rational design of Fe/Co-based diatomic catalysts for Li–S batteries by first-principles calculations

Fe/Co-based diatomic catalysts decorated on an N-doped graphene substrate are investigated by first-principles calculations to improve the electrochemical properties of Li–S batteries.Our results demonstrate that Fe CoN8@Gra not only possesses moderate adsorption energies towards Li2Snspecies,but also exhibits superior catalytic activity for both reduction and oxidation reactions of the sulfur cathode.Moreover,the metallic property of the diatomic catalysts can be well maintained after Li2Snadsorption,which could help the sulfur cathode to maintain high conductivity during the whole charge–discharge process.Given these exceptional properties,it is expected that Fe CoN8@Gra could be a promising diatomic catalyst for Li–S batteries and afford insights for further development of advanced Li–S batteries.

First-principles study of the bandgap renormalization and optical property ofβ-LiGaO_(2)

Theβ-LiGaO_(2)with an orthorhombic wurtzite-derived structure is a candidate ultrawide direct-bandgap semiconductor.In this work,using the non-adiabatic Allen-Heine-Cardona approach,we investigate the bandgap renormalization arising from electron-phonon coupling.We find a sizable zero-point motion correction of-0.362 eV to the gap atΓ,which is dominated by the contributions of long-wavelength longitudinal optical phonons.The bandgap ofβ-LiGaO_(2)decreases monotonically with increasing temperature.We investigate the optical spectra by comparing the model Bethe-Salpether equation method with the independent-particle approximation.The calculated optical spectra including electron-hole interactions exhibit strong excitonic effects,in qualitative agreement with the experiment.The contributing interband transitions and the binding energy for the excitonic states are analyzed.

Thermal transport properties of two-dimensional boron dichalcogenides from a first-principles and machine learning approach

The investigation of thermal transport is crucial to the thermal management of modern electronic devices.To obtain the thermal conductivity through solution of the Boltzmann transport equation,calculation of the anharmonic interatomic force constants has a high computational cost based on the current method of single-point density functional theory force calculation.The recent suggested machine learning interatomic potentials(MLIPs)method can avoid these huge computational demands.In this work,we study the thermal conductivity of two-dimensional MoS_(2)-like hexagonal boron dichalcogenides(H-B_(2)VI_(2);V I=S,Se,Te)with a combination of MLIPs and the phonon Boltzmann transport equation.The room-temperature thermal conductivity of H-B_(2)S_(2)can reach up to 336 W·m^(-1)·K^(-1),obviously larger than that of H-B_(2)Se_(2)and H-B_(2)Te_(2).This is mainly due to the difference in phonon group velocity.By substituting the different chalcogen elements in the second sublayer,H-B_(2)VIV I′have lower thermal conductivity than H-B_(2)VI_(2).The room-temperature thermal conductivity of B2STe is only 11%of that of H-B_(2)S_(2).This can be explained by comparing phonon group velocity and phonon relaxation time.The MLIP method is proved to be an efficient method for studying the thermal conductivity of materials,and H-B_(2)S_(2)-based nanodevices have excellent thermal conduction.

First-principles prediction of quantum anomalous Hall effect in two-dimensional Co2Te lattice

The quantum anomalous Hall effect(QAHE) has special quantum properties that are ideal for possible future spintronic devices. However, the experimental realization is rather challenging due to its low Curie temperature and small non-trivial bandgap in two-dimensional(2D) materials. In this paper, we demonstrate through first-principles calculations that monolayer Co2Te material is a promising 2D candidate to realize QAHE in practice. Excitingly, through Monte Carlo simulations, it is found that the Curie temperature of single-layer Co2Te can reach 573 K. The band crossing at the Fermi level in monolayer Co2Te is opened when spin–orbit coupling is considered, which leads to QAHE with a sizable bandgap of Eg= 96 me V, characterized by the non-zero Chern number(C = 1) and a chiral edge state. Therefore, our findings not only enrich the study of quantum anomalous Hall effect, but also broaden the horizons of the spintronics and topological nanoelectronics applications.

Initial micro-galvanic corrosion behavior between Mg_(2)Ca and α-Mg via quasi-in situ SEM approach and first-principles calculation

The initial micro-galvanic corrosion behavior of Mg-30wt%Ca alloy only containing Mg_(2)Ca andα-Mg was studied by immersion testing in a 0.9%Na Cl solution at 37°C.The quasi-in situ SEM and TEM results show that Mg_(2)Ca corroded easier thanα-Mg,indicating that Mg_(2)Ca acted as an anode.The work function(Φ)for Mg_(2)Ca calculated by first-principles is significantly lower compared to that forα-Mg.The Volta potential measured by a scanning Kelvin probe force microscope reveals that the Mg_(2)Ca had a relatively low Volta potential(ψ)value.The lowerΦandψvalues for Mg_(2)Ca indicate a lower electrochemical nobility,which is consistent with the experimental phenomenon.

First-principles study of the influences of oxygen defects upon the electronic properties of Nb-doped TiO_2 by GGA + U methods

The influence of oxygen defects upon the electronic properties of Nb-doped TiO2 has been studied by using the general gradient approximation （GGA）＋U method. Four independent models （i.e., an undoped anatase cell, an anatase cell with a Nb dopant at Ti site （NbTi）, an anatase cell with a Nb-dopant and an oxygen vacancy （NbTi＋Vo）, and an anatase cell with a Nb-dopant and an interstitial oxygen （NbTi＋Oi）） were considered. The density of states, effective mass, Bader charge, charge density, and electron localization function were calcul^ited. The results show that in the NbTi＋Vo cell both eg and t2g levels of Ti 3d orbits make contributions to the electronic conductivity, and the oxygen vacancies （Vo） collaborate with Nb-dopants to favor the high electrical conductivity by inducing the Nb-dopants to release more excess charges. In NbTi＋Oi, an unoccupied impurity level appears in the band gap, which served as an acceptor level and suppressed the electronic conductivity. The results qualitatively coincide with experimental results and possibly provide insights into the preparation of TCOs with desirable conductivity.

A comparison study on structure-function relationship of polysaccharides obtained from sea buckthorn berries using different methods:antioxidant and bile acid-binding capacity

In this study,the structural characters,antioxidant activities and bile acid-binding ability of sea buckthorn polysaccharides(HRPs)obtained by the commonly used hot water(HRP-W),pressurized hot water(HRP-H),ultrasonic(HRP-U),acid(HRP-C)and alkali(HRP-A)assisted extraction methods were investigated.The results demonstrated that extraction methods had significant effects on extraction yield,monosaccharide composition,molecular weight,particle size,triple-helical structure,and surface morphology of HRPs except for the major linkage bands.Thermogravimetric analysis showed that HRP-U with filamentous reticular microstructure exhibited better thermal stability.The HRP-A with the lowest molecular weight and highest arabinose content possessed the best antioxidant activities.Moreover,the rheological analysis indicated that HRPs with higher galacturonic acid content and molecular weight showed higher viscosity and stronger crosslinking network(HRP-C,HRP-W and HRP-U),which exhibited stronger bile acid binding capacity.The present findings provide scientific evidence in the preparation technology of sea buckthorn polysaccharides with good antioxidant and bile acid binding capacity which are related to the structure affected by the extraction methods.

Drilling-based measuring method for the c-φ parameter of rock and its field application

The technology of drilling tests makes it possible to obtain the strength parameter of rock accurately in situ. In this paper, a new rock cutting analysis model that considers the influence of the rock crushing zone(RCZ) is built. The formula for an ultimate cutting force is established based on the limit equilibrium principle. The relationship between digital drilling parameters(DDP) and the c-φ parameter(DDP-cφ formula, where c refers to the cohesion and φ refers to the internal friction angle) is derived, and the response of drilling parameters and cutting ratio to the strength parameters is analyzed. The drillingbased measuring method for the c-φ parameter of rock is constructed. The laboratory verification test is then completed, and the difference in results between the drilling test and the compression test is less than 6%. On this basis, in-situ rock drilling tests in a traffic tunnel and a coal mine roadway are carried out, and the strength parameters of the surrounding rock are effectively tested. The average difference ratio of the results is less than 11%, which verifies the effectiveness of the proposed method for obtaining the strength parameters based on digital drilling. This study provides methodological support for field testing of rock strength parameters.