N-Calculator与NUFER耦合模型和组合预测法在食物氮足迹研究中的应用

为改善传统N-Calculator模型的局限性,加强食物氮足迹核算结果与环境影响的联系,提高活性氮管理预见性,将N-Calculator模型与食物系统养分流动(Nutrient Flows in Food Chains,Environment and Resources Use,NUFER)模型进行耦合,以估算我国2001—2020年人均食物氮足迹,建立组合预测体系。结果显示:2001—2020年,我国人均食物氮足迹由16.04 kg N/a增至18.95 kg N/a;全国食物氮足迹由20.47 Mt N/a增至26.76 Mt N/a;居民饮食结构正由以植物源食物为主的低氮消费模式转向以动物源食物为主的高氮消费模式;食物生产过程产生的活性氮的最终归宿为大气(64.3%)、水体和深层土壤(35.7%);我国食物氮足迹与人均可支配收入、城市化率、动物源食物消费氮占比呈正相关性,与恩格尔系数呈负相关性;未来10 a我国人均食物氮足迹呈增长趋势,预测结果显示年均增幅为0.16 kg N/a。

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.

Design of high-performance ion-doped CoP systems for hydrogen evolution:From multi-level screening calculations to experiment

Rational design of high-performance electrocatalysts for hydrogen evolution reaction(HER)is vital for future renewable energy systems.The incorporation of foreign metal ions into catalysts can be an effective approach to optimize its performance.However,there is a lack of systematic theoretical studies to reveal the quantitative relationships at the electronic level.Here,we develop a multi-level screening methodology to search for highly stable and active dopants for CoP catalysts.The density functional theory(DFT)calculations and symbolic regression(SR)were performed to investigate the relationship between the adsorption free energy(ΔG_(H^(*)))and 10 electronic parameters.The mathematic formulas derived from SR indicate that the difference of work function(ΔΦ)between doped metal and the acceptor plays the most important role in regulatingΔG_(H^(*)),followed by the d-band center(d-BC)of doped system.The descriptor of HER can be expressed asΔG_(H^(*))=1.59×√|0.188ΔΦ+d BC+0.120|1/2-0.166 with a high determination coefficient(R^(2)=0.807).Consistent with the theoretical prediction,experimental results show that the Al-CoP delivers superior electrocatalytic HER activity with a low overpotential of75 m V to drive a current density of 10 mA cm^(-2),while the overpotentials for undoped CoP,Mo-CoP,and V-CoP are 206,134,and 83 m V,respectively.The current work proves that theΔΦis the most significant regulatory parameter ofΔG_(H^(*))for ion-doped electrocatalysts.This finding can drive the discovery of high-performance ion-doped electrocatalysts,which is crucial for electrocatalytic water splitting.

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.

Passivation of PEA^(+)to CsPbI_(3)(110)Surface States:From the First Principles Calculations

This work investigates the effect of passivation on the electronic properties of inorganic perovskite CsPbI_(3)materials by using first-principles calculations with density functional theory(DFT).The passivation effect after the addition of Phenylethylamine(PEA+)molecule to CsPbI_(3)(110)surface is studied.The results of density of states(DOS)calculations show that the CsPbI_(3)(110)surface model with I atom terminated reveals new electronic DOS peaks(surface states)near the Fermi level.These surface states are mainly due to the contribution of I-5p orbital and are harmful to the CsPbI_(3)-based solar cells because they reduce the photoelectric conversion efficiency.The surface states near the Fermi level are significantly reduced,and the decline rate reaches 38.8%with the addition with PEA+molecule to the CsPbI3(110)surface.

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.

Direct atomic-level insight into oxygen reduction reaction on size-dependent Pt-based electrocatalysts from density functional theory calculations

Developing novel oxygen reduction reaction(ORR)catalysts with high activity is urgent for proton exchange membrane fuel cells.Herein,we investigated a group of size-dependent Pt-based catalysts as promising ORR catalysts by density functional theory calculations,ranging from single-atom,nanocluster to bulk Pt catalysts.The results showed that the ORR overpotential of these Pt-based catalysts increased when its size enlarged to the nanoparticle scale or reduced to the single-atom scale,and the Pt_(38)cluster had the lowest ORR overpotential(0.46 V)compared with that of Pt_(111)(0.57 V)and single atom Pt(0.7 V).Moreover,we established a volcano curve relationship between the ORR overpotential and binding energy of O*(ΔE_(O*),confirming the intermediate species anchored on Pt38cluster with suitable binding energy located at top of volcano curve.The interaction between intermediate species and Pt-based catalysts were also investigated by the charge distribution and projected density of state and which further confirmed the results of volcano curve.

Benchmarking calculations of excitation energies and transition properties with spectroscopic accuracy of highly charged ions used for the fusion plasma and astrophysical plasma

Atomic radiative data such as excitation energies, transition wavelengths, radiative rates, and level lifetimes with high precision are the essential parameters for the abundance analysis, simulation, and diagnostics in fusion and astrophysical plasmas. In this work, we mainly focus on reviewing our two projects performed in the past decade. One is about the ions with Z■30 that are generally of astrophysical interest, and the other one is about the highly charged krypton(Z = 36)and tungsten(Z = 74) ions that are relevant in research of magnetic confinement fusion. Two different and independent methods, namely, multiconfiguration Dirac–Hartree–Fock(MCDHF) and the relativistic many-body perturbation theory(RMBPT) are usually used in our studies. As a complement/extension to our previous works for highly charged tungsten ions with open M-shell and open N-shell, we also mainly focus on presenting and discussing our complete RMBPT and MCDHF calculations for the excitation energies, wavelengths, electric dipole(E1), magnetic dipole(M1), electric quadrupole(E2), and magnetic quadrupole(M2) transition properties, and level lifetimes for the lowest 148 levels belonging to the 3l3configurations in Al-like W61+. We also summarize the uncertainties of our systematical theoretical calculations, by cross-checking/validating our datasets from our RMBPT and MCDHF calculations, and by detailed comparisons with available accurate observations and other theoretical calculations. The data are openly available in Science Data Bank at https://doi.org/10.57760/sciencedb.10569.

Random Green's Function Method for Large-Scale Electronic Structure Calculation

We report a linear-scaling random Green's function(rGF) method for large-scale electronic structure calculation. In this method, the rGF is defined on a set of random states and is efficiently calculated by projecting onto Krylov subspace. With the rGF method, the Fermi–Dirac operator can be obtained directly, avoiding the polynomial expansion to Fermi–Dirac function. To demonstrate the applicability, we implement the rGF method with the density-functional tight-binding method. It is shown that the Krylov subspace can maintain at small size for materials with different gaps at zero temperature, including H_(2)O and Si clusters. We find with a simple deflation technique that the rGF self-consistent calculation of H_(2)O clusters at T = 0 K can reach an error of~ 1 me V per H_(2)O molecule in total energy, compared to deterministic calculations. The rGF method provides an effective stochastic method for large-scale electronic structure simulation.

Structural,electronic,and Li-ion mobility properties of garnet-type Li_(7)La_(3)Zr_(2)O_(12) surface:An insight from first-principles calculations

Garnet-type Li_(7)La_(3)Zr_(2)O_(12)(LLZO) is a promising solid-state electrolyte for Li-ion batteries,but Li-dendrite's formation greatly limits the applications.In this paper,we systematically investigate the stability,electronic properties,and Li-ion mobility of the LLZO surface by the ifrst-principles calculations.We consider the(110) and(001) slab structures with different terminations in the t-and c-LLZO.Our results indicate that both(110) and(001) surfaces prefer to form Li-rich termination due to their low surface energies for either t-or c-LLZO.Moreover,with the decrease of Li contents the stability of Li-rich surfaces is improved initially and degrades later.Unfortunately,the localized surface states at the Fermi level can induce the formation of metallic Li on the Li-rich surfaces.In comparison,Li/La-termination has a relatively low metallic Li formation tendency due to its rather low diffusion barrier.In fact,Li-ion can spontaneously migrate along path II(Li3→Li2) on the Li/La-T(001) surface.In contrast,it is more difficult for Li-ion diffusion on the Li-T(001) surface,which has a minimum diffusion barrier of 0.50 eV.Interestingly,the minimum diffusion barrier decreases to 0.34 eV when removing four Li-ions from the Li-T(001) surface.Thus,our study suggests that by varying Li contents,the stability and Li-ion diffusion barrier of LLZO surfaces can be altered favorably.These advantages can inhibit the formation of metallic Li on the LLZO surfaces.

Emerging perovskite materials for supercapacitors:Structure,synthesis,modification,advanced characterization,theoretical calculation and electrochemical performance

As a new generation electrode materials for energy storage,perovskites have attracted wide attention because of their unique crystal structure,reversible active sites,rich oxygen vacancies,and good stability.In this review,the design and engineering progress of perovskite materials for supercapacitors(SCs)in recent years is summarized.Specifically,the review will focus on four types of perovskites,perovskite oxides,halide perovskites,fluoride perovskites,and multi-perovskites,within the context of their intrinsic structure and corresponding electrochemical performance.A series of experimental variables,such as synthesis,crystal structure,and electrochemical reaction mechanism,will be carefully analyzed by combining various advanced characterization techniques and theoretical calculations.The applications of these materials as electrodes are then featured for various SCs.Finally,we look forward to the prospects and challenges of perovskite-type SCs electrodes,as well as the future research direction.

Organic Compounds Possessing the Plastic Crystalline Phase: Calculation of Their Fusion Enthalpies

For the first time, for different organic and inorganic compounds possessing the plastic crystalline phase, a new semiempirical equation describing dependence of their fusion enthalpies on such physico-chemical quantities as normal melting temperature, surface tension, molar volume and critical molar volume is received on the base of the principle of corresponding states and the energy equipartition theorem. Moreover, the proposed equation allows one to take into account the particularities of one-particle molecular rotation in the plastic crystalline phase.

Utility and Application of a Versatile Analytical Method for MMF Calculation in AC Machines

A versatile analytical method(VAM) for calculating the harmonic components of the magnetomotive force(MMF) generated by diverse armature windings in AC machines has been proposed, and the versatility of this method has been established in early literature. However, its practical applications and significance in advancing the analysis of AC machines need further elaboration. This paper aims to complement VAM by augmenting its theory, offering additional insights into its conclusions, as well as demonstrating its utility in assessing armature windings and its application of calculating torque for permanent magnet synchronous machines(PMSM). This work contributes to advancing the analysis of AC machines and underscores the potential for improved design and performance optimization.

High-throughput calculation-based rational design of Fe-doped MoS_(2) nanosheets for electrocatalytic p H-universal overall water splitting

Electrocatalytic water splitting is crucial for H2generation via hydrogen evolution reaction(HER)but subject to the sluggish dynamics of oxygen evolution reaction(OER).In this work,single Fe atomdoped MoS_(2)nanosheets(SFe-DMNs)were prepared based on the high-throughput density functional theory(DFT)calculation screening.Due to the synergistic effect between Fe atom and MoS_(2)and optimized intermediate binding energy,the SFe-DMNs could deliver outstanding activity for both HER and OER.When assembled into a two-electrode electrolytic cell,the SFe-DMNs could achieve the current density of 50 mA cm^(-2)at a low cell voltage of 1.55 V under neutral condition.These results not only confirmed the effectiveness of high-throughput screening,but also revealed the excellent activity and thus the potential applications in fuel cells of SFe-DMNs.

Syntheses,crystal structures,and quantum chemistry calculation of two Ni(Ⅱ)coordination polymers

Two new coordination polymers,[Ni(Hpdc)(bib)(H_(2)O)]_(n)(1)and{[Ni(bib)_(3)](ClO_(4))_(2)}_(n)(2),were prepared by mixing Ni^(2+),3,5⁃pyrazoledicarboxylic acid(H3pdc)/p⁃nitrobenzoic acid and 1,4⁃bis(imidazol⁃1⁃ylmethyl)butane(bib)by a hydrothermal method,respectively.X⁃ray crystallography reveals a 2D network constructed by six⁃coordinated Ni(Ⅱ)centers,bib,and Hpdc2-ligands in complex 1,while a 2D network is built by Ni(Ⅱ)and bib ligands in 2.Furthermore,the quantum⁃chemical calculations have been performed on‘molecular fragments’extracted from the crystal structure of 1 using the PBE0/LANL2DZ method in Gaussian 16 and the VASP program.CCDC:2343794,1;2343798,2.

Hemoglobin loss method calculates blood loss during pancreaticoduodenectomy and predicts bleeding-related risk factors

BACKGROUND The common clinical method to evaluate blood loss during pancreaticoduoden-ectomy(PD)is visual inspection,but most scholars believe that this method is extremely subjective and inaccurate.Currently,there is no accurate,objective me-thod to evaluate the amount of blood loss in PD patients.We retrospectively analyzed the clinical data of 341 patients who underwent PD in Shandong Provincial Hospital from March 2017 to February 2019.According to different surgical methods,they were divided into an open PD(OPD)group and a laparoscopic PD(LPD)group.The differences and correlations between the in-traoperative estimation of blood loss(IEBL)obtained by visual inspection and the intraoperative calculation of blood loss(ICBL)obtained using the Hb loss method were analyzed.ICBL,IEBL and perioperative calculation of blood loss(PCBL)were compared between the two groups,and single-factor regression analysis was performed.RESULTS There was no statistically significant difference in the preoperative general patient information between the two groups(P>0.05).PD had an ICBL of 743.2(393.0,1173.1)mL and an IEBL of 100.0(50.0,300.0)mL(P<0.001).There was also a certain correlation between the two(r=0.312,P<0.001).Single-factor analysis of ICBL showed that a history of diabetes[95%confidence interval(CI):53.82-549.62;P=0.017]was an independent risk factor for ICBL.In addition,the single-factor analysis of PCBL showed that body mass index(BMI)(95%CI:0.62-76.75;P=0.046)and preoperative total bilirubin>200μmol/L(95%CI:7.09-644.26;P=0.045)were independent risk factors for PCBL.The ICBLs of the LPD group and OPD group were 767.7(435.4,1249.0)mL and 663.8(347.7,1138.2)mL,respectively(P>0.05).The IEBL of the LPD group 200.0(50.0,200.0)mL was slightly greater than that of the OPD group 100.0(50.0,300.0)mL(P>0.05).PCBL was greater in the LPD group than the OPD group[1061.6(612.3,1632.3)mL vs 806.1(375.9,1347.6)mL](P<0.05).CONCLUSION The ICBL in patients who underwent PD was greater than the IEBL,but there is a certain correlation between the two.The Hb loss method can be used to evaluate intraoperative blood loss.A history of diabetes,preoperative bilirubin>200μmol/L and high BMI increase the patient's risk of bleeding.

Partitioning Calculation Method of Short-Circuit Current for High Proportion DG Access to Distribution Network

Aiming at the problemthat the traditional short-circuit current calculationmethod is not applicable to Distributed Generation(DG)accessing the distribution network,the paper proposes a short-circuit current partitioning calculation method considering the degree of voltage drop at the grid-connected point of DG.Firstly,the output characteristics of DG in the process of low voltage ride through are analyzed,and the equivalent output model of DG in the fault state is obtained.Secondly,by studying the network voltage distribution law after fault in distribution networks under different DG penetration rates,the degree of voltage drop at the grid-connected point of DG is used as a partition index to partition the distribution network.Then,iterative computation is performed within each partition,and data are transferred between partitions through split nodes to realize the fast partition calculation of short-circuit current for high proportion DG access to distribution network,which solves the problems of long iteration time and large calculation error of traditional short-circuit current.Finally,a 62-node real distribution network model containing a high proportion of DG access is constructed onMATLAB/Simulink,and the simulation verifies the effectiveness of the short-circuit current partitioning calculation method proposed in the paper,and its calculation speed is improved by 48.35%compared with the global iteration method.

Temperature field calculation of rail flash welding

The forging stage of rail flash welding has a decisive influence on joint strength,and the study of the temperature distribution in the process has an important role in further improving joint strength.In this paper,three calculation methods for the temperature field are given.First,the finite element model of the temperature field before forging rail flash welding is established by using the transient heat module of Ansys software and verified by infrared temperature measurement.Second,the temperature distribution of different parts of the rail before flash welding is obtained by using infrared thermal imaging equipment.Third,Matlab software is used to calculate the temperature of the non-measured part.Finally,the temperature distribution function along the rail axis is fitted through the temperature measurement data.The temperature distribution before the top forging of the rail flash welding can be used to analyze the joint and heat-affected zone organization and properties effectively and to guide the parameter setting and industrial production.

Effect of Types and Orders of Electromagnetic Field Finite Element Meshes on Power Communication Harmonic Parameters Calculation Results of Tubular Hydrogenerators

In generator design field,waveform total harmonic distortion(THD)and telephone harmonic factor(THF)are parameters commonly used to measure the impact of generator no-load voltage harmonics on the power communication quality.Tubular hydrogenerators are considered the optimal generator for exploiting low-head,high-flow hydro resources,and they have seen increasingly widespread application in China's power systems recent years.However,owing to the compact and constrained internal space of such generators,their internal magnetic-field harmonics are pronounced.Therefore,accurate calculation of their THD and THF is crucial during the analysis and design stages to ensure the quality of power communication.Especially in the electromagnetic field finite element modeling analysis of such generators,the type and order of the finite element meshes may have a significant impact on the THD and THF calculation results,which warrants in-depth research.To address this,this study takes a real 34 MW large tubular hydrogenerator as an example,and establishes its electromagnetic field finite element model under no-load conditions.Two types of meshes,five mesh densities,and two mesh orders are analyzed to reveal the effect of electromagnetic field finite element mesh types and orders on the calculation results of THD and THF for such generators.

First-Principles Calculation of the Topological Nodal-Line Semimetal FeGe_(2)

The electronic and topological properties of FeGe2 with a tetragonal crystal structure were investigated via first-principles calculations.The results demonstrate that FeGe2 in this structure exhibits anti-ferromagnetism,with two bands crossing the Fermi level nesting each other at high-symmetry points in the Brillouin zone,forming a nodal ring where the nodes intersect in momentum space.Additionally,it possesses nontrivial topological surface states.Upon inclusion of SOC(spin-orbit coupling),there are no significant changes observed in the band structure,nodal features,or surface states,indicating the persistence of its topological nodal-line characteristics.