As a novel class of purely organic fluores-cent materials,multiple resonance thermal-ly activated delayed fluorescence(MR-TADF)compounds hold significant promise for next-generation display technologies.The efficiency...As a novel class of purely organic fluores-cent materials,multiple resonance thermal-ly activated delayed fluorescence(MR-TADF)compounds hold significant promise for next-generation display technologies.The efficiency of exciton utilization and the overall performance of organic light-emit-ting devices are closely linked to the singlet-triplet energy gap(ΔE_(ST))of MR-TADF emitters.Identifying an economic and accu-rate theoretical approach to predictΔE_(ST)would be beneficial for high-throughput screening and facilitate the inverse design of MR-TADF molecules.In this study,we evaluated the S_(1)state energy(E(S_(1))),T_(1)state ener-gy(E(T_(1))),andΔE_(ST)using three different physical interpretations:adiabatic excitation ener-gy,vertical absorption energy,and vertical emission energy.We employed the time-depen-dent density functional theory(TDDFT)and delta self-consistent field(ΔSCF)methods to calculate E(S_(1)),E(T_(1)),andΔE_(ST)for 20 MR-TADF molecules reported in the literature.We compared these calculated values with experimental data obtained from fluorescence spec-troscopy at room-temperature(or 77 K)and phosphorescence spectroscopy conducted at 77 K.Our findings indicate that the vertical absorption energy at the S0 state minimum,deter-mined by theΔSCF method,accurately predicts the S_(1)state energy.Similarly,the vertical absorption energy at the S0 state minimum,calculated using the TDDFT method,effectively predicts the T_(1)state energy.TheΔE_(ST)derived from the difference between these two excita-tion energies exhibited the smallest mean absolute error of only 0.039 eV compared to the ex-perimental values.This combination represents the most accurate and cost-effective method reported to date for predicting theΔE_(ST)of MR-TADF molecules,and can be integrated into AI-driven inverse design workflows for new emitters.展开更多
Ab initio modeling of dynamic structure factors(DSF)and related density response properties in the warm dense matter(WDM)regime is a challenging computational task.The DSF,convolved with a probing X-ray beam and instr...Ab initio modeling of dynamic structure factors(DSF)and related density response properties in the warm dense matter(WDM)regime is a challenging computational task.The DSF,convolved with a probing X-ray beam and instrument function,is measured in X-ray Thom-son scattering(XRTS)experiments,which allow the study of electronic structure properties at the microscopic level.Among the various ab initio methods,linear-response time-dependent density-functional theory(LR-TDDFT)is a key framework for simulating the DSF.The standard approach in LR-TDDFT for computing the DSF relies on the orbital representation.A significant drawback of this method is the unfavorable scaling of the number of required empty bands as the wavenumber increases,making LR-TDDFT impractical for modeling XRTS measurements over large energy scales,such as in backward scattering geometry.In this work,we consider and test an alternative approach to LR-TDDFT that employs the Liouville–Lanczos(LL)method for simulating the DSF of WDM.This approach does not require empty states and allows the DSF at large momentum transfer values and over a broad frequency range to be accessed.We compare the results obtained from the LL method with those from the solution of Dyson’s equation using the standard LR-TDDFT within the projector augmented-wave formalism for isochorically heated aluminum and warm dense hydrogen.Additionally,we utilize exact path integral Monte Carlo results for the imaginary-time density-density correlation function(ITCF)of warm dense hydrogen to rigorously benchmark the LL approach.We discuss the application of the LL method for calculating DSFs and ITCFs at different wavenumbers,the effects of pseudopotentials,and the role of Lorentzian smearing.The successful validation of the LL method under WDM conditions makes it a valuable addition to the ab initio simulation landscape,supporting experimental efforts and advancing WDM theory.展开更多
We report theoretical studies on the plasmon resonances in linear Au atomic chains by using ab initio time- dependent density functional theory. The dipole responses are investigated each as a function of chain length...We report theoretical studies on the plasmon resonances in linear Au atomic chains by using ab initio time- dependent density functional theory. The dipole responses are investigated each as a function of chain length. They converge into a single resonance in the longitudinal mode but split into two transverse modes. As the chain length increases, the longitudinal plasmon mode is redshifted in energy while the transverse modes shift in the opposite direction （blueshifts）. In addition, the energy gap between the two transverse modes reduces with chain length increasing. We find that there are unique characteristics, different from those of other metallic chains. These characteristics are crucial to atomic-scale engineering of single-molecule sensing, optical spectroscopy, and so on.展开更多
We propose a method for calculating the nonradiative decay rates for polyatomic molecules including anharmonic effects of the potential energy surface(PES)in the Franck-Condon region.The method combines the n-mode rep...We propose a method for calculating the nonradiative decay rates for polyatomic molecules including anharmonic effects of the potential energy surface(PES)in the Franck-Condon region.The method combines the n-mode repre-sentation method to construct the ab initio PES and the nearly exact time-dependent density matrix renormalization group method(TD-DMRG)to simulate quantum dynamics.In addition,in the framework of TD-DMRG,we further develop an algorithm to calculate the final-state-resolved rate coefficient which is very useful to analyze the contribution from each vibrational mode to the transition process.We use this method to study the internal conversion(IC)process of azulene after taking into account the anharmonicity of the ground state PES.The results show that even for this semi-rigid molecule,the intramode anharmonicity enhances the IC rate significantly,and after considering the two-mode coupling effect,the rate increases even further.The reason is that the anharmonicity enables the C-H vibrations to receive electronic energy while C-H vibrations do not contribute on the harmonic PES as the Huang-Rhys factor is close to 0.展开更多
Chiroptical properties including electronic circular dichroism（ECD） and optical rotatory dispersion（ORD） of artemisinin and artemether have been fully studied using quantum-chemical calculation based on time-depen...Chiroptical properties including electronic circular dichroism（ECD） and optical rotatory dispersion（ORD） of artemisinin and artemether have been fully studied using quantum-chemical calculation based on time-dependent density functional theory.Both theoretical ECD and ORD of these two compounds were in good match with the experimental data.ECD spectrum of artemether could be totally attributed to the peroxide group,and that of artemisinin was an overlay of contribution from δ-lactone and peroxide moieties,which leading to a positive maximum at 260 nm.Our results showed that peroxide group could produce a broad ECD band in the far-UV region originated from electron transitions of HOMO →LUMO,HOMO-1 →LUMO and HOMO-2 →LUMO in the case of artemether.This work provided a theoretical interpretation of the ECD behavior of peroxide bond.展开更多
Time-dependent density functional theory(TDDFT)method is used to investigate the details of the excited state intramolecular proton transfer(ESIPT)process and the mechanism for temperature effect on the Enol^(*)/Keto^...Time-dependent density functional theory(TDDFT)method is used to investigate the details of the excited state intramolecular proton transfer(ESIPT)process and the mechanism for temperature effect on the Enol^(*)/Keto^(*)emission ratio for the Me_(2)N-substited flavonoid(MNF)compound.The geometric structures of the S_(0) and S_(1) states are denoted as the Enol,Enol^(*),and Keto*.In addition,the absorption and fluorescence peaks are also calculated.It is noted that the calculated large Stokes shift is in good agreement with the experimental result.Furthermore,our results confirm that the ESIPT process happens upon photoexcitation,which is distinctly monitored by the formation and disappearance of the characteristic peaks of infrared(IR)spectra involved in the proton transfer and in the potential energy curves.Besides,the calculations of highest occupied molecular orbital(HOMO)and lowest unoccupied molecular orbital(LUMO)reveal that the electronegativity change of proton acceptor due to the intramolecular charge redistribution in the S_(1) state induces the ESIPT.Moreover,the thermodynamic calculation for the MNF shows that the Enol^(*)/Keto^(*)emission ratio decreasing with temperature increasing arises from the barrier lowering of ESIPT.展开更多
By using first-principles simulations based on time-dependent density functional theory, the chemical reaction of an HCl molecule encapsulated in C60 induced by femtosecond laser pulses is observed. The H atom starts ...By using first-principles simulations based on time-dependent density functional theory, the chemical reaction of an HCl molecule encapsulated in C60 induced by femtosecond laser pulses is observed. The H atom starts to leave the Cl atom and is reflected by the C60 wall. The coherent nuclear dynamic behaviors of bond breakage and recombination of the HCl molecule occurring in both polarized parallel and perpendicular to the H-Cl bond axis are investigated. The radial oscillation is also found in the two polarization directions of the laser pulse. The relaxation time of the H-Cl bond lengths in transverse polarization is slow in comparison with that in longitudinal polarization. Those results are important for studying the dynamics of the chemical bond at an atomic level.展开更多
Electronic properties, such as HOMO and LUMO energies, band gaps, ionization potential （IP） and electron affinity （EA） of 2,7- and 3,6-1inked carbazole trimers, two conjugated oligomcrs with different linkages of ...Electronic properties, such as HOMO and LUMO energies, band gaps, ionization potential （IP） and electron affinity （EA） of 2,7- and 3,6-1inked carbazole trimers, two conjugated oligomcrs with different linkages of carbazole, were studicd by the density functional theory with Becke-Lee-Young-Parr composite exchange correlation functional （B3LYP）. The absorption spectra of these compounds were also investigated by time-dependent density functional theory （TD-DFT） with 6-3 IG＊ basis set. The calculated results indicated that the HOMO and LUMO of the 2,7- and 3,6-1inked carbazole trimers are both slightly destabilized on going from methyl substitution to sec-butyl substitution. Both IP and EA exhibit their good hole-transporting but poor electronaccepting ability. The presence of alkyl groups on the nitrogen atoms does not affect the intra-chain electronic delocalization along the molecular frame. Thus no significant effect on the band gap and absorption spectra of compounds has been found.展开更多
The time-dependent density functional-based tight-bind （TD-DFTB） method is implemented on the multi-core and the graphical processing unit （GPU） system for excited state calcu-lations of large system with hundreds...The time-dependent density functional-based tight-bind （TD-DFTB） method is implemented on the multi-core and the graphical processing unit （GPU） system for excited state calcu-lations of large system with hundreds or thousands of atoms. Sparse matrix and OpenMP multithreaded are used for building the Hamiltonian matrix. The diagonal of the eigenvalue problem in the ground state is implemented on the GPUs with double precision. The GPU- based acceleration fully preserves all the properties, and a considerable total speedup of 8.73 can be achieved. A Krylov-space-based algorithm with the OpenMP parallel and CPU acceleration is used for finding the lowest eigenvalue and eigenvector of the large TDDFT matrix, which greatly reduces the iterations taken and the time spent on the excited states eigenvalue problem. The Krylov solver with the GPU acceleration of matrix-vector product can converge quickly to obtain the final result and a notable speed-up of 206 times can be observed for system size of 812 atoms. The calculations on serials of small and large systems show that the fast TD-DFTB code can obtain reasonable result with a much cheaper computational requirement compared with the first-principle results of CIS and full TDDFT calculation.展开更多
With the frame of the time-dependent local density approximation, an efficient description of the optical response of clusters has been used to study the photo-absorption cross section of Na2 and Na4 clusters. It is s...With the frame of the time-dependent local density approximation, an efficient description of the optical response of clusters has been used to study the photo-absorption cross section of Na2 and Na4 clusters. It is shown that our calculated results are in good agreement with the experiment. In addition, our calculated spectrum for the Na4 cluster is in better agreement with experiment than the GW absorption spectrum.展开更多
By using first-principles simulations based on time-dependent density functional theory,the chemical reaction of an HCl molecule encapsulated in C60induced by femtosecond laser pulses is observed.The H atom starts to ...By using first-principles simulations based on time-dependent density functional theory,the chemical reaction of an HCl molecule encapsulated in C60induced by femtosecond laser pulses is observed.The H atom starts to leave the Cl atom and is reflected by the C60wall.The coherent nuclear dynamic behaviors of bond breakage and recombination of the HCl molecule occurring in both polarized parallel and perpendicular to the H–Cl bond axis are investigated.The radial oscillation is also found in the two polarization directions of the laser pulse.The relaxation time of the H–Cl bond lengths in transverse polarization is slow in comparison with that in longitudinal polarization.Those results are important for studying the dynamics of the chemical bond at an atomic level.展开更多
Backfill is routinely adopted as a ground support measure for underground mines.However,ground stability enhancement by backfill has received limited research attention.This is likely to be because of the conventional...Backfill is routinely adopted as a ground support measure for underground mines.However,ground stability enhancement by backfill has received limited research attention.This is likely to be because of the conventional assumption that the fill material exhibits a significantly lower stiffness than the host rocks.Significantly,a recent pioneering work revealed the time-dependent ground stability around a backfilled stope with vertical walls through numerical modeling.In practice,underground stopes typically exhibit a higher or lower degree of inclination.This alters the stress state in peripheral rocks and may induce severe instability and dilution,particularly in stope-hanging walls.Hence,it is imperative to analyze the time-dependent ground stability of inclined backfilled stopes for backfill structure design.Therefore,comprehensive numerical simulations were performed using FLAC3D to address this knowledge deficiency by incorporating a coupled analysis of the backfill consolidation behavior and long-term creep deformation in surrounding rocks.The ground stability was evaluated based on the confinement effectiveness,strength-stress ratio,stress path relative to the yield surface,and time-dependent stress redistribution in the rocks.A parametric study revealed that the inclination angle of the backfilled stope reduced the confinement effectiveness in the host rocks when the wall creep was minor.This exacerbated the rock mass sloughing potential.However,a backfilled stope with a shallower dip angle achieved superior ground stability enhancement when the creep deformation was substantial,by applying a more significant compression on the backfill and effectively mobilizing its passive support performance during consolidation.Additional simulations were conducted to analyze the effects of stope height and width,mine depth,mechanical properties of rocks,backfill compressibility,and filling gap on the time-dependent stress redistribution and stability around the inclined backfilled stope.展开更多
In deep coal mining,surrounding rock is subjected to both high in-situ stress and intense mining disturbances,leading to significant time-dependent behavior.Accurately capturing this behavior is essential for predicti...In deep coal mining,surrounding rock is subjected to both high in-situ stress and intense mining disturbances,leading to significant time-dependent behavior.Accurately capturing this behavior is essential for predicting long-term roadway stability,necessitating the development of a reliable constitutive creep model and numerical simulation approach.In this study,creep experiments were conducted on pre-damaged rock with varying initial damage levels to investigate the time-dependent mechanical properties.Based on the experimental results,an accelerated-creep criterion was proposed,and an elastic-viscoplastic creep damage model(EVPCD)was established that simultaneously considers the effects of time-dependent damage and instantaneous damage caused by stress disturbances on rock creep behavior.Subsequently,the effectiveness of the proposed creep model was verified using experimental data,and the secondary development of the EVPCD model was completed based on the FLAC3D platform.Following this,a long-term stability analysis method of deep surrounding rock that accounts for excavation-and mining-induced disturbances was proposed.Using the main roadway of Xutuan Coal Mine as a case study,numerical simulations were carried out to investigate the time-dependent deformation and failure characteristics of the surrounding rock following excavation and mining disturbance.Combined with on-site monitoring of the surrounding rock damage areas,the results indicate that the EVPCD outperforms the CVISC and Nishihara models in predicting the time-dependent behavior of deep surrounding rock.展开更多
A multi-stage stress relaxation test was performed on a granodiorite sample to understand the deformation process prior to the macroscopic failure of brittle rocks,as well as the transient response during stress relax...A multi-stage stress relaxation test was performed on a granodiorite sample to understand the deformation process prior to the macroscopic failure of brittle rocks,as well as the transient response during stress relaxation.Distributed optical fiber sensing was used to measure strains across the sample surface by helically wrapping the single-mode fiber around the cylindrical sample.Close agreement was observed between the circumferential strains obtained from the optical fibers and the extensometer.The reconstructed full-field strain contours show strain heterogeneity from the crack closure phase,and the strains in the later deformation phase are dominantly localized within the former high-strain zone.The Gini coefficient was used to quantify the degree of strain localization and shows an initial increase during the crack closure phase,a decrease during the linear elastic phase,and a subsequent increase during the post-yielding phase.This behavior corresponds to a process of initial localization from an imperfect boundary condition,homogenization,and eventual relocalization prior to the macroscopic failure of the sample.The transient strain rate decay during the stress relaxation phase was quantified using the p-value in the“Omori-like"power law function.A higher initial stress at the onset of relaxation results in a lower p-value,indicating a slower strain rate decay.As the sample approaches macroscopic failure,the lowest p-value shifts from the most damaged zone to adjacent areas,suggesting stress redistribution or crack propagation in deformed crystalline rocks under stress relaxation conditions.展开更多
In-situ stress is a key parameter for underground mine design and rock stability analysis.The borehole overcoring technique is widely used for in-situ stress measurement,but the rheological recovery deformation of roc...In-situ stress is a key parameter for underground mine design and rock stability analysis.The borehole overcoring technique is widely used for in-situ stress measurement,but the rheological recovery deformation of rocks after stress relief introduces errors.To improve accuracy,this study proposes an in-situ stress solution theory that incorporates time-dependent stress relief effects.Triaxial stepwise loadingunloading rheological tests on granite and siltstone established quantitative relationships between instantaneous elastic recovery and viscoelastic recovery under different stress levels,confirming their impact on measurement accuracy.By integrating a dual-class elastic deformation recovery model,an improved in-situ stress solution theory was derived.Additionally,accounting for the nonlinear characteristics of rock masses,a determination method for time-dependent nonlinear mechanical parameters was proposed.Based on the CSIRO hollow inclusion strain cell,time-dependent strain correction equations and long-term confining pressure calibration equations were formulated.Finally,the proposed theory was successfully applied at one iron mine(736 m depth)in Xinjiang,China,and one coal mine(510 m depth)in Ningxia,China.Compared to classical theory,the calculated mean stress values showed accuracy improvements of 6.0%and 9.4%,respectively,validating the applicability and reliability of the proposed theory.展开更多
Spin excitations play a fundamental role in understanding magnetic properties of materials,and have significant technological implications for magnonic devices.However,accurately modeling these in transition-metal and...Spin excitations play a fundamental role in understanding magnetic properties of materials,and have significant technological implications for magnonic devices.However,accurately modeling these in transition-metal and rare-earth compounds remains a formidable challenge.Here,we present a fully first-principles approach for calculating spin-wave spectra based on time-dependent(TD)density-functional perturbation theory(DFPT),using nonempirical Hubbard functionals.This approach is implemented in a general noncollinear formulation,enabling the study of magnons in both collinear and noncollinear magnetic systems.Unlike methods that rely on empirical Hubbard U parameters to describe the ground state,and Heisenberg Hamiltonians for describing magnetic excitations,the methodology developed here probes directly the dynamical spin susceptibility(efficiently evaluated with TDDFPT throught the Liouville-Lanczos approach),and treats the linear variation of the Hubbard augmentation(in itself calculated non-empirically)in full at a self-consistent level.Furthermore,the method satisfies the Goldstone condition without requiring empirical rescaling of the exchange-correlation kernel or explicit enforcement of sum rules,in contrast to existing state-of-the-art techniques.We benchmark the novel computational scheme on prototypical transition-metal monoxides NiO and MnO,showing remarkable agreement with experiments and highlighting the fundamental role of these newly implemented Hubbard corrections.The method holds great promise for describing collective spin excitations in complex materials containing localized electronic states.展开更多
First-principles electron dynamics calculations can be applied in the investigation of a wide range of ultrafast phenomena in attosecond physics.They offer unique microscopic insight into light-induced ultrafast pheno...First-principles electron dynamics calculations can be applied in the investigation of a wide range of ultrafast phenomena in attosecond physics.They offer unique microscopic insight into light-induced ultrafast phenomena in both gas and condensed phases of matter,and thus,they are apowerful tool to develop our understanding of the physics of attosecond phenomena.We specifically review techniques employing time-dependent density functional theory(TDDFT)for investigating attosecond and strong-field phenomena.First,we describe this theoretical framework that enables the modeling of perturbative and non-perturbative electron dynamics in materials,including atoms,molecules,and solids.We then discuss its application to attosecond experiments,focusing on the reconstruction of attosecond beating by interference of two-photon transitions(RABBIT)measurements.Wealso briefly review first-principles calculations of optical properties of solids with TDDFT in the linear response regime and their extension to calculations of transient optical properties of solids in non-equilibrium phases,by simulating experimental pump-probe setups.We further demonstrate the application of TDDFT simulation to high-order harmonic generation in solids.First-principles calculations have predictive power,and hence they can be utilized to design future experiments to explore nonequilibrium and nonlinear ultrafast phenomena in matter and characterize and control metastable light-induced quantum states.展开更多
The behavior of interacting electrons in a perfect crystal under macroscopic external electric and magnetic fields is studied. Effective Maxwell equations for the macroscopic electric and magnetic fields are derived s...The behavior of interacting electrons in a perfect crystal under macroscopic external electric and magnetic fields is studied. Effective Maxwell equations for the macroscopic electric and magnetic fields are derived starting from time-dependent density functional theory. Effective permittivity and permeability coefficients are obtained.展开更多
Time-dependent density-functional theory(TDDFT)has been successfully applied to predict excited-state properties of isolated and periodic systems.However,it cannot address a system coupled to an environment or whose n...Time-dependent density-functional theory(TDDFT)has been successfully applied to predict excited-state properties of isolated and periodic systems.However,it cannot address a system coupled to an environment or whose number of electrons is not conserved.To tackle these problems,TDDFT needs to be extended to accommodate open systems.This paper provides a comprehensive account of the recent developments of TDDFT for open systems(TDDFT-OS),including both theoretical and practical aspects.The practicality and accuracy of a latest TDDFT-OS method is demonstrated with two numerical examples:the time-dependent electron transport through a series of quasi-one-dimensional atomic chains,and the real-time electronic dynamics on a two-dimensional graphene surface.The advancement of TDDFT-OS may lead to promising applications in various fields of chemistry,including energy conversion and heterogeneous catalysis.展开更多
Catastrophic failure in engineering structures of island reefs would occur when the tertiary creep initiates in coral reef limestone with a transition from short-to long-term load.Due to the complexity of biological s...Catastrophic failure in engineering structures of island reefs would occur when the tertiary creep initiates in coral reef limestone with a transition from short-to long-term load.Due to the complexity of biological structures,the underlying micro-behaviors involving time-dependent deformation are poorly understood.For this,an abnormal phenomenon was observed where the axial and lateral creep deformations were mutually independent by a series of triaxial tests under constant stress and strain rate conditions.The significantly large lateral creep deformation implies that the creep process cannot be described in continuum mechanics regime.Herein,it is hypothesized that sliding mechanism of crystal cleavages dominates the lateral creep deformation in coral reef limestone.Then,approaches of polarizing microscope(PM)and scanning electronic microscope(SEM)are utilized to validate the hypothesis.It shows that the sliding behavior of crystal cleavages combats with conventional creep micro-mechanisms at certain condition.The former is sensitive to time and strain rate,and is merely activated in the creep regime.展开更多
基金support provided by the National Natural Science Foundation of China(No.22273043).
文摘As a novel class of purely organic fluores-cent materials,multiple resonance thermal-ly activated delayed fluorescence(MR-TADF)compounds hold significant promise for next-generation display technologies.The efficiency of exciton utilization and the overall performance of organic light-emit-ting devices are closely linked to the singlet-triplet energy gap(ΔE_(ST))of MR-TADF emitters.Identifying an economic and accu-rate theoretical approach to predictΔE_(ST)would be beneficial for high-throughput screening and facilitate the inverse design of MR-TADF molecules.In this study,we evaluated the S_(1)state energy(E(S_(1))),T_(1)state ener-gy(E(T_(1))),andΔE_(ST)using three different physical interpretations:adiabatic excitation ener-gy,vertical absorption energy,and vertical emission energy.We employed the time-depen-dent density functional theory(TDDFT)and delta self-consistent field(ΔSCF)methods to calculate E(S_(1)),E(T_(1)),andΔE_(ST)for 20 MR-TADF molecules reported in the literature.We compared these calculated values with experimental data obtained from fluorescence spec-troscopy at room-temperature(or 77 K)and phosphorescence spectroscopy conducted at 77 K.Our findings indicate that the vertical absorption energy at the S0 state minimum,deter-mined by theΔSCF method,accurately predicts the S_(1)state energy.Similarly,the vertical absorption energy at the S0 state minimum,calculated using the TDDFT method,effectively predicts the T_(1)state energy.TheΔE_(ST)derived from the difference between these two excita-tion energies exhibited the smallest mean absolute error of only 0.039 eV compared to the ex-perimental values.This combination represents the most accurate and cost-effective method reported to date for predicting theΔE_(ST)of MR-TADF molecules,and can be integrated into AI-driven inverse design workflows for new emitters.
基金supported by the Center for Advanced Systems Understanding(CASUS),financed by Germany’s Federal Ministry of Education and Research(BMBF)and the Saxon State Government out of the State Budget approved by the Saxon State Parliamentfunding from the European Research Council(ERC)under the European Union’s Horizon 2022 research and innovation programme(Grant Agreement No.101076233,“PREXTREME”)funding from the European Union’s Just Transition Fund(JTF)within the project Röntgenlaser-Optimierung der Laserfusion(ROLF),Contract No.5086999001,co-financed by the Saxon State Government out of the State Budget approved by the Saxon State Parliament.
文摘Ab initio modeling of dynamic structure factors(DSF)and related density response properties in the warm dense matter(WDM)regime is a challenging computational task.The DSF,convolved with a probing X-ray beam and instrument function,is measured in X-ray Thom-son scattering(XRTS)experiments,which allow the study of electronic structure properties at the microscopic level.Among the various ab initio methods,linear-response time-dependent density-functional theory(LR-TDDFT)is a key framework for simulating the DSF.The standard approach in LR-TDDFT for computing the DSF relies on the orbital representation.A significant drawback of this method is the unfavorable scaling of the number of required empty bands as the wavenumber increases,making LR-TDDFT impractical for modeling XRTS measurements over large energy scales,such as in backward scattering geometry.In this work,we consider and test an alternative approach to LR-TDDFT that employs the Liouville–Lanczos(LL)method for simulating the DSF of WDM.This approach does not require empty states and allows the DSF at large momentum transfer values and over a broad frequency range to be accessed.We compare the results obtained from the LL method with those from the solution of Dyson’s equation using the standard LR-TDDFT within the projector augmented-wave formalism for isochorically heated aluminum and warm dense hydrogen.Additionally,we utilize exact path integral Monte Carlo results for the imaginary-time density-density correlation function(ITCF)of warm dense hydrogen to rigorously benchmark the LL approach.We discuss the application of the LL method for calculating DSFs and ITCFs at different wavenumbers,the effects of pseudopotentials,and the role of Lorentzian smearing.The successful validation of the LL method under WDM conditions makes it a valuable addition to the ab initio simulation landscape,supporting experimental efforts and advancing WDM theory.
基金Project supported by the National Natural Science Foundation of China (Grant Nos.11074176 and 10976019)the Research Fund for the Doctoral Program of Higher Education of China (Grant No.20100181110080)
文摘We report theoretical studies on the plasmon resonances in linear Au atomic chains by using ab initio time- dependent density functional theory. The dipole responses are investigated each as a function of chain length. They converge into a single resonance in the longitudinal mode but split into two transverse modes. As the chain length increases, the longitudinal plasmon mode is redshifted in energy while the transverse modes shift in the opposite direction （blueshifts）. In addition, the energy gap between the two transverse modes reduces with chain length increasing. We find that there are unique characteristics, different from those of other metallic chains. These characteristics are crucial to atomic-scale engineering of single-molecule sensing, optical spectroscopy, and so on.
基金supported by the National Natural Science Foundation of China through the Project "Science Center for Luminescence from Molecular Aggregates(SCELMA)" (No.21788102)the Ministry of Science and Technology of China through the National Key R&D Plan (No.2017YFA0204501)supported by the National Natural Science Foundation of China (No.22003029)
文摘We propose a method for calculating the nonradiative decay rates for polyatomic molecules including anharmonic effects of the potential energy surface(PES)in the Franck-Condon region.The method combines the n-mode repre-sentation method to construct the ab initio PES and the nearly exact time-dependent density matrix renormalization group method(TD-DMRG)to simulate quantum dynamics.In addition,in the framework of TD-DMRG,we further develop an algorithm to calculate the final-state-resolved rate coefficient which is very useful to analyze the contribution from each vibrational mode to the transition process.We use this method to study the internal conversion(IC)process of azulene after taking into account the anharmonicity of the ground state PES.The results show that even for this semi-rigid molecule,the intramode anharmonicity enhances the IC rate significantly,and after considering the two-mode coupling effect,the rate increases even further.The reason is that the anharmonicity enables the C-H vibrations to receive electronic energy while C-H vibrations do not contribute on the harmonic PES as the Huang-Rhys factor is close to 0.
基金supported by the Fundamental Research Funds for the Central Institutes of China(No.2012ZD03)
文摘Chiroptical properties including electronic circular dichroism（ECD） and optical rotatory dispersion（ORD） of artemisinin and artemether have been fully studied using quantum-chemical calculation based on time-dependent density functional theory.Both theoretical ECD and ORD of these two compounds were in good match with the experimental data.ECD spectrum of artemether could be totally attributed to the peroxide group,and that of artemisinin was an overlay of contribution from δ-lactone and peroxide moieties,which leading to a positive maximum at 260 nm.Our results showed that peroxide group could produce a broad ECD band in the far-UV region originated from electron transitions of HOMO →LUMO,HOMO-1 →LUMO and HOMO-2 →LUMO in the case of artemether.This work provided a theoretical interpretation of the ECD behavior of peroxide bond.
基金Project supported by the National Basic Research Program of China(Grant No.2013CB922204)the National Natural Science Foundation of China(Grant Nos.11574115 and 11704146)the Natural Science Foundation of Jilin Province,China(Grant No.20150101063JC)
文摘Time-dependent density functional theory(TDDFT)method is used to investigate the details of the excited state intramolecular proton transfer(ESIPT)process and the mechanism for temperature effect on the Enol^(*)/Keto^(*)emission ratio for the Me_(2)N-substited flavonoid(MNF)compound.The geometric structures of the S_(0) and S_(1) states are denoted as the Enol,Enol^(*),and Keto*.In addition,the absorption and fluorescence peaks are also calculated.It is noted that the calculated large Stokes shift is in good agreement with the experimental result.Furthermore,our results confirm that the ESIPT process happens upon photoexcitation,which is distinctly monitored by the formation and disappearance of the characteristic peaks of infrared(IR)spectra involved in the proton transfer and in the potential energy curves.Besides,the calculations of highest occupied molecular orbital(HOMO)and lowest unoccupied molecular orbital(LUMO)reveal that the electronegativity change of proton acceptor due to the intramolecular charge redistribution in the S_(1) state induces the ESIPT.Moreover,the thermodynamic calculation for the MNF shows that the Enol^(*)/Keto^(*)emission ratio decreasing with temperature increasing arises from the barrier lowering of ESIPT.
基金Project supported by the National Natural Science Foundation of China （Grant Nos. 11074176 and 10976019） and the Doctoral Program of Higher Education of China （Grant No. 20100181110080）.
文摘By using first-principles simulations based on time-dependent density functional theory, the chemical reaction of an HCl molecule encapsulated in C60 induced by femtosecond laser pulses is observed. The H atom starts to leave the Cl atom and is reflected by the C60 wall. The coherent nuclear dynamic behaviors of bond breakage and recombination of the HCl molecule occurring in both polarized parallel and perpendicular to the H-Cl bond axis are investigated. The radial oscillation is also found in the two polarization directions of the laser pulse. The relaxation time of the H-Cl bond lengths in transverse polarization is slow in comparison with that in longitudinal polarization. Those results are important for studying the dynamics of the chemical bond at an atomic level.
基金The project was supported by the Scientific Research Fund of Hunan Provincial Education Department (No. 05A002)the Prominent Mid-youth Science and Technology Foundation of Hunan Province (No. 04JJ1010)
文摘Electronic properties, such as HOMO and LUMO energies, band gaps, ionization potential （IP） and electron affinity （EA） of 2,7- and 3,6-1inked carbazole trimers, two conjugated oligomcrs with different linkages of carbazole, were studicd by the density functional theory with Becke-Lee-Young-Parr composite exchange correlation functional （B3LYP）. The absorption spectra of these compounds were also investigated by time-dependent density functional theory （TD-DFT） with 6-3 IG＊ basis set. The calculated results indicated that the HOMO and LUMO of the 2,7- and 3,6-1inked carbazole trimers are both slightly destabilized on going from methyl substitution to sec-butyl substitution. Both IP and EA exhibit their good hole-transporting but poor electronaccepting ability. The presence of alkyl groups on the nitrogen atoms does not affect the intra-chain electronic delocalization along the molecular frame. Thus no significant effect on the band gap and absorption spectra of compounds has been found.
文摘The time-dependent density functional-based tight-bind （TD-DFTB） method is implemented on the multi-core and the graphical processing unit （GPU） system for excited state calcu-lations of large system with hundreds or thousands of atoms. Sparse matrix and OpenMP multithreaded are used for building the Hamiltonian matrix. The diagonal of the eigenvalue problem in the ground state is implemented on the GPUs with double precision. The GPU- based acceleration fully preserves all the properties, and a considerable total speedup of 8.73 can be achieved. A Krylov-space-based algorithm with the OpenMP parallel and CPU acceleration is used for finding the lowest eigenvalue and eigenvector of the large TDDFT matrix, which greatly reduces the iterations taken and the time spent on the excited states eigenvalue problem. The Krylov solver with the GPU acceleration of matrix-vector product can converge quickly to obtain the final result and a notable speed-up of 206 times can be observed for system size of 812 atoms. The calculations on serials of small and large systems show that the fast TD-DFTB code can obtain reasonable result with a much cheaper computational requirement compared with the first-principle results of CIS and full TDDFT calculation.
基金The project supported by National Natural Science Foundation of China under Grant Nos. 10405025, 10575012, 10435020, and 10535010
文摘With the frame of the time-dependent local density approximation, an efficient description of the optical response of clusters has been used to study the photo-absorption cross section of Na2 and Na4 clusters. It is shown that our calculated results are in good agreement with the experiment. In addition, our calculated spectrum for the Na4 cluster is in better agreement with experiment than the GW absorption spectrum.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11074176 and 10976019)the Doctoral Program of Higher Education of China(Grant No.20100181110080)
文摘By using first-principles simulations based on time-dependent density functional theory,the chemical reaction of an HCl molecule encapsulated in C60induced by femtosecond laser pulses is observed.The H atom starts to leave the Cl atom and is reflected by the C60wall.The coherent nuclear dynamic behaviors of bond breakage and recombination of the HCl molecule occurring in both polarized parallel and perpendicular to the H–Cl bond axis are investigated.The radial oscillation is also found in the two polarization directions of the laser pulse.The relaxation time of the H–Cl bond lengths in transverse polarization is slow in comparison with that in longitudinal polarization.Those results are important for studying the dynamics of the chemical bond at an atomic level.
基金funding support from the National Natural Science Foundation of China(Nos.52304101 and 52204153)the China Postdoctoral Science Foundation(No.2023MD734215)+2 种基金the Youth Talent Support Program of Xi’an Association for Science and Technology(No.959202413070)the Key Research and Development Program of Shaanxi(No.2023-LL-QY-07)the Key Research and Development Program of Zhejiang(No.2023C03182).
文摘Backfill is routinely adopted as a ground support measure for underground mines.However,ground stability enhancement by backfill has received limited research attention.This is likely to be because of the conventional assumption that the fill material exhibits a significantly lower stiffness than the host rocks.Significantly,a recent pioneering work revealed the time-dependent ground stability around a backfilled stope with vertical walls through numerical modeling.In practice,underground stopes typically exhibit a higher or lower degree of inclination.This alters the stress state in peripheral rocks and may induce severe instability and dilution,particularly in stope-hanging walls.Hence,it is imperative to analyze the time-dependent ground stability of inclined backfilled stopes for backfill structure design.Therefore,comprehensive numerical simulations were performed using FLAC3D to address this knowledge deficiency by incorporating a coupled analysis of the backfill consolidation behavior and long-term creep deformation in surrounding rocks.The ground stability was evaluated based on the confinement effectiveness,strength-stress ratio,stress path relative to the yield surface,and time-dependent stress redistribution in the rocks.A parametric study revealed that the inclination angle of the backfilled stope reduced the confinement effectiveness in the host rocks when the wall creep was minor.This exacerbated the rock mass sloughing potential.However,a backfilled stope with a shallower dip angle achieved superior ground stability enhancement when the creep deformation was substantial,by applying a more significant compression on the backfill and effectively mobilizing its passive support performance during consolidation.Additional simulations were conducted to analyze the effects of stope height and width,mine depth,mechanical properties of rocks,backfill compressibility,and filling gap on the time-dependent stress redistribution and stability around the inclined backfilled stope.
基金funded by the National Natural Science Foundation of China(Nos.52004098,U24B2041,and 52274079)the Key Research and Development Program of Henan Province(No.251111320400)+1 种基金the Key Research Project Plan for Higher Education Institutions in Henan Province(Nos.24A570006 and 25A570002)the Scientific and Technological Research Project in Henan Province(No.242102320061).
文摘In deep coal mining,surrounding rock is subjected to both high in-situ stress and intense mining disturbances,leading to significant time-dependent behavior.Accurately capturing this behavior is essential for predicting long-term roadway stability,necessitating the development of a reliable constitutive creep model and numerical simulation approach.In this study,creep experiments were conducted on pre-damaged rock with varying initial damage levels to investigate the time-dependent mechanical properties.Based on the experimental results,an accelerated-creep criterion was proposed,and an elastic-viscoplastic creep damage model(EVPCD)was established that simultaneously considers the effects of time-dependent damage and instantaneous damage caused by stress disturbances on rock creep behavior.Subsequently,the effectiveness of the proposed creep model was verified using experimental data,and the secondary development of the EVPCD model was completed based on the FLAC3D platform.Following this,a long-term stability analysis method of deep surrounding rock that accounts for excavation-and mining-induced disturbances was proposed.Using the main roadway of Xutuan Coal Mine as a case study,numerical simulations were carried out to investigate the time-dependent deformation and failure characteristics of the surrounding rock following excavation and mining disturbance.Combined with on-site monitoring of the surrounding rock damage areas,the results indicate that the EVPCD outperforms the CVISC and Nishihara models in predicting the time-dependent behavior of deep surrounding rock.
基金support of her postdoctoral research at the GFZ Helmholtz Centre for Geosciences.P.Pan acknowledges the financial support of the National Natural Science Foundation of China(Grant No.52339001)H.Hofmann and Y.Ji acknowledge the financial support of the Helmholtz Association's Initiative and Networking Fund for the Helmholtz Young Investigator Group ARES(contract number VH-NG-1516).
文摘A multi-stage stress relaxation test was performed on a granodiorite sample to understand the deformation process prior to the macroscopic failure of brittle rocks,as well as the transient response during stress relaxation.Distributed optical fiber sensing was used to measure strains across the sample surface by helically wrapping the single-mode fiber around the cylindrical sample.Close agreement was observed between the circumferential strains obtained from the optical fibers and the extensometer.The reconstructed full-field strain contours show strain heterogeneity from the crack closure phase,and the strains in the later deformation phase are dominantly localized within the former high-strain zone.The Gini coefficient was used to quantify the degree of strain localization and shows an initial increase during the crack closure phase,a decrease during the linear elastic phase,and a subsequent increase during the post-yielding phase.This behavior corresponds to a process of initial localization from an imperfect boundary condition,homogenization,and eventual relocalization prior to the macroscopic failure of the sample.The transient strain rate decay during the stress relaxation phase was quantified using the p-value in the“Omori-like"power law function.A higher initial stress at the onset of relaxation results in a lower p-value,indicating a slower strain rate decay.As the sample approaches macroscopic failure,the lowest p-value shifts from the most damaged zone to adjacent areas,suggesting stress redistribution or crack propagation in deformed crystalline rocks under stress relaxation conditions.
基金supported by the National Science and Technology Major Project of the Ministry of Science and Technology of China(No.2024ZD1700201)the National Natural Science Foundation of China(Nos.U2034206,51974014 and 51574014)+1 种基金the Guangdong Basic and Applied Basic Research Foundation(No.2024A1515011631)the National Key Research and Development Project of China(No.2022YFC3004601)。
文摘In-situ stress is a key parameter for underground mine design and rock stability analysis.The borehole overcoring technique is widely used for in-situ stress measurement,but the rheological recovery deformation of rocks after stress relief introduces errors.To improve accuracy,this study proposes an in-situ stress solution theory that incorporates time-dependent stress relief effects.Triaxial stepwise loadingunloading rheological tests on granite and siltstone established quantitative relationships between instantaneous elastic recovery and viscoelastic recovery under different stress levels,confirming their impact on measurement accuracy.By integrating a dual-class elastic deformation recovery model,an improved in-situ stress solution theory was derived.Additionally,accounting for the nonlinear characteristics of rock masses,a determination method for time-dependent nonlinear mechanical parameters was proposed.Based on the CSIRO hollow inclusion strain cell,time-dependent strain correction equations and long-term confining pressure calibration equations were formulated.Finally,the proposed theory was successfully applied at one iron mine(736 m depth)in Xinjiang,China,and one coal mine(510 m depth)in Ningxia,China.Compared to classical theory,the calculated mean stress values showed accuracy improvements of 6.0%and 9.4%,respectively,validating the applicability and reliability of the proposed theory.
基金support by the NCCR MARVEL,a National Centre of Competence in Research,funded by the Swiss National Science Foundation(Grant number 205602)the Fellowship from the EPFL QSE Center“Many-body neural simulations of quantum materials”(Grant number 10060)supported by a grant from the Swiss National Supercomputing Centre(CSCS)under project ID s1073 and mr33(Piz Daint).
文摘Spin excitations play a fundamental role in understanding magnetic properties of materials,and have significant technological implications for magnonic devices.However,accurately modeling these in transition-metal and rare-earth compounds remains a formidable challenge.Here,we present a fully first-principles approach for calculating spin-wave spectra based on time-dependent(TD)density-functional perturbation theory(DFPT),using nonempirical Hubbard functionals.This approach is implemented in a general noncollinear formulation,enabling the study of magnons in both collinear and noncollinear magnetic systems.Unlike methods that rely on empirical Hubbard U parameters to describe the ground state,and Heisenberg Hamiltonians for describing magnetic excitations,the methodology developed here probes directly the dynamical spin susceptibility(efficiently evaluated with TDDFPT throught the Liouville-Lanczos approach),and treats the linear variation of the Hubbard augmentation(in itself calculated non-empirically)in full at a self-consistent level.Furthermore,the method satisfies the Goldstone condition without requiring empirical rescaling of the exchange-correlation kernel or explicit enforcement of sum rules,in contrast to existing state-of-the-art techniques.We benchmark the novel computational scheme on prototypical transition-metal monoxides NiO and MnO,showing remarkable agreement with experiments and highlighting the fundamental role of these newly implemented Hubbard corrections.The method holds great promise for describing collective spin excitations in complex materials containing localized electronic states.
基金supported by JSPS KAKENHI Grant Numbers JP20K14382 and JP21H01842the Cluster of Excellence 'CUI: Advanced Imaging of Matter'- EXC 2056 - project ID 390715994+5 种基金SFB-925 "Light induced dynamics and control of correlated quantum systems" – project 170620586 of the Deutsche Forschungsgemeinschaft (DFG)the Max Planck-New York City Center for Non-Equilibrium Quantum PhenomenaThis work was also supported by MEXT Promotion of Development of a Joint Usage/ Research System Project: Coalition of Universities for Research Excellence Program (CURE) Grant Number JPMXP1323015474We also acknowledge support from the Marie Sk{\l}odowska- Curie Doctoral Network TIMES, grant No. 101118915, and SPARKLE grant No. 101169225the Italian Ministry of University and Research (MUR) under the PRIN 2022 Grant No 2022PX279E_003Next Generation EUPartenariato Esteso NQSTI - Spoke 2 (THENCE-PE00000023). The Flatiron Institute is a division of the Simons Foundation. This work used computational resources of the HPC systems at the Max Planck Computing and Data Facility (MPCDF), and the Fujitsu PRIMERGY CX400M1/CX2550M5 (Oakbridge-CX) at the Information Technology Center, the University of Tokyo through the HPCI System Research Project (Project ID:hp220112).
文摘First-principles electron dynamics calculations can be applied in the investigation of a wide range of ultrafast phenomena in attosecond physics.They offer unique microscopic insight into light-induced ultrafast phenomena in both gas and condensed phases of matter,and thus,they are apowerful tool to develop our understanding of the physics of attosecond phenomena.We specifically review techniques employing time-dependent density functional theory(TDDFT)for investigating attosecond and strong-field phenomena.First,we describe this theoretical framework that enables the modeling of perturbative and non-perturbative electron dynamics in materials,including atoms,molecules,and solids.We then discuss its application to attosecond experiments,focusing on the reconstruction of attosecond beating by interference of two-photon transitions(RABBIT)measurements.Wealso briefly review first-principles calculations of optical properties of solids with TDDFT in the linear response regime and their extension to calculations of transient optical properties of solids in non-equilibrium phases,by simulating experimental pump-probe setups.We further demonstrate the application of TDDFT simulation to high-order harmonic generation in solids.First-principles calculations have predictive power,and hence they can be utilized to design future experiments to explore nonequilibrium and nonlinear ultrafast phenomena in matter and characterize and control metastable light-induced quantum states.
文摘The behavior of interacting electrons in a perfect crystal under macroscopic external electric and magnetic fields is studied. Effective Maxwell equations for the macroscopic electric and magnetic fields are derived starting from time-dependent density functional theory. Effective permittivity and permeability coefficients are obtained.
基金supported by the National Natural Science Foundation of China(21103157,21233007,and 21322305)the Fundamental Research Funds for Central Universities(2340000034 and 2340000025)the Strategic Priority Research Program(B)of the CAS(XDB01020000)
文摘Time-dependent density-functional theory(TDDFT)has been successfully applied to predict excited-state properties of isolated and periodic systems.However,it cannot address a system coupled to an environment or whose number of electrons is not conserved.To tackle these problems,TDDFT needs to be extended to accommodate open systems.This paper provides a comprehensive account of the recent developments of TDDFT for open systems(TDDFT-OS),including both theoretical and practical aspects.The practicality and accuracy of a latest TDDFT-OS method is demonstrated with two numerical examples:the time-dependent electron transport through a series of quasi-one-dimensional atomic chains,and the real-time electronic dynamics on a two-dimensional graphene surface.The advancement of TDDFT-OS may lead to promising applications in various fields of chemistry,including energy conversion and heterogeneous catalysis.
基金supported by the National Natural Science Foundation of China(Grant Nos.41877267,41877260)the Priority Research Program of the Chinese Academy of Science(Grant No.XDA13010201).
文摘Catastrophic failure in engineering structures of island reefs would occur when the tertiary creep initiates in coral reef limestone with a transition from short-to long-term load.Due to the complexity of biological structures,the underlying micro-behaviors involving time-dependent deformation are poorly understood.For this,an abnormal phenomenon was observed where the axial and lateral creep deformations were mutually independent by a series of triaxial tests under constant stress and strain rate conditions.The significantly large lateral creep deformation implies that the creep process cannot be described in continuum mechanics regime.Herein,it is hypothesized that sliding mechanism of crystal cleavages dominates the lateral creep deformation in coral reef limestone.Then,approaches of polarizing microscope(PM)and scanning electronic microscope(SEM)are utilized to validate the hypothesis.It shows that the sliding behavior of crystal cleavages combats with conventional creep micro-mechanisms at certain condition.The former is sensitive to time and strain rate,and is merely activated in the creep regime.
