An accurate and simultaneous ab initio prediction for both light nuclei and nuclear matter has been a longstanding challenge in nuclear physics, due to the significant uncertainties associated with the three-nucleon f...An accurate and simultaneous ab initio prediction for both light nuclei and nuclear matter has been a longstanding challenge in nuclear physics, due to the significant uncertainties associated with the three-nucleon forces.In this Letter, we develop the relativistic quantum Monte Carlo methods for the nuclear ab initio problem, and calculate the ground-state energies of A ≤ 4 nuclei using the two-nucleon Bonn force with an unprecedented high accuracy. The present relativistic results significantly outperform the nonrelativistic results with only twonucleon forces. We demonstrate that both light nuclei and nuclear matter can be well described simultaneously in the relativistic ab initio calculations, even in the absence of three-nucleon forces, and a correlation between the properties of light A ≤ 4 nuclei and the nuclear saturation is revealed. This provides a quantitative understanding of the connection between the light nuclei and nuclear matter saturation properties.展开更多
Understanding the adsorption behavior of heavy metals and metalloids on clay minerals is essential for remediating heavy metal-contaminated soils.The adsorption of heavy metals and metalloids on illite(001)and sodium ...Understanding the adsorption behavior of heavy metals and metalloids on clay minerals is essential for remediating heavy metal-contaminated soils.The adsorption of heavy metals and metalloids on illite(001)and sodium montmorillonite(Na-MMT)(001)surfaces was investigated using first-principles calculations in this study,especially As atom and H_(3)AsO_(3) molecule.The adsorption energies of the As atom were−1.94 eV on the illite(001)and−0.56 eV on the Na-MMT(001),whereas,the adsorption energies of the H_(3)AsO_(3) molecule were−1.40 eV on illite(001)and−1.01 eV on Na-MMT(001).The above results indicate that the adsorption was more energetically favorable on illite(001).Additionally,compared to Na-MMT(001),there were more significant interactions between the atoms/molecules on the illite(001).After As atom and H_(3)AsO_(3) molecule adsorption,the electrons were transferred from mineral surface atoms to the adsorbates on both illite(001)and Na-MMT(001)surfaces.Moreover,the adsorption of As atom on illite(001)and Na-MMT(001)surfaces were more energy favorable compared to Hg,Cd,and Cr atoms.Overall,this work provides new insights into the adsorption behavior of As atoms and As molecules on illite and Na-MMT.The results indicate that illite rich soils are more prone to contamination by arsenic compared to soils primarily composed of Na-MMT minerals.展开更多
The electrical conductivity of minerals under extreme conditions is governed by their variations in composition and structure.Constitution water,which is present in various polymorphic phases of olivine,can significan...The electrical conductivity of minerals under extreme conditions is governed by their variations in composition and structure.Constitution water,which is present in various polymorphic phases of olivine,can significantly enhance electrical conductivity under mantle pressure-temperature conditions,therefore playing a key role in proton transport.Despite this,the conductive mechanisms in hydrous olivine,particularly in hydrous ringwoodite,and the dynamic behavior of hydrogen at elevated temperatures,remain poorly understood.In this study,we investigated the proton conduction mechanisms in hydrous ringwoodite through first-principles calculations.Several hydrous models were considered,and ab initio molecular dynamics(AIMD)simulations were employed to simulate hydrous configurations at high temperatures.Calculations based on density functional perturbation theory(DFPT)and vibrational density of states(VDOS)analyses were conducted to probe the stability of hydrous structures,and investigate the dynamic behavior of internal hydrogen.Our results indicate that hydrogen trapped in Mg^(2+)and Fe^(3+)defects exhibits significantly higher mobility than hydrogen trapped in Si^(4+)defects.At elevated temperatures,we observed the ionization of hydrogen from cationic defects,leading to high and highly anisotropic proton conductivity along the[100]crystallographic direction.This thermal ionization-induced anisotropic conductivity is consistent with experimental observations of olivine single crystals.Finally,the conductivity of the 0.79 wt%hydrous ringwoodite structure was found to range from 10^(-0.3)to 10^(0.4)S/m,the 1.19 wt%structure ranged from 10^(0.4)to 10^(0.9)S/m in the transition region,and the 1.62 wt%structure exhibited conductivity ranging from 10^(0.7)to 10^(1.2)S/m.These results are in excellent agreement with prior experimental data,providing further insight into the proton conduction mechanisms of hydrous olivine under extreme mantle conditions.展开更多
Hydrated ions play essential roles in diverse chemical and biological processes,yet accurately characterizing their hydration structures remains challenging due to the delicate interplay of ion–water and water–water...Hydrated ions play essential roles in diverse chemical and biological processes,yet accurately characterizing their hydration structures remains challenging due to the delicate interplay of ion–water and water–water interactions.Here,we use ab initio molecular dynamics(AIMD)simulations based on the strongly constrained and appropriately normed(SCAN)exchange–correlation functional to systematically investigate the hydration structures of eight representative ions(Mg^(2+),Ca^(2+),Li^(+),Na^(+),K^(+),F^(-),Cl^(-),Br^(-))in aqueous solution.Compared to the widely used Perdew–Burke–Ernzerhof(PBE)functional,SCAN substantially improves the description of solvent water by weakening the hydrogen-bond network and enhancing structural disorder,yielding results in closer agreement with experiments.SCAN modifies ionic hydration shells in an ion-specific manner,governed by ionic size and charge,and reproduces experimental hydration geometries especially well for intermediate-size monovalent ions(Na^(+),Cl^(-)).Moreover,SCAN consistently reduces the overpolarization of water molecules near ions.These improvements lead to more accurate and physically consistent hydration structures,highlighting SCAN's utility for modeling complex aqueous systems and offering guidance for future studies of ionic solvation.展开更多
High multipole electromagnetic transitions are rare in nature.The highest-multipole transition observed in atomic nuclei is the electric hexacontatetrapole E6 transition from the T_(1/2)=2.54(2)-min J^(π)=1_(9/2)-iso...High multipole electromagnetic transitions are rare in nature.The highest-multipole transition observed in atomic nuclei is the electric hexacontatetrapole E6 transition from the T_(1/2)=2.54(2)-min J^(π)=1_(9/2)-isomer to the 7/2^(-)ground state in^(53)Fe with an angular momentum change of six units.In the present work,we performed ab initio calculations for this unique case by employing chiral effective field theory(EFT)forces.The in-medium similarity renormalization group is used to derive the valence-space effective Hamiltonian and multipolar transition operators.Bare nucleon charges were used in all the multipolar transition rate calculations,providing good agreement with the experimental data.The valence space takes the full fp shell.In^(53)Fe,the low-lying states were dominated by the 0f_(7/2)component.Two different versions of the chiral EFT two-plus three-nucleon interaction were used to test the dependence on the interaction used.We also tested the convergence of the transition rate calculations against the harmonic oscillator parameter hΩand basis truncations e_(max)and E_(3max)for twoand three-nucleon forces,respectively.展开更多
Thermodynamic optimization of the AF-BeF_(2)(A=K,Rb,and Cs),KF-CsF,and RbF-CsF systems was performed within the framework of phase diagrams calculation.The model parameters were optimized based on experimental data an...Thermodynamic optimization of the AF-BeF_(2)(A=K,Rb,and Cs),KF-CsF,and RbF-CsF systems was performed within the framework of phase diagrams calculation.The model parameters were optimized based on experimental data and theoretically calculated values.The results show that the thermodynamically calculated values for the AF-BeF_(2)(A=K,Rb,and Cs),KF-CsF,and RbF-CsF systems agree well with the experimental data.Next,a set of reliable and self-consistent thermodynamic databases was built,and the liquidus projections and invariant points of the sub-ternary systems of the KF-RbF-CsF-BeF_(2)system were calculated.Furthermore,the melting temperature with the corresponding composition was predicted using the phase diagrams calculation technique,and the radial distribution functions,coordination numbers,angular distribution functions,and diffusion coefficients of the quaternary KF-RbF-CsF-BeF_(2)system were calculated using ab initio molecular dynamics.The results show that the quaternary KF-RbF-CsF-BeF_(2)system with the proportion 3.50-28.92-21.78-45.80 mol%or 1.80-35.42-52.40-10.38 mol%is one of the most promising candidate coolants for molten salt reactors in terms of thermodynamics and kinetics.This work provides direct guidelines for the screening and optimization of molten salts in the nuclear energy field.展开更多
As the size of transistors shrinks and power density increases,thermal simulation has become an indispensable part of the device design procedure.However,existing works for advanced technology transistors use simplifi...As the size of transistors shrinks and power density increases,thermal simulation has become an indispensable part of the device design procedure.However,existing works for advanced technology transistors use simplified empirical models to calculate effective thermal conductivity in the simulations.In this work,we present a dataset of size-dependent effective thermal conductivity with electron and phonon properties extracted from ab initio computations.Absolute in-plane and cross-plane thermal conductivity data of eight semiconducting materials(Si,Ge,GaN,AlN,4H-SiC,GaAs,InAs,BAs)and four metallic materials(Al,W,TiN,Ti)with the characteristic length ranging from 5 nm to 50 nm have been provided.Besides the absolute value,normalized effective thermal conductivity is also given,in case it needs to be used with updated bulk thermal conductivity in the future.展开更多
We calculate the electrical and thermal conductivity of hydrogen for a wide range of densities and temperatures by using molecular dynamics simulations informed by density functional theory.On the basis of the corresp...We calculate the electrical and thermal conductivity of hydrogen for a wide range of densities and temperatures by using molecular dynamics simulations informed by density functional theory.On the basis of the corresponding extended ab initio data set,we construct interpolation formulas covering the range from low-density,high-temperature to high-density,low-temperature plasmas.Our conductivity model repro-duces the well-known limits of the Spitzer and Ziman theory.We compare with available experimental data andfind very good agreement.The new conductivity model can be applied,for example,in dynamo simulations for magneticfield generation in gas giant planets,brown dwarfs,and stellar envelopes.展开更多
The energy bands,electronic structures of CuN3 and AgN3 crystallines were investigated by periodic ab initio method.The charge density projection shows that there are overlaps of isodensities between the terminal nitr...The energy bands,electronic structures of CuN3 and AgN3 crystallines were investigated by periodic ab initio method.The charge density projection shows that there are overlaps of isodensities between the terminal nitrogen and metallic ion,indicating that the metals and the azides are combined by covalent bonds.The crystal lattice energies are-781.05 and-840.83 kJ/mol for CuN3 and AgN3 respectively.These results approach the data obtained by Gray′s approximate method.The frontier crystal orbital mainly consists of the atomic orbital of azide′s terminal nitrogen.The energy gap for AgN3 is smaller than that of CuN3,and the highest occupied crystal orbitals of AgN-3 consist of both the atomic orbitals of the terminal nitrogen in azide and the silver ion,which facilitates the electron to leap from terminal nitrogen in azide to metallic ion directly.Hence silver azide is slightly more sensitive than copper azide.The elastic coefficients C11,C22 and C33 of CuN3 are predicted to be 96.52,96.86 and 154.06 GPa,C11 and C22 of AgN3 are 303.29 and 138.80 GPa.展开更多
基金supported in part by the National Natural Science Foundation of China (Grant Nos. 12141501, 123B2080, 12435006, 12475117, and 11935003)the National Key Laboratory of Neutron Science and Technology (Grant No. NST202401016)+2 种基金the National Key R&D Program of China (Grant No. 2024YFE0109803)the High-performance Computing Platform of Peking Universitythe funding support from the State Key Laboratory of Nuclear Physics and Technology, Peking University (Grant No. NPT2023ZX03)。
文摘An accurate and simultaneous ab initio prediction for both light nuclei and nuclear matter has been a longstanding challenge in nuclear physics, due to the significant uncertainties associated with the three-nucleon forces.In this Letter, we develop the relativistic quantum Monte Carlo methods for the nuclear ab initio problem, and calculate the ground-state energies of A ≤ 4 nuclei using the two-nucleon Bonn force with an unprecedented high accuracy. The present relativistic results significantly outperform the nonrelativistic results with only twonucleon forces. We demonstrate that both light nuclei and nuclear matter can be well described simultaneously in the relativistic ab initio calculations, even in the absence of three-nucleon forces, and a correlation between the properties of light A ≤ 4 nuclei and the nuclear saturation is revealed. This provides a quantitative understanding of the connection between the light nuclei and nuclear matter saturation properties.
基金Project(22376221)supported by the National Natural Science Foundation of ChinaProject(2024JJ2074)supported by the Natural Science Foundation of Hunan Province,ChinaProject(2023QNRC001)supported by the Young Elite Scientists Sponsorship Program by CAST。
文摘Understanding the adsorption behavior of heavy metals and metalloids on clay minerals is essential for remediating heavy metal-contaminated soils.The adsorption of heavy metals and metalloids on illite(001)and sodium montmorillonite(Na-MMT)(001)surfaces was investigated using first-principles calculations in this study,especially As atom and H_(3)AsO_(3) molecule.The adsorption energies of the As atom were−1.94 eV on the illite(001)and−0.56 eV on the Na-MMT(001),whereas,the adsorption energies of the H_(3)AsO_(3) molecule were−1.40 eV on illite(001)and−1.01 eV on Na-MMT(001).The above results indicate that the adsorption was more energetically favorable on illite(001).Additionally,compared to Na-MMT(001),there were more significant interactions between the atoms/molecules on the illite(001).After As atom and H_(3)AsO_(3) molecule adsorption,the electrons were transferred from mineral surface atoms to the adsorbates on both illite(001)and Na-MMT(001)surfaces.Moreover,the adsorption of As atom on illite(001)and Na-MMT(001)surfaces were more energy favorable compared to Hg,Cd,and Cr atoms.Overall,this work provides new insights into the adsorption behavior of As atoms and As molecules on illite and Na-MMT.The results indicate that illite rich soils are more prone to contamination by arsenic compared to soils primarily composed of Na-MMT minerals.
基金supported by the National Natural Science Foundation of China(41930112,91755215)the Sichuan Students'Innovation and entrepreneurship training program(s202310616086).
文摘The electrical conductivity of minerals under extreme conditions is governed by their variations in composition and structure.Constitution water,which is present in various polymorphic phases of olivine,can significantly enhance electrical conductivity under mantle pressure-temperature conditions,therefore playing a key role in proton transport.Despite this,the conductive mechanisms in hydrous olivine,particularly in hydrous ringwoodite,and the dynamic behavior of hydrogen at elevated temperatures,remain poorly understood.In this study,we investigated the proton conduction mechanisms in hydrous ringwoodite through first-principles calculations.Several hydrous models were considered,and ab initio molecular dynamics(AIMD)simulations were employed to simulate hydrous configurations at high temperatures.Calculations based on density functional perturbation theory(DFPT)and vibrational density of states(VDOS)analyses were conducted to probe the stability of hydrous structures,and investigate the dynamic behavior of internal hydrogen.Our results indicate that hydrogen trapped in Mg^(2+)and Fe^(3+)defects exhibits significantly higher mobility than hydrogen trapped in Si^(4+)defects.At elevated temperatures,we observed the ionization of hydrogen from cationic defects,leading to high and highly anisotropic proton conductivity along the[100]crystallographic direction.This thermal ionization-induced anisotropic conductivity is consistent with experimental observations of olivine single crystals.Finally,the conductivity of the 0.79 wt%hydrous ringwoodite structure was found to range from 10^(-0.3)to 10^(0.4)S/m,the 1.19 wt%structure ranged from 10^(0.4)to 10^(0.9)S/m in the transition region,and the 1.62 wt%structure exhibited conductivity ranging from 10^(0.7)to 10^(1.2)S/m.These results are in excellent agreement with prior experimental data,providing further insight into the proton conduction mechanisms of hydrous olivine under extreme mantle conditions.
基金supported by the National Natural Science Foundation of China(Grant Nos.12535001,11935002,and 11525520)the National Key Research and Development Program of China(Grant No.2021YFA1400500)。
文摘Hydrated ions play essential roles in diverse chemical and biological processes,yet accurately characterizing their hydration structures remains challenging due to the delicate interplay of ion–water and water–water interactions.Here,we use ab initio molecular dynamics(AIMD)simulations based on the strongly constrained and appropriately normed(SCAN)exchange–correlation functional to systematically investigate the hydration structures of eight representative ions(Mg^(2+),Ca^(2+),Li^(+),Na^(+),K^(+),F^(-),Cl^(-),Br^(-))in aqueous solution.Compared to the widely used Perdew–Burke–Ernzerhof(PBE)functional,SCAN substantially improves the description of solvent water by weakening the hydrogen-bond network and enhancing structural disorder,yielding results in closer agreement with experiments.SCAN modifies ionic hydration shells in an ion-specific manner,governed by ionic size and charge,and reproduces experimental hydration geometries especially well for intermediate-size monovalent ions(Na^(+),Cl^(-)).Moreover,SCAN consistently reduces the overpolarization of water molecules near ions.These improvements lead to more accurate and physically consistent hydration structures,highlighting SCAN's utility for modeling complex aqueous systems and offering guidance for future studies of ionic solvation.
基金supported by the National Key R&D Program of China(Nos.2024YFA1610900 and 2023YFA1606401)the National Natural Science Foundation of China(Nos.12335007 and 12035001)the United Kingdom Science and Technology Facilities Council(No.ST/V001108/1)。
文摘High multipole electromagnetic transitions are rare in nature.The highest-multipole transition observed in atomic nuclei is the electric hexacontatetrapole E6 transition from the T_(1/2)=2.54(2)-min J^(π)=1_(9/2)-isomer to the 7/2^(-)ground state in^(53)Fe with an angular momentum change of six units.In the present work,we performed ab initio calculations for this unique case by employing chiral effective field theory(EFT)forces.The in-medium similarity renormalization group is used to derive the valence-space effective Hamiltonian and multipolar transition operators.Bare nucleon charges were used in all the multipolar transition rate calculations,providing good agreement with the experimental data.The valence space takes the full fp shell.In^(53)Fe,the low-lying states were dominated by the 0f_(7/2)component.Two different versions of the chiral EFT two-plus three-nucleon interaction were used to test the dependence on the interaction used.We also tested the convergence of the transition rate calculations against the harmonic oscillator parameter hΩand basis truncations e_(max)and E_(3max)for twoand three-nucleon forces,respectively.
基金supported by the National Natural Science Foundation of China(Nos.12205364 and 12375282)Guangdong Provincial Natural Science Foundation(Nos.2024A1515012570 and 2024A1515010885)the Fundamental Research funds for the Central Universities,Sun Yat sen University。
文摘Thermodynamic optimization of the AF-BeF_(2)(A=K,Rb,and Cs),KF-CsF,and RbF-CsF systems was performed within the framework of phase diagrams calculation.The model parameters were optimized based on experimental data and theoretically calculated values.The results show that the thermodynamically calculated values for the AF-BeF_(2)(A=K,Rb,and Cs),KF-CsF,and RbF-CsF systems agree well with the experimental data.Next,a set of reliable and self-consistent thermodynamic databases was built,and the liquidus projections and invariant points of the sub-ternary systems of the KF-RbF-CsF-BeF_(2)system were calculated.Furthermore,the melting temperature with the corresponding composition was predicted using the phase diagrams calculation technique,and the radial distribution functions,coordination numbers,angular distribution functions,and diffusion coefficients of the quaternary KF-RbF-CsF-BeF_(2)system were calculated using ab initio molecular dynamics.The results show that the quaternary KF-RbF-CsF-BeF_(2)system with the proportion 3.50-28.92-21.78-45.80 mol%or 1.80-35.42-52.40-10.38 mol%is one of the most promising candidate coolants for molten salt reactors in terms of thermodynamics and kinetics.This work provides direct guidelines for the screening and optimization of molten salts in the nuclear energy field.
基金Project supported by the National Key R&D Project from Ministry of Science and Technology of China(Grant No.2022YFA1203100)the National Natural Science Foundation of China(Grant No.52122606)the funding from Shanghai Polytechnic University.
文摘As the size of transistors shrinks and power density increases,thermal simulation has become an indispensable part of the device design procedure.However,existing works for advanced technology transistors use simplified empirical models to calculate effective thermal conductivity in the simulations.In this work,we present a dataset of size-dependent effective thermal conductivity with electron and phonon properties extracted from ab initio computations.Absolute in-plane and cross-plane thermal conductivity data of eight semiconducting materials(Si,Ge,GaN,AlN,4H-SiC,GaAs,InAs,BAs)and four metallic materials(Al,W,TiN,Ti)with the characteristic length ranging from 5 nm to 50 nm have been provided.Besides the absolute value,normalized effective thermal conductivity is also given,in case it needs to be used with updated bulk thermal conductivity in the future.
基金supported by the Priority Program SPP 1992 of the German Science Foundation(DFG)The Diversity of Exoplanets under project number 362460292.
文摘We calculate the electrical and thermal conductivity of hydrogen for a wide range of densities and temperatures by using molecular dynamics simulations informed by density functional theory.On the basis of the corresponding extended ab initio data set,we construct interpolation formulas covering the range from low-density,high-temperature to high-density,low-temperature plasmas.Our conductivity model repro-duces the well-known limits of the Spitzer and Ziman theory.We compare with available experimental data andfind very good agreement.The new conductivity model can be applied,for example,in dynamo simulations for magneticfield generation in gas giant planets,brown dwarfs,and stellar envelopes.
文摘The energy bands,electronic structures of CuN3 and AgN3 crystallines were investigated by periodic ab initio method.The charge density projection shows that there are overlaps of isodensities between the terminal nitrogen and metallic ion,indicating that the metals and the azides are combined by covalent bonds.The crystal lattice energies are-781.05 and-840.83 kJ/mol for CuN3 and AgN3 respectively.These results approach the data obtained by Gray′s approximate method.The frontier crystal orbital mainly consists of the atomic orbital of azide′s terminal nitrogen.The energy gap for AgN3 is smaller than that of CuN3,and the highest occupied crystal orbitals of AgN-3 consist of both the atomic orbitals of the terminal nitrogen in azide and the silver ion,which facilitates the electron to leap from terminal nitrogen in azide to metallic ion directly.Hence silver azide is slightly more sensitive than copper azide.The elastic coefficients C11,C22 and C33 of CuN3 are predicted to be 96.52,96.86 and 154.06 GPa,C11 and C22 of AgN3 are 303.29 and 138.80 GPa.
