The Probability Density Function Related to Shallow Cumulus Entrainment Rate and Its Influencing Factors in a Large-Eddy Simulation

The process of entrainment-mixing between cumulus clouds and the ambient air is important for the development of cumulus clouds.Accurately obtaining the entrainment rate(λ)is particularly important for its parameterization within the overall cumulus parameterization scheme.In this study,an improved bulk-plume method is proposed by solving the equations of two conserved variables simultaneously to calculateλof cumulus clouds in a large-eddy simulation.The results demonstrate that the improved bulk-plume method is more reliable than the traditional bulk-plume method,becauseλ,as calculated from the improved method,falls within the range ofλvalues obtained from the traditional method using different conserved variables.The probability density functions ofλfor all data,different times,and different heights can be well-fitted by a log-normal distribution,which supports the assumed stochastic entrainment process in previous studies.Further analysis demonstrate that the relationship betweenλand the vertical velocity is better than other thermodynamic/dynamical properties;thus,the vertical velocity is recommended as the primary influencing factor for the parameterization ofλin the future.The results of this study enhance the theoretical understanding ofλand its influencing factors and shed new light on the development ofλparameterization.

Sensitivity impacts owing to the variations in the type of zero-range pairing forces on the fission properties using the density functional theory

Using the Skyrme density functional theory,potential energy surfaces of^(240)Pu with constraints on the axial quadrupole and octupole deformations(q_(20)and q_(30))were calculated.The volume-like and surface-like pairing forces,as well as a combination of these two forces,were used for the Hartree–Fock–Bogoliubov approximation.Variations in the least-energy fission path,fission barrier,pairing energy,total kinetic energy,scission line,and mass distribution of the fission fragments based on the different forms of the pairing forces were analyzed and discussed.The fission dynamics were studied based on the timedependent generator coordinate method plus the Gaussian overlap approximation.The results demonstrated a sensitivity of the mass and charge distributions of the fission fragments on the form of the pairing force.Based on the investigation of the neutron-induced fission of^(239)Pu,among the volume,mixed,and surface pairing forces,the mixed pairing force presented a good reproduction of the experimental data.

Density functional theory study of B- and Si-doped carbons and their adsorption interactions with sulfur compounds

Understanding the adsorption interactions between carbon materials and sulfur compounds has far-reaching impacts,in addition to their well-known important role in energy storage and conversion,such as lithium-ion batteries.In this paper,properties of intrinsic B or Si single-atom doped,and B-Si codoped graphene(GR)and graphdiyne(GDY)were investigated by using density functional theory-based calculations,in which the optimal doping configurations were explored for potential applications in adsorbing sulfur compounds.Results showed that both B or Si single-atom doping and B-Si codoping could substantially enhance the electron transport properties of GR and GDY,improving their surface activity.Notably,B and Si atoms displayed synergistic effects for the codoped configurations,where B-Si codoped GR/GDY exhibited much better performance in the adsorption of sulfurcontaining chemicals than single-atom doped systems.In addition,results demonstrated that,after B-Si codoping,the adsorption energy and charge transfer amounts of GDY with sulfur compounds were much larger than those of GR,indicating that B-Si codoped GDY might be a favorable material for more effectively interacting with sulfur reagents.

Local thermal conductivity of inhomogeneous nano-fluidic films:A density functional theory perspective

Combining the mean field Pozhar-Gubbins(PG)theory and the weighted density approximation,a novel method for local thermal conductivity of inhomogeneous fluids is proposed.The correlation effect that is beyond the mean field treatment is taken into account by the simulation-based empirical correlations.The application of this method to confined argon in slit pore shows that its prediction agrees well with the simulation results,and that it performs better than the original PG theory as well as the local averaged density model(LADM).In its further application to the nano-fluidic films,the influences of fluid parameters and pore parameters on the thermal conductivity are calculated and investigated.It is found that both the local thermal conductivity and the overall thermal conductivity can be significantly modulated by these parameters.Specifically,in the supercritical states,the thermal conductivity of the confined fluid shows positive correlation to the bulk density as well as the temperature.However,when the bulk density is small,the thermal conductivity exhibits a decrease-increase transition as the temperature is increased.This is also the case in which the temperature is low.In fact,the decrease-increase transition in both the small-bulk-density and low-temperature cases arises from the capillary condensation in the pore.Furthermore,smaller pore width and/or stronger adsorption potential can raise the critical temperature for condensation,and then are beneficial to the enhancement of the thermal conductivity.These modulation behaviors of the local thermal conductivity lead immediately to the significant difference of the overall thermal conductivity in different phase regions.

Solvent effects on Diels-Alder reaction in ionic liquids:A reaction density functional study

Extensive experimental studies have been performed on the Diels-Alder(DA)reactions in ionic liquids(ILs),which demonstrate that the IL environment can significantly influence the reaction rates and selectivity.However,the underlying microscopic mechanism remains ambiguous.In this work,the multiscale reaction density functional theory is applied to explore the effect of 1-butyl-3-methylimidazolium hexafluorophosphate([BMIM][PF_(6)])solvent on the reaction of cyclopentadiene(CP)with acrolein,methyl acrylate,or acrylonitrile.By analyzing the free energy landscape during the reaction,it is found that the polarization effect has a relatively small influence,while the solvation effect makes both the activation free energy and reaction free energy decrease.In addition,the rearrangement of local solvent structure shows that the cation spatial distribution responds more evidently to the reaction than the anion,and this indicates that the cation plays a dominant role in the solvation effect and so as to affect the reaction rates and selectivity of the DA reactions.

Probability Density Analysis of Nonlinear Random Ship Rolling

Ship rolling in random waves is a complicated nonlinear motion that contributes substantially to ship instability and capsizing.The finite element method(FEM)is employed in this paper to solve the Fokker Planck(FP)equations numerically for homoclinic and heteroclinic ship rolling under random waves described as periodic and Gaussian white noise excitations.The transient joint probability density functions(PDFs)and marginal PDFs of the rolling responses are also obtained.The effects of stimulation strength on ship rolling are further investigated from a probabilistic standpoint.The homoclinic ship rolling has two rolling states,the connection between the two peaks of the PDF is observed when the periodic excitation amplitude or the noise intensity is large,and the PDF is remarkably distributed in phase space.These phenomena increase the possibility of a random jump in ship motion states and the uncertainty of ship rolling,and the ship may lose stability due to unforeseeable facts or conditions.Meanwhile,only one rolling state is observed when the ship is in heteroclinic rolling.As the periodic excitation amplitude grows,the PDF concentration increases and drifts away from the beginning location,suggesting that the ship rolling substantially changes in a cycle and its stability is low.The PDF becomes increasingly uniform and covers a large region as the noise intensity increases,reducing the certainty of ship rolling and navigation safety.The current numerical solutions and analyses may be applied to evaluate the stability of a rolling ship in irregular waves and capsize mechanisms.

Effects of zinc on χ-Fe_(5)C_(2) for carbon dioxide hydrogenation to olefins:Insights from experimental and density function theory calculations

Production of light olefins from CO_(2), the primary greenhouse gases, is of great importance to mitigate the adverse effects of CO_(2) emission on environment and to supply the value-added products from nonpetroleum resource. However, development of robust catalyst with controllable selectivity and stability remains a challenge. Herein, we report that Zn-promoted Fe catalyst can boost the stable and selective production of light olefins from CO_(2). Specifically, the Zn-promoted Fe exhibits a highly stable activity and olefin selectivity over 200 h time-on-stream compared to the unpromoted Fe catalyst, primarily owing to the preservation of active χ-Fe_(5)C_(2) phase. Structural characterizations of the spent catalysts suggest that Zn substantially regulates the content of iron carbide on the surface and suppresses the reoxidation of bulk iron carbide during the reaction. DFT calculations confirm that adsorption of surface carbon atoms and graphene-like carbonaceous species are not thermochemically favored on Zn-promoted Fe catalyst. Carbon deposition by CAC coupling reactions of two surface carbon atoms and dehydrogenation of CH intermediate are also inhibited. Furthermore, the effects of Zn on antioxidation of iron carbide were also investigated. Zn favored the hydrogenation of surface adsorbed oxygen atoms to H_(2)O and the desorption of H_(2)O, which reduces the possibility of surface carbide being oxidized by the chemisorbed oxygen.

Investigation of the structural, electronic and mechanical properties of CaO–SiO_(2) compound particles in steel based on density functional theory

CaO–SiO_(2)compounds compromise one of the most common series of oxide particles in liquid steels, which could significantly affect the service performance of the steels as crack initiation sites. However, the structural, electronic, and mechanical properties of the compounds in CaO–SiO_(2)system are still not fully clarified due to the difficulties in the experiments. In this study, a thorough investigation of these properties of CaO–SiO_(2)compound particles in steels was conducted based on first-principles density functional theory. Corresponding phases were determined by thermodynamic calculation, including gamma dicalcium silicate(γ-C2S), alpha-prime(L) dicalcium silicate(αL′-C2S), alpha-prime(H) dicalcium silicate(αH′-C2S), alpha dicalcium silicate(α-C2S), rankinite(C3S2), hatrurite(C3S), wollastonite(CS), and pseudowollastonite(Ps-CS). The results showed that the calculated crystal structures of the eight phases agree well with the experimental results. All the eight phases are stable according to the calculated formation energies, and γ-C2S is the most stable. O atom contributes the most to the reactivity of these phases. The Young’s modulus of the eight phases is in the range of 100.63–132.04 GPa. Poisson’s ratio is in the range of0.249–0.281. This study provided further understanding concerning the CaO–SiO_(2)compound particles in steels and fulfilled the corresponding property database, paving the way for inclusion engineering and design in terms of fracture-resistant steels.

Density Functional Theory Study of Marine Polybrominated Diphenyl Ethers in Anaerobic Degradation

Polybrominated diphenyl ethers(PBDEs)are a kind of serious pollutants in the ocean.Biodegradation is considered as an economical and safe way for PBDEs removal and reductive debromination dominates the initial pathway of anaerobic degradation.On the basis of experimental study,Octa-BDE 197,Hepta-BDE 183,Hexa-BDE 153,Penta-BDE 99 and Tetra-BDE 47 were selected as the initial degradation objects,and their debromination degradation were studied using density functional theory.The structures were optimized by Gaussian 09 program.Furthermore,the molecular orbitals and charge distribution were analyzed.All C-Br bond dissociation energies at different positions including ortho,meta and para bromine atoms were calculated and the sequence of debromination was obtained.There is a close relationship between molecular structure,charge,molecular orbital and C-Br bond.All PBDEs exhibited similar debromination pathways with preferential removal of meta and para bromines.

Ultra-deep Removal of Metal Ions from Coal Tar by Complexation:Experimental Studies and Density Functional Theory Simulations

As one of the important aspects of upgrading coal tar,the ultra-deep removal of metal ions via the complexation method was investigated by screening four complexing agents and performing density functional theory(DFT)simulations.Analysis of the compositions and contents of the metallic compounds in the coal tar revealed that the main components were iron and calcium naphthenates.Direct filtration reduced the mechanical impurity content from 0.24%to 0.0752%,indicating that most of the large particles could be easily removed.Among the four complexing agents,namely,acetic acid,oxalic acid,citric acid,and ethylenediaminetetraacetic acid,oxalic acid exhibited the best demetallization performance.The DFT simulations suggested that the high performance of oxalic acid originated from its 1:1 coordination mode,rigid dicarboxyl structure,and greater binding energy.

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.

Hydrogen evolution reaction between small-sized Zr_(n)(n=2–5)clusters and water based on density functional theory

Density functional theory(DFT)is used to calculate the most stable structures of Zr_(n)(n=2-5)clusters as well as the adsorption energy values of Zr_(n)(n=2-5)clusters after adsorbing single water molecule.The results reveal that there is a significant linear relationship between the adsorption energy values and the energy gaps of the Zr_(n)(n=2-5)clusters.Furthermore,the calculations of the reaction paths between Zr_(n)(n=2-5)and single water molecule show that water molecule can react with Zr_(n)(n=2-5)clusters to dissociate,producing hydrogen,and O atoms mix with the clusters to generate Zr_(n)O(n=2-5),all of which are exothermic reactions.According to the released energy,the Zr4 cluster is the most efficient in Zr_(n)(n=2-5)clusters reacting with single water molecule.The natural population analysis(NPA)and density of states(DOS)demonstrate the production of hydrogen and orbital properties in different energy ranges,respectively,jointly forecasting that Zr_(n)O(n=2-5)will probably continue to react with more water molecules.Our findings contribute to better understanding of Zr's chemical reactivity,which can conduce to the development of effective Zr-based catalysts and hydrogen-production methods.

Level density of odd-A nuclei at saddle point

Based on the covariant density functional theory,by employing the core–quasiparticle coupling(CQC)model,the nuclear level density of odd-A nuclei at the saddle point is achieved.The total level density is calculated via the convolution of the intrinsic level density and the collective level density.The intrinsic level densities are obtained in the finite-temperature covariant density functional theory,which takes into account the nuclear deformation and pairing self-consistently.For saddle points on the free energy surface in the(β_(2),γ)plane,the entropy and the associated intrinsic level density are compared with those of the global minima.By introducing a quasiparticle to the two neighboring even–even core nuclei,whose properties are determined by the five-dimensional collective Hamiltonian model,the collective levels of the odd-A nuclei are obtained via the CQC model.The total level densities of the^(234-240)U agree well with the available experimental data and Hilaire’s result.Furthermore,the ratio of the total level densities at the saddle points to those at the global minima and the ratio of the total level densities to the intrinsic level densities are discussed separately.

Calculation of microscopic nuclear level densities based on covariant density functional theory

In this study,a microscopic method for calculating the nuclear level density(NLD)based on the covariant density functional theory(CDFT)is developed.The particle-hole state density is calculated by a combinatorial method using single-particle level schemes obtained from the CDFT,and the level densities are then obtained by considering collective effects such as vibration and rotation.Our results are compared with those of other NLD models,including phenomenological,microstatisti-cal and nonrelativistic Hartree–Fock–Bogoliubov combinatorial models.This comparison suggests that the general trends among these models are essentially the same,except for some deviations among the different NLD models.In addition,the NLDs obtained using the CDFT combinatorial method with normalization are compared with experimental data,including the observed cumulative number of levels at low excitation energies and the measured NLDs.The CDFT combinatorial method yields results that are in reasonable agreement with the existing experimental data.

Spatial Agglomeration and Diffusion of Population Based on a Regional Density Function Approach:A Case Study of Shandong Province in China

Population density functions have long been used to describe the spatial structure of regional population distributions.Several studies have been conducted to examine the population distribution in Shandong Province,China,but few have applied regional density functions to the analysis.Therefore,based on the 2000,2010,and 2020 population censuses,this study used monocentric and polycentric regional density functions to study the characteristics of population agglomeration and diffusion in Shandong.This is followed by an in-depth discussion based on population growth rate data and hot-and cold-spot analyses.The results showed that the Shandong Province population was spatially unevenly distributed.Population growth rates were higher in urban centers and counties,with more significant changes in population size in the eastern coastal areas than in the inland areas.As verified in this study,the logarithmic form of the single-center regional density function R2 was greater than 0.8,which was in line with the population spatial structure of Shandong Province.During the study period,the estimated population density of the regional center and the absolute value of the regional population density gradient both increased,indicating a clear and increasing trend of centripetal agglomeration of regional centers over the study period.Overall,the R2 value of the multicenter region density function was higher than that of the single-center region density function.The polycentric regional density function showed that the population density gradient of some centers had a downward trend,which reflected the spatial development trend of outward diffusion in these centers.Meanwhile,the variation in the estimated population density and the population density gradient exhibited differences in the central population distribution patterns at different levels.

Active MoS_(2)-based electrode for green ammonia synthesis

Nitrogen electro-reduction under mild conditions is one promising alternative approach of the energyconsuming Haber-Bosch process for the artificial ammonia synthesis.One critical aspect to unlocking this technology is to discover the catalysts with high selectivity and efficiency.In this work,the N_(2)-to-NH_(3)conversion on the functional MoS_(2)is fully investigated by density functional theory calculations since the layered MoS_(2)provides the ideal platform for the elaborating copies of the nitrogenase found in nature,wherein the functionalization is achieved via basal-adsorption,basal-substitution or edge-substitution of transition metal elements.Our results reveal that the edge-functionalization is a feasible strategy for the activity promotion;however,the basal-adsorption and basal-substitution separately suffer from the electrochemical instability and the NRR inefficiency.Specifically,MoS_(2)functionalized via edge W-substitution exhibits an exceptional activity.The energetically favored reaction pathway is through the distal pathway and a limiting potential is less than 0.20 V.Overall,this work escalates the rational design of the high-effective catalysts for nitrogen fixation and provides the explanation why the predicated catalyst have a good performance,paving the guidance for the experiments.

Mg/Fe site-specific dual-doping to boost the performance of cobalt-free nickle-rich layered oxide cathode for high-energy lithium-ion batteries

Layer-type LiNi0.9Mn0.1O2is promising to be the primary cathode material for lithium-ion batteries(LIBs)due to its excellent electrochemical performance.Unfortunately,the cathode with high nickel content suffers from severely detrimental structural transformation that causes rapid capacity attenuation.Herein,site-specific dual-doping with Fe and Mg ions is proposed to enhance the structural stability of LiNi0.9Mn0.1O2.The Fe3+dopants are inserted into transition metal sites(3b)and can favorably provide additional redox potential to compensate for charge and enhance the reversibility of anionic redox.The Mg ions are doped into the Li sites(3a)and serve as O_(2)^(-)-Mg^(2+)-O_(2)^(-)pillar to reinforce the electrostatic cohesion between the two adjacent transition-metal layers,which further suppress the cracking and the generation of harmful phase transitions,ultimately improving the cyclability.The theoretical calculations,including Bader charge and crystal orbital Hamilton populations(COHP)analyses,confirm that the doped Fe and Mg can form stable bonds with oxygen and the electrostatic repulsion of O_(2)^(-)-O_(2)^(-)can be effectively suppressed,which effectively mitigates oxygen anion loss at the high delithiation state.This dual-site doping strategy offers new avenues for understanding and regulating the crystalline oxygen redox and demonstrates significant potential for designing high-performance cobalt-free nickel-rich cathodes.

Edge and lithium concentration effects on intercalation kinetics for graphite anodes

Graphite interfaces are an important part of the anode in lithium-ion batteries(LIBs),significantly influencing Li intercalation kinetics.Graphite anodes adopt different stacking sequences depending on the concentration of the intercalated Li ions.In this work,we performed first-principles calculations to comprehensively address the energetics and dynamics of Li intercalation and Li vacancy diffusion near the no n-basal edges of graphite,namely the armchair and zigzag-edges,at high Li concentration.We find that surface effects persist in stage-Ⅱ that bind Li strongly at the edge sites.However,the pronounced effect previously identified at the zigzag edge of pristine graphite is reduced in LiC_(12),penetrating only to the subsurface site,and eventually disappearing in LiC_(6).Consequently,the distinctive surface state at the zigzag edge significantly impacts and restrains the charging rate at the initial lithiation of graphite anodes,whilst diminishes with an increasing degree of lithiation.Longer diffusion time for Li hopping to the bulk site from either the zigzag edge or the armchair edge in LiC_(6) was observed during high state of charge due to charge repulsion.Effectively controlling Li occupation and diffusion kinetics at this stage is also crucial for enhancing the charge rate.

Graphene effectively activating "dead" water molecules between manganese dioxide layers in potassium-ion battery

Aqueous potassium-ion batteries(APIBs),recognized as safe and reliable new energy devices,are considered as one of the alternatives to traditional batteries.Layered MnO_(2),serving as the main cathode,exhibits a lower specific capacity in aqueous electrolytes compared to organic systems and operates through a different reaction mechanism.The application of highly conductive graphene may effectively enhance the capacity of APIBs but could complicate the potassium storage environment.In this study,a MnO_(2) cathode pre-intercalated with K~+ions and grown on graphene(KMO@rGO) was developed using the microwave hydrothermal method for APIBs.KMO@rGO achieved a specific capacity of 90 mA h g^(-1) at a current density of 0.1 A g^(-1),maintaining a capacity retention rate of>90% after 5000 cycles at 5 A g^(-1).In-situ and exsitu characterization techniques revealed the energy-storage mechanism of KMO@rGO:layered MnO_(2)traps a large amount of "dead" water molecules during K~+ions removal.However,the introduction of graphene enables these water molecules to escape during K~+ ions insertion at the cathode.The galvanostatic intermittent titration technique and density functional theory confirmed that KMO@rGO has a higher K~+ions migration rate than MnO_(2).Therefore,the capacity of this cathode depends on the interaction between dead water and K~+ions during the energy-storage reaction.The optimal structural alignment between layered MnO_(2) and graphene allows electrons to easily flow into the external circuit.Rapid charge compensation forces numerous low-solvent K~+ions to displace interlayer dead water,enhancing the capacity.This unique reaction mechanism is unprecedented in other aqueous battery studies.

Machine learning with active pharmaceutical ingredient/polymer interaction mechanism:Prediction for complex phase behaviors of pharmaceuticals and formulations

The high throughput prediction of the thermodynamic phase behavior of active pharmaceutical ingredients(APIs)with pharmaceutically relevant excipients remains a major scientific challenge in the screening of pharmaceutical formulations.In this work,a developed machine-learning model efficiently predicts the solubility of APIs in polymers by learning the phase equilibrium principle and using a few molecular descriptors.Under the few-shot learning framework,thermodynamic theory(perturbed-chain statistical associating fluid theory)was used for data augmentation,and computational chemistry was applied for molecular descriptors'screening.The results showed that the developed machine-learning model can predict the API-polymer phase diagram accurately,broaden the solubility data of APIs in polymers,and reproduce the relationship between API solubility and the interaction mechanisms between API and polymer successfully,which provided efficient guidance for the development of pharmaceutical formulations.