Pressure Tuned Structural, Electronic and Elastic Properties of U 3 Si 2 C 2 : A First Principles Study

: U 3 Si 2 C 2 is expected to be a new nuclear fuel as a ternary compound of uranium, silicon and carbon. However, the relevant research on U 3 Si 2 C 2 under accident conditions is rarely reported. Hence it is necessary to explore the service behavior of the potential U-Si-C ternary nuclear fuel in extreme environments. In this work, the structural characteristics, electronic behaviors and mechanical properties of U 3 Si 2 C 2 , such as stable crystalline structures, density of states, charge distributions, electron localization function, electronic thermal conductivity and elastic modulus under extreme high pressure are calculated by density functional theory. The calculation results show that the lattice volume sharply increases when the external stress reached 9.8 GPa. Ionic and metallic nature coexist as to the bonding characteristics of U 3 Si 2 C 2 , and the ionic takes the dominant position in bonding. The toughness of U 3 Si 2 C 2 is predicted to be better compared to U 3 Si 2 . Our theoretical investigation may help with the application of U 3 Si 2 C 2 -based fuel and the design of ternary uranium fuels.


Methodology
The DFT theoretical calculations in this paper are carried out using Vienna ab initio simulations package (VASP) [37,38]. During the calculation process, the projectedaugmented-wave potential (PAW) [39] is adopted with the cutoff energy of 520 eV. The scheme by Perdew-Burke-Ernzerhof (PBE) [40] is functional with generalized gradient approximation (GGA) can be used to describe the exchange and correlation interactions between electrons. The strong on-site Coulomb interaction among the localized U 5f electrons is eliminated by Hubbard U approximation (GGA + U) [41], which has been proved effective in our previous works on U-Si and U-Si-Al compounds [19,20,42]. A 8 × 8 × 4 γ-centered Monkhorst-Pack [43] k-point grid is adopted for the unit cell for Brillouin zone sampling. The optimizations of the structural parameters under different pressures from 0 to 24 GPa are performed through conjugate-gradient algorithm [44] in the condition that the force on the atoms is less than 0.001 eV/Å and the total energy difference is smaller than 1.0 × 10 −6 eV/cell. The theoretical phonon spectrum is obtained according to the density functional perturbation theory (DFPT) [45,46] with 2 × 2 × 1 supercell, and the phonon band structures are calculated by the Phonopy package [47]. The ab-initio molecular dynamics (AIMD) simulations for 56-atom supercell under NPT ensemble from 300 to 2100 K with a time step of 1 fs can be employed to determine the structural stability under extreme temperatures.
The mechanical anisotropy is evaluated by universal anisotropic index A U [48] in this work, which is described as Equation (1): For crystals with tetragonal symmetry, [49][50][51], the B V and B R are the bulk modulus following Voigt and Reuss approximations respectively, G V and G R are the corresponding shear modulus. These four variables can be expressed as Equations (2)- (7): (2) In which S ij represents the elastic compliance constants, and C ij is the elastic constants. The mechanical anisotropy is also evaluated by three-dimensional (3D) Young's modulus. The directional dependent Young's moduli for tetragonal crystal system can be calculated by Equation (8): l 1 , l 2 , and l 3 donate the cosines between the given vector and a, b, and c axis. In the tetragonal crystal, 6 independent compliance constants S ij are available. The Young's modulus E and Poisson's ratio ν [49] are expressed as follows:

Ground State Structural Parameters
The original structure of U 3 Si 2 C 2 is obtained from Pöttgen et al. [18]. The U 3 Si 2 C 2 compound crystallizes in the structure of tetragonal, belonging to the space group of I4/mmm, which has 14 atoms in the conventional unit cell and an exact 3:2:2 stoichiometry, as shown in Figure 1. The types of uranium atoms can be divided into two (U1 and U2) [19] by their nearest neighbor atoms. The structural optimization of U 3 Si 2 C 2 are implemented with GGA and GGA + U methods, respectively. In order to study the ground-state properties of U 3 Si 2 C 2 , the common DFT + U method is adopted for structural optimization, which has been widely used for other nuclear fuels such as UO 2 [52], U-Si [36,42], and U-C [32,53]. The influence of Hubbard U parameter values from 0 to 4 eV on the lattice constants (a and c) and the volumes (V) of U 3 Si 2 C 2 unit cells are listed in Table 1. The values of a, c and V calculated by using PAW-GGA potential are 3.645 Å, 16.790 Å, and 222.24 Å 3 . The volume obtained from PAW-GGA potential underestimates severely compared with the experimental [30] and DFT + U [31] results. However, as shown in Table 1, the optimized volume using GGA + U (U = 3.5 eV) is in better agreement with the experimental data. Therefore, 3.5 eV is chosen as the calculation parameter of U in following studies.

Stability of Ground State Structure
The dynamical stability of U 3 Si 2 C 2 has been confirmed by phonon dispersion curves and its phonon density of states. As shown in Figure 2, the absence of imaginary frequency appearing in the Brillouin zone demonstrated the structure is dynamically stable at 0 K. C sublattice occupied higher vibration frequency bands in the phonon dispersion, the value of which are in the range of 10.0-12.0 THz and 18.6-20.4 THz. The vibration frequency of phonon dispersion is determined by the relative atomic mass; hence the low vibration frequency bands (below 5 THz) are mainly originated from the vibrations of the metal U sublattice. A maximum peak can be observed at the frequency of 9.03 THz in the phonon density of states (p-DOS) in Figure 2, indicating that the lattice vibration of the U 3 Si 2 C 2 crystal is mainly concentrated on this frequency. The mechanically stability of U3Si2C2 can be judged by elastic constants. For the tetragonal (I) symmetry, there exists six independent elastic constants, namely C11, C12, C13, C33, C44, and C66 [54], and the calculated elastic constants are listed in Error! Reference source not found.. All elastic constants obey the elastic stability criteria list as the Equation (11), proving this specific crystal is mechanically stable.  In order to evaluate the thermodynamically stability of U3Si2C2, the average formation energy (ΔE) has been calculated as the following Equation (12): where Etotal is the total energy of the calculated U3Si2C2 unit cell, NU, NSi, and NC represent the number of atoms in the cell, and EU, ESi, and EC denote the energy of single atom of bulk U (Cmcm), Si (Fd3 _ m), and C (R3 _ m). In our calculation, the average ΔE have a negative value of −0.252 eV/atom, which indicates that the structure shows thermodynamical stability.
Besides, as an important performance metric of nuclear materials, the thermal stability has also been evaluated by Ab initio molecular dynamics (AIMD) simulations over supercells composite of 56 atoms. The simulation with a total duration of 10,000 fs is performed from 300 K to 2200 K with a time step of 1 fs. The peritectic reaction observed at about 1700 °C in the former research [55][56][57] illustrates that the phase of U3Si2C2 decompose into UC and a liquid phase of composition close to USi2 by taking the phase diagram into consideration. As shown in Error! Reference source not found.a, the average energy of U3Si2C2 is almost unchanged below 1973 K, and which emerge a steep drop at the time of 7.0 ps and 5.3 ps in the case of 1973 K and 2200 K. In Error! Reference source not The mechanically stability of U 3 Si 2 C 2 can be judged by elastic constants. For the tetragonal (I) symmetry, there exists six independent elastic constants, namely C 11 , C 12 , C 13 , C 33 , C 44 , and C 66 [54], and the calculated elastic constants are listed in Table 2. All elastic constants obey the elastic stability criteria list as the Equation (11), proving this specific crystal is mechanically stable. In order to evaluate the thermodynamically stability of U 3 Si 2 C 2 , the average formation energy (∆E) has been calculated as the following Equation (12): where E total is the total energy of the calculated U 3 Si 2 C 2 unit cell, N U , N Si , and N C represent the number of atoms in the cell, and E U , E Si , and E C denote the energy of single atom of bulk U (Cmcm), Si (Fd3m), and C (R3m). In our calculation, the average ∆E have a negative value of −0.252 eV/atom, which indicates that the structure shows thermodynamical stability. Besides, as an important performance metric of nuclear materials, the thermal stability has also been evaluated by Ab initio molecular dynamics (AIMD) simulations over supercells composite of 56 atoms. The simulation with a total duration of 10,000 fs is performed from 300 K to 2200 K with a time step of 1 fs. The peritectic reaction observed at about 1700 • C in the former research [55][56][57] illustrates that the phase of U 3 Si 2 C 2 decompose into UC and a liquid phase of composition close to USi 2 by taking the phase diagram into consideration. As shown in Figure 3a, the average energy of U 3 Si 2 C 2 is almost unchanged below 1973 K, and which emerge a steep drop at the time of 7.0 ps and 5.3 ps in the case of 1973 K and 2200 K. In Figure 3b, the bond length of C-Si stabilizes at 1.9 Å below 1973 K, while the curves rise sharply and oscillate continuously at 1973 K and 2200 K. The curves of energy and bond length indicate that the structure of U 3 Si 2 C 2 is unstable in the environ-ment of exceeding 1973 K, reflected in the destruction of the C-Si bond. The instability of C-Si bond above 1973 K can also be illustrated by the absences of the peak at around 1.9 Å in the radial distribution function curves shown in Figure 3c. found.b, the bond length of C-Si stabilizes at 1.9 Å below 1973 K, while the curves rise sharply and oscillate continuously at 1973 K and 2200 K. The curves of energy and bond length indicate that the structure of U3Si2C2 is unstable in the environment of exceeding 1973 K, reflected in the destruction of the C-Si bond. The instability of C-Si bond above 1973 K can also be illustrated by the absences of the peak at around 1.9 Å in the radial distribution function curves shown in Error! Reference source not found.c.

Structural Parameters and Electronic Behaviors
The structural optimization from 0 GPa to 24 GPa have been carried out by adding isotropic external stress to stress tensor in VASP code [37,38]. Error! Reference source not found. shows the variations of lattice parameters (a and c), cell volume (V) and c/a as a function of pressure. The lattice constants and volumes decrease as the pressure increases at the beginning, moreover, the decrease of lattice parameter c is faster than that of parameter a. However, the lattice parameter c increased dramatically from 17.44 Å to 18.02 Å when the external stress reached 9.8 GPa, resulting in a sudden increase in the cell volume. To explore the mechanism for the increase of lattice parameters, the ratio of the interatomic distance for the U1-U2, U2-C, C-Si, and Si-Si atom pairs along the c direction under high pressure to that under 0 GPa have been calculated (Error! Reference source not found.). The value of distance for the Si-C pair decrease slower compared to other atom pairs when the external pressure is below 9.8 GPa. While the external pressure exceeds 9.8 GPa, the Si-C distance reducing by 0.7% shows obviously opposite trend compared to the other atom pairs, which may be explained by the formation of covalent bond between Si and C atoms. When the external pressure reaches 9.8 GPa, the distance of U1-U2 almost changes to the original length (0 GPa), and the Si-Si distance is even 0.058 Å higher. In order to study the sharply variation of c-axis and volume mutations, the density of states (DOS) under different external pressures, including total density of states (TDOS) and the partial density of states (PDOS) of U3Si2C2, have been calculated to analyze its electronic structure.

Structural Parameters and Electronic Behaviors
The structural optimization from 0 GPa to 24 GPa have been carried out by adding isotropic external stress to stress tensor in VASP code [37,38]. The lattice parameters of the optimized cells under different pressures are listed in the supplementary materials (Table S1 and Figure S1). Figure 4 shows the variations of lattice parameters (a and c), cell volume (V) and c/a as a function of pressure. The lattice constants and volumes decrease as the pressure increases at the beginning, moreover, the decrease of lattice parameter c is faster than that of parameter a. However, the lattice parameter c increased dramatically from 17.44 Å to 18.02 Å when the external stress reached 9.8 GPa, resulting in a sudden increase in the cell volume. To explore the mechanism for the increase of lattice parameters, the ratio of the interatomic distance for the U1-U2, U2-C, C-Si, and Si-Si atom pairs along the c direction under high pressure to that under 0 GPa have been calculated ( Figure 5). The value of distance for the Si-C pair decrease slower compared to other atom pairs when the external pressure is below 9.8 GPa. While the external pressure exceeds 9.8 GPa, the Si-C distance reducing by 0.7% shows obviously opposite trend compared to the other atom pairs, which may be explained by the formation of covalent bond between Si and C atoms. When the external pressure reaches 9.8 GPa, the distance of U1-U2 almost changes to the original length (0 GPa), and the Si-Si distance is even 0.058 Å higher. In order to study the sharply variation of c-axis and volume mutations, the density of states (DOS) under different external pressures, including total density of states (TDOS) and the partial density of states (PDOS) of U 3 Si 2 C 2 , have been calculated to analyze its electronic structure.   The calculated DOS curves are shown in Error! Reference source not found.. It can be seen that the state density curves at 0 GPa in Error! Reference source not found.a pass through the Fermi energy level, indicating U3Si2C2 has metallic properties. The curve for the total density of states displays slight asymmetry, indicating magnetic moment existing in the lattice which is mainly originating from the 5f electrons of uranium atoms (at both U1 and U2 sites) as shown in the partial density of states. Strong interactions between the s orbitals of Si atoms and C atoms can be observed around the peaks at the energy of −10.5 eV and −6.7 eV. The intense peak of 5f states at the Fermi level suggests delocalized nature of uranium 5f electrons. The value of magnetic moments of atoms at U1 and U2 sites are   The calculated DOS curves are shown in Error! Reference source not found.. It can be seen that the state density curves at 0 GPa in Error! Reference source not found.a pass through the Fermi energy level, indicating U3Si2C2 has metallic properties. The curve for the total density of states displays slight asymmetry, indicating magnetic moment existing in the lattice which is mainly originating from the 5f electrons of uranium atoms (at both U1 and U2 sites) as shown in the partial density of states. Strong interactions between the s orbitals of Si atoms and C atoms can be observed around the peaks at the energy of −10.5 eV and −6.7 eV. The intense peak of 5f states at the Fermi level suggests delocalized nature of uranium 5f electrons. The value of magnetic moments of atoms at U1 and U2 sites are The calculated DOS curves are shown in Figure 6, and more information can be seen in supplementary Figure S2. It can be seen that the state density curves at 0 GPa in Figure 6a pass through the Fermi energy level, indicating U 3 Si 2 C 2 has metallic properties. The curve for the total density of states displays slight asymmetry, indicating magnetic moment existing in the lattice which is mainly originating from the 5f electrons of uranium atoms (at both U1 and U2 sites) as shown in the partial density of states. Strong interactions between the s orbitals of Si atoms and C atoms can be observed around the peaks at the energy of −10.5 eV and −6.7 eV. The intense peak of 5f states at the Fermi level suggests delocalized nature of uranium 5f electrons. The value of magnetic moments of atoms at U1 and U2 sites are −0.61 µ B and 2.08 µ B , respectively. The TDOS and PDOS at 8 GPa vary slightly compared to that at 0 GPa; however, when the pressure increases to 12 GPa, the density of states of U1 atoms near the Fermi level is close to 0, indicating that its metallicity is sharply weakened (showing in Figure 6c). Besides, the spin-down peak of U1-5f orbital in locating at −2 eV disappeared, and a spin-up peak locating at −1.6 eV appeared. Therefore, the magnetic moment of the U1 atom turns positive under the pressure of 12 GPa, and U1 atom has the same direction of magnetic moments as the U2 atom. In conclusion, the c-axis and volume mutation mentioned above may be caused by charge redistributions due to the external pressure, and the specific mechanism and impact will be reported in future work. −0.61 μB and 2.08 μB, respectively. The TDOS and PDOS at 8 GPa vary slightly compared to that at 0 GPa; however, when the pressure increases to 12 GPa, the density of states of U1 atoms near the Fermi level is close to 0, indicating that its metallicity is sharply weakened (showing in Error! Reference source not found.c). Besides, the spin-down peak of U1-5f orbital in locating at −2eV disappeared, and a spin-up peak locating at −1.6eV appeared. Therefore, the magnetic moment of the U1 atom turns positive under the pressure of 12 GPa, and U1 atom has the same direction of magnetic moments as the U2 atom. In conclusion, the c-axis and volume mutation mentioned above may be caused by charge redistributions due to the external pressure, and the specific mechanism and impact will be reported in future work. The projected band-structures of U atoms (divided into U1 site and U2 site) of U3Si2C2 at 0 GPa and 12 GPa have been calculated as shown in Error! Reference source not found.. With the pressure increasing to 12 GPa, the contribution of uranium atoms around the Fermi level gradually decreases, and the result is consistent with the previous PDOS analysis. Moreover, the contribution around the Fermi level of U1 atoms in the spin-down energy band under 12 GPa significantly reduced compared to U2 atoms. In addition, the band structures at Γ point without external pressure shows that the spin-up electrons have almost no contribution at Fermi surface, while an electron pocket formed by the spindown electrons which are mainly contributed by the U1 atoms appears. Furthermore, multiple-electron pockets located at Fermi surface formed by spin-up electrons can be observed at Γ point under 12 GPa, which are dominated by both U1 and U2 atoms, while the spin-down electron pocket are mainly contributed by U2 atoms. The projected band-structures of U atoms (divided into U1 site and U2 site) of U 3 Si 2 C 2 at 0 GPa and 12 GPa have been calculated as shown in Figure 7. With the pressure increasing to 12 GPa, the contribution of uranium atoms around the Fermi level gradually decreases, and the result is consistent with the previous PDOS analysis. Moreover, the contribution around the Fermi level of U1 atoms in the spin-down energy band under 12 GPa significantly reduced compared to U2 atoms. In addition, the band structures at Γ point without external pressure shows that the spin-up electrons have almost no contribution at Fermi surface, while an electron pocket formed by the spin-down electrons which are mainly contributed by the U1 atoms appears. Furthermore, multiple-electron pockets located at Fermi surface formed by spin-up electrons can be observed at Γ point under 12 GPa, which are dominated by both U1 and U2 atoms, while the spin-down electron pocket are mainly contributed by U2 atoms. Properties related to electron localization obtaining from the electron localization function (ELF) [58] showing in Error! Reference source not found. is calculated to reveal atomic bonding of U3Si2C2. In Error! Reference source not found.a,b, the formation of C-Si covalent bonds can be revealed that the charge is mainly concentrated between C and Si atoms, and the electrons of Si atoms are more delocalized as to that of C atoms. The distribution of electrons at atomic gaps suggests a metallic behavior of the compound. The crown-shaped charge accumulations between U and Si implies the U−Si bond is biased towards the ionic bond, and similar structure has been reported in U3Si2. In our previous research, the Si-Si bond is covalent in U3Si2 [19,20,42], however, no similar Si-Si bond can be observed in U3Si2C2 structure. At the external pressure of 12 GPa, the degree of electronic localization at the U1 site displays more anisotropic, as shown in Error! Reference source not found.c,d. The differences concerning the charge of U3Si2C2 structure under high pressures may cause some variations on the properties associating with the electronic structure closely, such as thermal conductivity and elastic constants. Properties related to electron localization obtaining from the electron localization function (ELF) [58] showing in Figure 8 is calculated to reveal atomic bonding of U 3 Si 2 C 2 . In Figure 8a,b, the formation of C-Si covalent bonds can be revealed that the charge is mainly concentrated between C and Si atoms, and the electrons of Si atoms are more delocalized as to that of C atoms. The distribution of electrons at atomic gaps suggests a metallic behavior of the compound. The crown-shaped charge accumulations between U and Si implies the U−Si bond is biased towards the ionic bond, and similar structure has been reported in U 3 Si 2 . In our previous research, the Si-Si bond is covalent in U 3 Si 2 [19,20,42], however, no similar Si-Si bond can be observed in U 3 Si 2 C 2 structure. At the external pressure of 12 GPa, the degree of electronic localization at the U1 site displays more anisotropic, as shown in Figure 8c,d. The differences concerning the charge of U 3 Si 2 C 2 structure under high pressures may cause some variations on the properties associating with the electronic structure closely, such as thermal conductivity and elastic constants.

Thermal and Mechanical Properties
Thermal conductivity is critical to nuclear fuel because it determines the temperature gradients between the surface and centerline of fuel, which may cause thermal stress and make the nuclear fuel deform and generate cracks in the cladding. The electronic thermal conductivity (Ke/τ) is obtained by solving Boltzmann transport equations which can be achieved in BoltzTraP2 [59] as shown in Error! Reference source not found.. The curves illustrate that the electronic thermal conductivity showed positive correlated with temperature, besides, the applied pressure can also influence the electronic thermal conductivity. However, when the applied pressure is more than 10 GPa, the slope of the curve remains basically unchanged, implying that the applied pressure may not be the major factor to affect the electronic thermal conductivity in the condition that the pressure exceeds 10 GPa.

Thermal and Mechanical Properties
Thermal conductivity is critical to nuclear fuel because it determines the temperature gradients between the surface and centerline of fuel, which may cause thermal stress and make the nuclear fuel deform and generate cracks in the cladding. The electronic thermal conductivity (Ke/τ) is obtained by solving Boltzmann transport equations which can be achieved in BoltzTraP2 [59] as shown in Figure 9. The curves illustrate that the electronic thermal conductivity showed positive correlated with temperature, besides, the applied pressure can also influence the electronic thermal conductivity. However, when the applied pressure is more than 10 GPa, the slope of the curve remains basically unchanged, implying that the applied pressure may not be the major factor to affect the electronic thermal conductivity in the condition that the pressure exceeds 10 GPa. Crystals 2021, 11, x FOR PEER REVIEW 11 of 15 The mechanical anisotropy properties of the nuclear materials may cause the deformation and fractures imposed with external stress. In this work, the mechanical anisotropy is evaluated by universal anisotropic index A U and three-dimensional (3D) Young's modulus (E). The universal anisotropic index A U is used to evaluate the mechanical anisotropy, which is suitable for describing single crystals of various types [53]. When A U = 0, it is shown that the crystal is isotropic, and the difference between its value and 0 can be used to evaluate the degree of anisotropy. The calculated Young's modulus (E), bulk modulus (B), shear modulus (G), and Poisson's ratio under different pressures are also shown in Error! Reference source not found. to evaluate the mechanical properties of U3Si2C2. The Poisson's ratio is associated with the properties of interatomic bonding, the value of Poisson's ratio is small (ν = 0.1) for covalent materials, and for ionic crystals it reaches 0.25, whereas for metallic materials it is 0.33 typically [60]. According to the calculated Poisson's ratio, ionic and metallic nature coexist as to the bonding characteristics of U3Si2C2, moreover, the ionic takes the dominant position in bonding. In our previous research, bonding behaviors of U3Si2 (ν = 0.20) [42] is mainly reflected on ionic nature, and U3Si2C2 shows better toughness compared to U3Si2, which may provide theoretical support for the application of U3Si2C2 in nuclear reactors. As the pressure increases to 9.8 GPa, the atomic distance between U1 and U2 decreases from 3.249 Å to 3.098 Å, which may lead to the strength of the metal bond between uranium atoms. When the external pressure applied from 0 GPa to 8 GPa, the Young's modulus, shear modulus and bulk modulus of U3Si2C2 increase 37.8, 34.5 and 14.2 GPa respectively, then decrease sharply as the external pressure exceeds 9.8 GPa, and the changing trend of which seems to be consistent with the variation of the crystal lattice constant. As shown in Error! Reference source not found., the three-dimensional Young's modulus distribution of U3Si2C2 at different pressures and the general anisotropy index A U indicate that with the increase of hydrostatic pressure, the maximum Young's modulus increases gradually. The maximum Young's modulus of U3Si2C2 under different pressures are all appear along the [110] direction. The maximum value of 3D Young's modulus is 252 GPa under 0 GPa pressure, and raising to 317 GPa, 289.26 GPa, and 319.59 GPa, respectively, under the external pressure 8 GPa, 12 GPa, and 20 GPa. Under high-pressure environment, the value of A U for U3Si2C2 increases to higher than 0.8, while that remains basically unchanged when the pressure exceeds 8 GPa. Hence as the pressure increases, the degree of isotropy for U3Si2C2 is reduced slightly. The mechanical anisotropy properties of the nuclear materials may cause the deformation and fractures imposed with external stress. In this work, the mechanical anisotropy is evaluated by universal anisotropic index A U and three-dimensional (3D) Young's modulus (E). The universal anisotropic index A U is used to evaluate the mechanical anisotropy, which is suitable for describing single crystals of various types [53]. When A U = 0, it is shown that the crystal is isotropic, and the difference between its value and 0 can be used to evaluate the degree of anisotropy. The calculated Young's modulus (E), bulk modulus (B), shear modulus (G), and Poisson's ratio under different pressures are also shown in Table 3 to evaluate the mechanical properties of U 3 Si 2 C 2 . The Poisson's ratio is associated with the properties of interatomic bonding, the value of Poisson's ratio is small (ν = 0.1) for covalent materials, and for ionic crystals it reaches 0.25, whereas for metallic materials it is 0.33 typically [60]. According to the calculated Poisson's ratio, ionic and metallic nature coexist as to the bonding characteristics of U 3 Si 2 C 2 , moreover, the ionic takes the dominant position in bonding. In our previous research, bonding behaviors of U 3 Si 2 (ν = 0.20) [42] is mainly reflected on ionic nature, and U 3 Si 2 C 2 shows better toughness compared to U 3 Si 2 , which may provide theoretical support for the application of U 3 Si 2 C 2 in nuclear reactors. As the pressure increases to 9.8 GPa, the atomic distance between U1 and U2 decreases from 3.249 Å to 3.098 Å, which may lead to the strength of the metal bond between uranium atoms. When the external pressure applied from 0 GPa to 8 GPa, the Young's modulus, shear modulus and bulk modulus of U 3 Si 2 C 2 increase 37.8, 34.5 and 14.2 GPa respectively, then decrease sharply as the external pressure exceeds 9.8 GPa, and the changing trend of which seems to be consistent with the variation of the crystal lattice constant. As shown in Figure 10, the three-dimensional Young's modulus distribution of U 3 Si 2 C 2 at different pressures and the general anisotropy index A U indicate that with the increase of hydrostatic pressure, the maximum Young's modulus increases gradually. The maximum Young's modulus of U 3 Si 2 C 2 under different pressures are all appear along the [110] direction. The maximum value of 3D Young's modulus is 252 GPa under 0 GPa pressure, and raising to 317 GPa, 289.26 GPa, and 319.59 GPa, respectively, under the external pressure 8 GPa, 12 GPa, and 20 GPa. Under high-pressure environment, the value of A U for U 3 Si 2 C 2 increases to higher than 0.8, while that remains basically unchanged when the pressure exceeds 8 GPa. Hence as the pressure increases, the degree of isotropy for U 3 Si 2 C 2 is reduced slightly.

Conclusions
In our present research, the structural characteristics, electronic behaviors as well as mechanical properties of U3Si2C2 compounds under high pressure have been calculated. The lattice constants and volumes decrease as the pressure increases; however, the lattice parameter c increased dramatically when the external stress reached 9.8 GPa, resulting in a sudden increase in the cell volume. U1 and U2 have opposite magnetic moment directions, while the magnetic moment direction of U1 atoms reverse under high pressures. Ionic and metallic nature coexist as to the bonding characteristics of U3Si2C2, and the ionic takes the dominant position in bonding. The temperature as well as the pressure influences the electronic thermal conductivity of U3Si2C2, while the electronic thermal conductivity remains basically unchanged with the applied pressure exceeds 10 GPa. Lastly, the maximum value of 3D Young's modulus is 252 GPa under 0 GPa, and the maximum Young's modulus of U3Si2C2 under different pressures are all appear along the [110] direction. Our theoretical investigation may provide support for the application of U-Si-Cbased fuels.

Conclusions
In our present research, the structural characteristics, electronic behaviors as well as mechanical properties of U 3 Si 2 C 2 compounds under high pressure have been calculated. The lattice constants and volumes decrease as the pressure increases; however, the lattice parameter c increased dramatically when the external stress reached 9.8 GPa, resulting in a sudden increase in the cell volume. U1 and U2 have opposite magnetic moment directions, while the magnetic moment direction of U1 atoms reverse under high pressures. Ionic and metallic nature coexist as to the bonding characteristics of U 3 Si 2 C 2 , and the ionic takes the dominant position in bonding. The temperature as well as the pressure influences the electronic thermal conductivity of U 3 Si 2 C 2 , while the electronic thermal conductivity remains basically unchanged with the applied pressure exceeds 10 GPa. Lastly, the maximum value of 3D Young's modulus is 252 GPa under 0 GPa, and the maximum Young's modulus of U 3 Si 2 C 2 under different pressures are all appear along the [110] direction. Our theoretical investigation may provide support for the application of U-Si-C-based fuels.