The ground and the lowest-lying triplet excited state geometries, electronic structures, and spectroscopic properties of three mixed-ligand Ru(II) complexes [Ru(terpy)(phen)X]+ (terpy=2,2',6',2″-terpyridine...The ground and the lowest-lying triplet excited state geometries, electronic structures, and spectroscopic properties of three mixed-ligand Ru(II) complexes [Ru(terpy)(phen)X]+ (terpy=2,2',6',2″-terpyridine, phen=l,10-phenanthroline, and X=-C-=CH (1), X=Cl (2), X-CN (3)) were investigated theoretically using the density functional theory method. The ground and excited state geometries have been fully optimized at the B3LYP/LanL2DZ and UB3LYP/LanL2DZ levels, respectively. The absorption and emission spectra of the com- plexes in CHaCN solutions were calculated by time-dependent density functional theory with the PCM solvent model. The calculated bond lengths of Ru-C, Ru-N, and Ru-Cl in the ground state agree well with the corresponding experimental results. The highest occupied molecular orbital were dominantly localized on the Ru atom and monodentate X ligand for 1 and 2, Ru atom and terpy ligand for a, while the lowest unoccupied molecular orbital were π*(terpy) type orbital. Therefore, the lowest-energy absorptions of 1 and 2 at 688 and 631 nln are attributed to a dyz (Ru)+Tr/p(X)--π* (terpy) transition with MLCT/XLCT (metal-to-ligand charge transfer/X ligand to terpy ligand charge transfer) character, whereas that of 3 at 529 nm is related to a dyz (Ru)+π(terpy)-π* (terpy) transition with MLCT and ILCT transition character. The calculated phosphorescence of three complexes at 1011 nm (1), 913 nm (2), and 838 nm (3) have similar transition properties to that of the lowest-lying absorption. It is shown that the lowest lying absorptions and emissions transition character of these Ru(II) complexes can be tuned by changing the electron-withdrawing ability of the monodentate ligand.展开更多
To investigate the non-covalent interaction between cyclodextrins (CD) and lithium ion, a stoichiometry of α-CD, β-CD, heptakis(2,6-di-O-methyl)-β-CD (DM-β-CD), or heptakis(2,3,6-tri-O-methyl)-β-CD (TM-...To investigate the non-covalent interaction between cyclodextrins (CD) and lithium ion, a stoichiometry of α-CD, β-CD, heptakis(2,6-di-O-methyl)-β-CD (DM-β-CD), or heptakis(2,3,6-tri-O-methyl)-β-CD (TM-β-CD) was mixed with lithium salt, respectively, and then incubated at room temperature for 10 min to reach the equilibrium. In posi- tive mode, the electrospray ionization mass spectrometry (ESI-MS) results demonstrated that lithium ion can conjugate to α-, β-, DM-β- or TM-β-CD and form 1:1 stoichiometric non-covalent complexes. The binding of the complexes was further confirmed by collision- induced dissociation. The dissociation constants Kdl of four complexes (Li+α-CD, Li+β- CD, Li+DM-β-CD, and Li+TM-β-CD) were determined by mass spectrometric titration. The results showed Kdl were 18.7, 26.7, 33.6, 30.5 μmol/L for the complexes of Li+ with α-CD, β-CD, DM-β-CD, and TM-β-CD, respectively. Kdl for the Li+ complexes of/3-CD is smaller than that of DM-β-CD due to its steric effect of the partial substituted -CH3. The Kdl for the Li+ complexes of DM-β-CD is nearly in agreement with that of TM-β-CD, indicating Li+ is more likely to locate in the small rim of DM-β-CD's hydrophobic cavity. The DFT results showed through electrostatic interaction, one Li+ can strongly conjugate to four neighboring oxygen atoms. For the (α-CD+Li)+ complex, one Li+ may also situate the small rim of α-CD's hydrophobic cavity to form a non-specific host-guest complex.展开更多
The nitrogen reduction reaction(NRR)using new and efficient electrocatalysts is a promising al‐ternative to the traditional Haber‐Bosch process.Nevertheless,it remains a challenge to design efficient catalysts with ...The nitrogen reduction reaction(NRR)using new and efficient electrocatalysts is a promising al‐ternative to the traditional Haber‐Bosch process.Nevertheless,it remains a challenge to design efficient catalysts with improved catalytic performance.Herein,various O‐functional MXenes were investigated as NRR catalysts by a combination of density functional theory calculations and least absolute shrinkage and selection operator(LASSO)regression.Nb_(3)C_(2)O_(X) has been regarded as a promising catalyst for the NRR because of its stability,activity,and selectivity.The poten‐tial‐determining step is*NH_(2) hydrogenation to*NH3 with a limiting potential of-0.45 V.Further‐more,via LASSO regression,the descriptors and equations fitting the relationship between the properties of O‐functional MXenes and NRR activity have been proposed.This work not only pro‐vides a rational design strategy for catalysts but also provides machine learning data for further investigation.展开更多
KSSOLV(Kohn-Sham Solver)is a MATLAB(Matrix Laboratory)toolbox for solving the Kohn-Sham density functional theory(KS-DFT)with the plane-wave basis set.In the KS-DFT calculations,the most expensive part is commonly the...KSSOLV(Kohn-Sham Solver)is a MATLAB(Matrix Laboratory)toolbox for solving the Kohn-Sham density functional theory(KS-DFT)with the plane-wave basis set.In the KS-DFT calculations,the most expensive part is commonly the diagonalization of Kohn-Sham Hamiltonian in the self-consistent field(SCF)scheme.To enable a personal computer to perform medium-sized KS-DFT calculations that contain hundreds of atoms,we present a hybrid CPU-GPU implementation to accelerate the iterative diagonalization algorithms implemented in KSSOLV by using the MATLAB built-in Parallel Computing Toolbox.We compare the performance of KSSOLV-GPU on three types of GPU,including RTX3090,V100,and A100,with conventional CPU implementation of KSSOLV respectively and numerical results demonstrate that hybrid CPU-GPU implementation can achieve a speedup of about 10 times compared with sequential CPU calculations for bulk silicon systems containing up to 128 atoms.展开更多
A new potential energy surface (PES) for the atmospheric formation of sulfuric acid from OH+SO2 is investigated using density functional theory and high-level ab initio molecular orbital theory. A pathway focused o...A new potential energy surface (PES) for the atmospheric formation of sulfuric acid from OH+SO2 is investigated using density functional theory and high-level ab initio molecular orbital theory. A pathway focused on the new PES assumes the reaction to take place between the radical complex SO3.HO2 and H2O. The unusual stability of SO3.HO2 is the principal basis of the new pathway, which has the same final outcome as the current reaction mechanism in the literature but it avoids the production and complete release of SO3. The entire reaction pathway is composed of three consecutive elementary steps: (1) HOSO2+O2-+SO3.HO2, (2) SO3.HO2+H20-+SO3·H2O·HO2, (3) SO3.H20.HO2-+H2SO4+HO2. All three steps have small energy barriers, under 10 kcal/rnol, and are exotherrnic, and the new pathway is there- fore favorable both kinetically and therrnodynarnically. As a key step of the reactions, step (3), HO2 serves as a bridge molecule for low-barrier hydrogen transfer in the hydrolysis of SO3. Two significant atmospheric implications are expected frorn the present study. First, SO3 is not released from the oxidation of SO2 by OH radical in the atmosphere. Second, the conversion of SO2 into sulfuric acid is weakly dependent on the humidity of air.展开更多
Nanosystems play an important role in many applications.Due to their complexity,it is challenging to accurately characterize their structure and properties.An important means to reach such a goal is computational simu...Nanosystems play an important role in many applications.Due to their complexity,it is challenging to accurately characterize their structure and properties.An important means to reach such a goal is computational simulation,which is grounded on ab initio electronic structure calculations.Low scaling and accurate electronic-structure algorithms have been developed in recent years.Especially,the efficiency of hybrid density functional calculations for periodic systems has been significantly improved.With electronic structure information,simulation methods can be developed to directly obtain experimentally comparable data.For example,scanning tunneling microscopy images can be effectively simulated with advanced algorithms.When the system we are interested in is strongly coupled to environment,such as the Kondo effect,solving the hierarchical equations of motion turns out to be an effective way of computational characterization.Furthermore,the first principles simulation on the excited state dynamics rapidly emerges in recent years,and nonadiabatic molecular dynamics method plays an important role.For nanosystem involved chemical processes,such as graphene growth,multiscale simulation methods should be developed to characterize their atomic details.In this review,we review some recent progresses in methodology development for computational characterization of nanosystems.Advanced algorithms and software are essential for us to better understand of the nanoworld.展开更多
The structural and thermodynamic properties of TiAI intermetallics under high pressure have been investigated by ab initio plane-wave pseudopotential density functional theory method. It is found that the ratio of lat...The structural and thermodynamic properties of TiAI intermetallics under high pressure have been investigated by ab initio plane-wave pseudopotential density functional theory method. It is found that the ratio of lattice parameter c to a keeps almost constant with a value of 1.02 under the pressure from 0 to 20 GPa, which agrees well with the experimental results. With the pressure increasing from 20 to 45 GPa the values of c/a decrease almost linearly from 1.02 to 0.99. These calculated results indicate under low pressure the variation rate for a-axis is almost the same to that for c-axis, but under higher pressure the variation for a-axis is smaller than along e-axis. Through the quasi-harmonic Debye model, the equation of state (EOS) of TiAI intermetallies, as well as the thermal expansion and heat capacity at various pressures and temperatures are also studied.展开更多
With the application of the special properties of strongly stationary m-dependent series,this paper is concerned with the empirical likelihood confidence intervals of density func- tion under m-dependent series.The li...With the application of the special properties of strongly stationary m-dependent series,this paper is concerned with the empirical likelihood confidence intervals of density func- tion under m-dependent series.The limit distribution of empirical likelihood ratio statistics is given out,and the empirical likelihood confidence intervals of parameters can be constructed.A simulation study is conducted to show the finite sample performance of the empirical likelihood based method.展开更多
文摘The ground and the lowest-lying triplet excited state geometries, electronic structures, and spectroscopic properties of three mixed-ligand Ru(II) complexes [Ru(terpy)(phen)X]+ (terpy=2,2',6',2″-terpyridine, phen=l,10-phenanthroline, and X=-C-=CH (1), X=Cl (2), X-CN (3)) were investigated theoretically using the density functional theory method. The ground and excited state geometries have been fully optimized at the B3LYP/LanL2DZ and UB3LYP/LanL2DZ levels, respectively. The absorption and emission spectra of the com- plexes in CHaCN solutions were calculated by time-dependent density functional theory with the PCM solvent model. The calculated bond lengths of Ru-C, Ru-N, and Ru-Cl in the ground state agree well with the corresponding experimental results. The highest occupied molecular orbital were dominantly localized on the Ru atom and monodentate X ligand for 1 and 2, Ru atom and terpy ligand for a, while the lowest unoccupied molecular orbital were π*(terpy) type orbital. Therefore, the lowest-energy absorptions of 1 and 2 at 688 and 631 nln are attributed to a dyz (Ru)+Tr/p(X)--π* (terpy) transition with MLCT/XLCT (metal-to-ligand charge transfer/X ligand to terpy ligand charge transfer) character, whereas that of 3 at 529 nm is related to a dyz (Ru)+π(terpy)-π* (terpy) transition with MLCT and ILCT transition character. The calculated phosphorescence of three complexes at 1011 nm (1), 913 nm (2), and 838 nm (3) have similar transition properties to that of the lowest-lying absorption. It is shown that the lowest lying absorptions and emissions transition character of these Ru(II) complexes can be tuned by changing the electron-withdrawing ability of the monodentate ligand.
文摘To investigate the non-covalent interaction between cyclodextrins (CD) and lithium ion, a stoichiometry of α-CD, β-CD, heptakis(2,6-di-O-methyl)-β-CD (DM-β-CD), or heptakis(2,3,6-tri-O-methyl)-β-CD (TM-β-CD) was mixed with lithium salt, respectively, and then incubated at room temperature for 10 min to reach the equilibrium. In posi- tive mode, the electrospray ionization mass spectrometry (ESI-MS) results demonstrated that lithium ion can conjugate to α-, β-, DM-β- or TM-β-CD and form 1:1 stoichiometric non-covalent complexes. The binding of the complexes was further confirmed by collision- induced dissociation. The dissociation constants Kdl of four complexes (Li+α-CD, Li+β- CD, Li+DM-β-CD, and Li+TM-β-CD) were determined by mass spectrometric titration. The results showed Kdl were 18.7, 26.7, 33.6, 30.5 μmol/L for the complexes of Li+ with α-CD, β-CD, DM-β-CD, and TM-β-CD, respectively. Kdl for the Li+ complexes of/3-CD is smaller than that of DM-β-CD due to its steric effect of the partial substituted -CH3. The Kdl for the Li+ complexes of DM-β-CD is nearly in agreement with that of TM-β-CD, indicating Li+ is more likely to locate in the small rim of DM-β-CD's hydrophobic cavity. The DFT results showed through electrostatic interaction, one Li+ can strongly conjugate to four neighboring oxygen atoms. For the (α-CD+Li)+ complex, one Li+ may also situate the small rim of α-CD's hydrophobic cavity to form a non-specific host-guest complex.
文摘The nitrogen reduction reaction(NRR)using new and efficient electrocatalysts is a promising al‐ternative to the traditional Haber‐Bosch process.Nevertheless,it remains a challenge to design efficient catalysts with improved catalytic performance.Herein,various O‐functional MXenes were investigated as NRR catalysts by a combination of density functional theory calculations and least absolute shrinkage and selection operator(LASSO)regression.Nb_(3)C_(2)O_(X) has been regarded as a promising catalyst for the NRR because of its stability,activity,and selectivity.The poten‐tial‐determining step is*NH_(2) hydrogenation to*NH3 with a limiting potential of-0.45 V.Further‐more,via LASSO regression,the descriptors and equations fitting the relationship between the properties of O‐functional MXenes and NRR activity have been proposed.This work not only pro‐vides a rational design strategy for catalysts but also provides machine learning data for further investigation.
基金supported by the National Natural Science Foundation of China (No.21688102,No.21803066,and No.22003061)the Chinese Academy of Sciences Pioneer Hundred Talents Program (KJ2340000031,KJ2340007002)+7 种基金the National Key Research and Development Program of China(2016YFA0200604)the Anhui Initiative in Quantum Information Technologies (AHY090400)the Strategic Priority Research of Chinese Academy of Sciences(XDC01040100)CAS Project for Young Scientists in Basic Research (YSBR-005)the Fundamental Research Funds for the Central Universities (WK2340000091,WK2060000018)the Hefei National Laboratory for Physical Sciences at the Microscale (SK2340002001)the Research Start-Up Grants (KY2340000094)the Academic Leading Talents Training Program(KY2340000103) from University of Science and Technology of China
文摘KSSOLV(Kohn-Sham Solver)is a MATLAB(Matrix Laboratory)toolbox for solving the Kohn-Sham density functional theory(KS-DFT)with the plane-wave basis set.In the KS-DFT calculations,the most expensive part is commonly the diagonalization of Kohn-Sham Hamiltonian in the self-consistent field(SCF)scheme.To enable a personal computer to perform medium-sized KS-DFT calculations that contain hundreds of atoms,we present a hybrid CPU-GPU implementation to accelerate the iterative diagonalization algorithms implemented in KSSOLV by using the MATLAB built-in Parallel Computing Toolbox.We compare the performance of KSSOLV-GPU on three types of GPU,including RTX3090,V100,and A100,with conventional CPU implementation of KSSOLV respectively and numerical results demonstrate that hybrid CPU-GPU implementation can achieve a speedup of about 10 times compared with sequential CPU calculations for bulk silicon systems containing up to 128 atoms.
基金partially funded by National Science Foundation of the United States(No.1012994)by California State University,Fullerton
文摘A new potential energy surface (PES) for the atmospheric formation of sulfuric acid from OH+SO2 is investigated using density functional theory and high-level ab initio molecular orbital theory. A pathway focused on the new PES assumes the reaction to take place between the radical complex SO3.HO2 and H2O. The unusual stability of SO3.HO2 is the principal basis of the new pathway, which has the same final outcome as the current reaction mechanism in the literature but it avoids the production and complete release of SO3. The entire reaction pathway is composed of three consecutive elementary steps: (1) HOSO2+O2-+SO3.HO2, (2) SO3.HO2+H20-+SO3·H2O·HO2, (3) SO3.H20.HO2-+H2SO4+HO2. All three steps have small energy barriers, under 10 kcal/rnol, and are exotherrnic, and the new pathway is there- fore favorable both kinetically and therrnodynarnically. As a key step of the reactions, step (3), HO2 serves as a bridge molecule for low-barrier hydrogen transfer in the hydrolysis of SO3. Two significant atmospheric implications are expected frorn the present study. First, SO3 is not released from the oxidation of SO2 by OH radical in the atmosphere. Second, the conversion of SO2 into sulfuric acid is weakly dependent on the humidity of air.
基金supported by the Ministry of Science and Technology(No.2016YFA0200604)。
文摘Nanosystems play an important role in many applications.Due to their complexity,it is challenging to accurately characterize their structure and properties.An important means to reach such a goal is computational simulation,which is grounded on ab initio electronic structure calculations.Low scaling and accurate electronic-structure algorithms have been developed in recent years.Especially,the efficiency of hybrid density functional calculations for periodic systems has been significantly improved.With electronic structure information,simulation methods can be developed to directly obtain experimentally comparable data.For example,scanning tunneling microscopy images can be effectively simulated with advanced algorithms.When the system we are interested in is strongly coupled to environment,such as the Kondo effect,solving the hierarchical equations of motion turns out to be an effective way of computational characterization.Furthermore,the first principles simulation on the excited state dynamics rapidly emerges in recent years,and nonadiabatic molecular dynamics method plays an important role.For nanosystem involved chemical processes,such as graphene growth,multiscale simulation methods should be developed to characterize their atomic details.In this review,we review some recent progresses in methodology development for computational characterization of nanosystems.Advanced algorithms and software are essential for us to better understand of the nanoworld.
基金Support from Ph. D. Program Foundation (B2009-59)the National Science Foundations of China under Grant No. 50802024+1 种基金Henan Educational Committee under Grant No. 2011A140007Young Key Teacher by Henan Polytechnic University
文摘The structural and thermodynamic properties of TiAI intermetallics under high pressure have been investigated by ab initio plane-wave pseudopotential density functional theory method. It is found that the ratio of lattice parameter c to a keeps almost constant with a value of 1.02 under the pressure from 0 to 20 GPa, which agrees well with the experimental results. With the pressure increasing from 20 to 45 GPa the values of c/a decrease almost linearly from 1.02 to 0.99. These calculated results indicate under low pressure the variation rate for a-axis is almost the same to that for c-axis, but under higher pressure the variation for a-axis is smaller than along e-axis. Through the quasi-harmonic Debye model, the equation of state (EOS) of TiAI intermetallies, as well as the thermal expansion and heat capacity at various pressures and temperatures are also studied.
基金11th Five-Year Plan Social Science Project of Office of Education of Jilin Province(No.2007235)
文摘With the application of the special properties of strongly stationary m-dependent series,this paper is concerned with the empirical likelihood confidence intervals of density func- tion under m-dependent series.The limit distribution of empirical likelihood ratio statistics is given out,and the empirical likelihood confidence intervals of parameters can be constructed.A simulation study is conducted to show the finite sample performance of the empirical likelihood based method.