For non-stationary complex dynamic systems,a standardized algorithm is developed to compute time correlation functions,addressing the limitations of traditional methods reliant on the stationary assumption.The propose...For non-stationary complex dynamic systems,a standardized algorithm is developed to compute time correlation functions,addressing the limitations of traditional methods reliant on the stationary assumption.The proposed algorithm integrates two-point and multi-point time correlation functions into a unified framework.Further,it is verified by a practical application in complex financial systems,demonstrating its potential in various complex dynamic systems.展开更多
Dear Editor,This letter addresses distributed optimization for resource allocation problems with time-varying objective functions and time-varying constraints.Inspired by the distributed average tracking(DAT)approach,...Dear Editor,This letter addresses distributed optimization for resource allocation problems with time-varying objective functions and time-varying constraints.Inspired by the distributed average tracking(DAT)approach,a distributed control protocol is proposed for optimal resource allocation.The convergence to a time-varying optimal solution within a predefined time is proved.Two numerical examples are given to illustrate the effectiveness of the proposed approach.展开更多
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.展开更多
A series of blue and blue‑green Ir(Ⅲ)complexes has been investigated theoretically to explore their electronic structures,photophysical properties,efficiency roll‑off effect,and thermal activation delayed fluorescenc...A series of blue and blue‑green Ir(Ⅲ)complexes has been investigated theoretically to explore their electronic structures,photophysical properties,efficiency roll‑off effect,and thermal activation delayed fluorescence(TADF)properties.All calculations were performed using density functional theory(DFT)and time‑dependent density functional theory(TDDFT).Calculations for electronic structures,frontier molecular orbital characteristics(which determine the efficiency roll‑off effect of the complexes),and photophysical properties were conducted using the Gaussian 09 software package.The calculation of spin‑orbit coupling matrix elements<T|HSOC|S>,which determine the TADF properties of the complexes,was performed using the ORCA software package.The calculation results show that the auxiliary ligand tetraphenylimidodiphosphinate(tpip),a strong electron‑withdrawing group,can mitigate the efficiency roll‑off effect of the complex.Furthermore,TADF is observed in one of the designed complexes,(F_(3)Phppy)_(2)Ir(tpip),where F_(3)Phppy=2‑[4‑(2,4,6‑trifluorophenyl)phenyl]pyridine.展开更多
An earthquake of Ms= 6, 9 occurred at the Gonghe, Qinghai Province, China on April 26, 1990. Three larger aftershocks took place at the same region, Ms= 5. 0 on May 7, 1990, Ms= 6. 0 on Jan. 3, 1994 and Ms= 5. 7on Feb...An earthquake of Ms= 6, 9 occurred at the Gonghe, Qinghai Province, China on April 26, 1990. Three larger aftershocks took place at the same region, Ms= 5. 0 on May 7, 1990, Ms= 6. 0 on Jan. 3, 1994 and Ms= 5. 7on Feb. 16, 1994. The long-period recordings of the main shock from China Digital Seismograph Network (CDSN) are deconvolved for the source time functions by the correspondent0 recordings of the three aftershocks asempirical Green's functions (EGFs). No matter which aftershock is taken as EGF, the relative source time functions (RSTFs) Obtained are nearly identical. The RSTFs suggest the Ms= 6. 9 event consists of at least two subevents with approximately equal size whose occurrence times are about 30 s apart, the first one has a duration of 12 s and a rise time of about 5 s, and the second one has a duration of 17 s and a rise time of about & s. COmParing the RSTFs obtained from P- and SH-phases respectively, we notice that those from SH-phases are a slightly more complex than those from p-phases, implying other finer subevents exist during the process of the main shock. It is interesting that the results from the EGF deconvolution of long-Period way form data are in good agreement with the results from the moment tensor inversion and from the EGF deconvolution of broadband waveform data. Additionally, the two larger aftershocks are deconvolved for their RSTFs. The deconvolution results show that the processes of the Ms= 6. 0 event on Jan. 3, 1994 and the Ms= 5. 7 event on Feb. 16,1994 are quite simple, both RSTFs are single impulses.The RSTFs of the Ms= 6. 9 main shock obtained from different stations are noticed to be azimuthally dependent, whose shapes are a slightly different with different stations. However, the RSTFs of the two smaller aftershocks are not azimuthally dependent. The integrations of RSTFs over the processes are quite close to each other, i. e., the scalar seismic moments estimated from different stations are in good agreement. Finally the scalar seismic moments of the three aftershocks are compared. The relative scalar seismic moment Of the three aftershocks deduced from the relative scalar seismic moments of the Ms=6. 9 main shock are very close to those inverted directly from the EGF deconvolution. The relative scalar seismic moment of the Ms =6. 9 main shock calculated using the three aftershocks as EGF are 22 (the Ms= 6. 0 aftershock being EGF), 26 (the Ms= 5. 7 aftershock being EGF) and 66 (the Ms= 5. 5 aftershock being EGF), respectively. Deducingfrom those results, the relative scalar sesimic moments of the Ms= 6. 0 to the Ms= 5. 7 events, the Ms= 6. 0 tothe Ms= 5. 5 events and the Ms= 5. 7 to the Ms= 5. 5 events are 1. 18, 3. 00 and 2. 54, respectively. The correspondent relative scalar seismic moments calculated directly from the waveform recordings are 1. 15, 3. 43, and 3. 05.展开更多
The surface subsidence is a common environmental hazard in mined-out area. Based on careful analysis of the regularity of surface subsidence in mined-out area, we proposed a new time function based on Harris curve mod...The surface subsidence is a common environmental hazard in mined-out area. Based on careful analysis of the regularity of surface subsidence in mined-out area, we proposed a new time function based on Harris curve model in consideration of the shortage of current surface subsidence time functions. By analyzing the characteristics of the new time function, we found that it could meet the dynamic process, the velocity change process and the acceleration change process during surface subsidence. Then its rationality had been verified through project cases. The results show that the proposed time function model can give a good reflection of the regularity of surface subsidence in mined-out area and can accurately predict surface subsidence. And the prediction data of the model are a little greater than measured data on condition of proper measured data quantity, which is safety in the engineering. This model provides a new method for the analysis of surface subsidence in mined-out area and reference for future prediction, and it is valuable to engineering application.展开更多
Plant capacity for water storage leads to time lags between basal stem sap flow and transpiration in various woody plants. Internal water storage depends on the sizes of woody plants. However, the changes and its infl...Plant capacity for water storage leads to time lags between basal stem sap flow and transpiration in various woody plants. Internal water storage depends on the sizes of woody plants. However, the changes and its influencing factors in time lags of basal stem flow during the development of herbaceous plants including crops remain unclear. A field experiment was conducted in an arid region of Northwest China to examine the time lag characteristics of sap flow in seed-maize and to calibrate the transpiration modeling. Cross-correlation analysis was used to estimate the time lags between stem sap flow and meteorological driving factors including solar radiation（R_s） and vapor pressure deficit of the air（VPD_（air））. Results indicate that the changes in seed-maize stem sap flow consistently lagged behind the changes in R_s and preceded the changes in VPD_（air） both on hourly and daily scales, suggesting that light-mediated stomatal closures drove sap flow responses. The time lag in the maize＇s sap flow differed significantly during different growth stages and the difference was potentially due to developmental changes in capacitance tissue and/or xylem during ontogenesis. The time lags between stem sap flow and R_s in both female plants and male plants corresponded to plant use of stored water and were independent of total plant water use. Time lags of sap flow were always longer in male plants than in female plants. Theoretically, dry soil may decrease the speed by which sap flow adjusts ahead of shifts in VPD_（air） in comparison with wet soil and also increase the speed by which sap flow adjusts to R_s. However, sap flow lags that were associated with R_s before irrigation and after irrigation in female plants did not shift. Time series analysis method provided better results for simulating seed-maize sap flow with advantages of allowing for fewer variables to be included. This approach would be helpful in improving the accuracy of estimation for canopy transpiration and conductance using meteorological measurements.展开更多
A large earthquake (Mw=7.6) occurred in Jiji (Chi-Chi), Taiwan, China on September 20, 1999, and was followed by many moderate-size shocks in the following days. Two of the largest aftershocks with the magnitudes of M...A large earthquake (Mw=7.6) occurred in Jiji (Chi-Chi), Taiwan, China on September 20, 1999, and was followed by many moderate-size shocks in the following days. Two of the largest aftershocks with the magnitudes of Mw=6.1 and Mw=6.2, respectively, were used as empirical Green's functions (EGFs) to obtain the source time functions (STFs) of the main shock from long-period waveform data of the Global Digital Seismograph Network (GDSN) including IRIS, GEOSCOPE and CDSN. For the Mw=6.1 aftershock of September 22, there were 97 pairs of phases clear enough from 78 recordings of 26 stations; for the Mw=6.2 aftershock of September 25, there were 81 pairs of phases clear enough from 72 recordings of 24 stations. For each station, 2 types of STFs were retrieved, which are called P-STF and S-STF due to being from P and S phases, respectively. Totally, 178 STF individuals were obtained for source-process analysis of the main shock. It was noticed that, in general, STFs from most of the stations had similarities except that those in special azimuths looked different or odd due to the mechanism difference between the main shock and the aftershocks; and in detail, the shapes of the STFs varied with azimuth. Both of them reflected the stability and reliability of the retrieved STFs. The comprehensive analysis of those STFs suggested that this event consisted of two sub-events, the total duration time was about 26 s, and on the average, the second event was about 7 s later than the first one, and the moment-rate amplitude of the first event was about 15% larger than that of the second one.展开更多
Based on the existing continuous borehole strain observation,the multiquadric function fitting method was used to deal with time series data. The impact of difference kernel function parameters was discussed to obtain...Based on the existing continuous borehole strain observation,the multiquadric function fitting method was used to deal with time series data. The impact of difference kernel function parameters was discussed to obtain a valuable fitting result,from which the physical connotation of the original data and its possible applications were analyzed.Meanwhile,a brief comparison was made between the results of multiquadric function fitting and polynomial fitting.展开更多
For a general normed vector space,a special optimal value function called a maximal time function is considered.This covers the farthest distance function as a special case,and has a close relationship with the smalle...For a general normed vector space,a special optimal value function called a maximal time function is considered.This covers the farthest distance function as a special case,and has a close relationship with the smallest enclosing ball problem.Some properties of the maximal time function are proven,including the convexity,the lower semicontinuity,and the exact characterizations of its subdifferential formulas.展开更多
The Region-Time-Length Algorithm (RTL algorithm) is introduced and improved in the paper. Compared with the original definition, the influence of rupture length on RTL function is emphasized and the weights of epicent...The Region-Time-Length Algorithm (RTL algorithm) is introduced and improved in the paper. Compared with the original definition, the influence of rupture length on RTL function is emphasized and the weights of epicentral distance function, time function, and rupture length function are ensured to be equal. The retrospective examinations of RTL algorithm in North China have indicated that the anomalies obtained by the improved RTL algorithm show the short or intermediate-short term precursory features in most cases. There are two types of RTL anomalous patterns before the main shock. For the I-type, the variation pattern of the VRTL, numerical values of the VRTL(x, y, z, t) function, is complete and most of them have shown the changing pattern of rising from 0 turning dropping or dropping from 0 turning rising. For the II-type, the variation pattern of VRTL is not complete, which increases or decreases quickly from 0 and there is no evident turning, the main shock generally occurs in the short period around the peak VRTL. The rising of VRTL indicates an increase of seismic activity relative to the background level, which means the enhancement of seismic activity, while the dropping of VRTL indicates the decrease of seismic activity relative to the background level, which represents the seismic quiescence to a certain extent. According to statistical examination results of RTL algorithm in North China, the methods to distinguish the intermediate and short-term anomalies and to estimate the occurrence time of the coming main shock are given in the paper. For both I and II-type RTL anomalies, the R-value, i.e., the forecasting score, is about 0.6 and 0.3 for the 3 months forecasting period and about 0.7 and 0.4 for the 6 months forecasting period. The preliminary discussion is also made for the influences of characteristic time-span t0, characteristic distance r0, and threshold magnitude M0 on computation of VRTL, as well as some other significant problems in application.展开更多
The lumped time distribution functions were proposed, which can be used for describing the dynamicsystems with two or more than two states of the end of growing polymer chain during chain addition polymerization.Numer...The lumped time distribution functions were proposed, which can be used for describing the dynamicsystems with two or more than two states of the end of growing polymer chain during chain addition polymerization.Numerical analysis of the lumped time distribution functions was carried out. The method for calculating molecularweight distribution of polymer in the stable free radical polymerization and more general cases was developed basedon the lumped time distribution functions.展开更多
This paper applies a new formulation to do moment tensor inversion for earthquakes in the Kushiro area of Japan. Comparing with conventional moment tensor inversion method, the new one takes the effect of source time ...This paper applies a new formulation to do moment tensor inversion for earthquakes in the Kushiro area of Japan. Comparing with conventional moment tensor inversion method, the new one takes the effect of source time function into consideration. For the inversion, best solution is obtained by minimizing the difference between the observed seismograms and the synthetic ones. And the best-fitting focal depth is determined from the variance reduction. The results indicate that half duration of source time function is proportional to the magnitude of earthquakes. Large earthquakes have long half duration, whereas that of moderate-small earthquakes is comparatively shorter. The focal mechanisms of all three earthquakes are of thrust fault type, which is mainly ascribed to the collision of the North American plate with the Eurasia plate in the late Cretaceous or Paleogene.展开更多
Finite time stability and stabilization are studied for hy-brid dynamic systems. By combining multiple Lyapunov function and finite time Lyapunov function, a sufficient condition of finite time stability is given for ...Finite time stability and stabilization are studied for hy-brid dynamic systems. By combining multiple Lyapunov function and finite time Lyapunov function, a sufficient condition of finite time stability is given for the system. Compared with the previ-ous works, our results have less conservativeness. Furthermore, based on the state partition of continuous and resetting parts of system, a hybrid feedback controller is constructed, which stabi-lizes the closed-loop systems in finite time. Finally, a numerical example is provided to demonstrate the effectiveness of the pro-posed method.展开更多
Escape time algorithm is an effective theoretical algorithm of constructing fractal graphics. The key of this algorithm lies in the construction of escape time function. A new escape time function is presented based o...Escape time algorithm is an effective theoretical algorithm of constructing fractal graphics. The key of this algorithm lies in the construction of escape time function. A new escape time function is presented based on the research of escape time algorithm. An accelerated escape time algorithm is carried out in this paper. The experiments have demonstrated that the new algorithm is not only as precise as the old, but also faster when it is used to construct Julia set.展开更多
In this study,application of the spectral representation method for generation of endurance time excitation functions is introduced.Using this method,the intensifying acceleration time series is generated so that its ...In this study,application of the spectral representation method for generation of endurance time excitation functions is introduced.Using this method,the intensifying acceleration time series is generated so that its acceleration response spectrum in any desired time duration is compatible with a time-scaled predefined acceleration response spectrum.For this purpose,simulated stationary acceleration time series is multiplied by the time dependent linear modulation function,then using a simple iterative scheme,it is forced to match a target acceleration response spectrum.It is shown that the generated samples have excellent conformity in low frequency,which is useful for nonlinear endurance time analysis.In the second part of this study,it is shown that this procedure can be extended to generate a set of spatially correlated endurance time excitation functions.This makes it possible to assess the performance of long structures under multi-support seismic excitation using endurance time analysis.展开更多
It has been a period of time since the concept of scenario earthquake was proposed, but this concept has rarely been used in seismic safety evaluation in China since then. Meanwhile, because of the uncertainties of ma...It has been a period of time since the concept of scenario earthquake was proposed, but this concept has rarely been used in seismic safety evaluation in China since then. Meanwhile, because of the uncertainties of magnitudes-distances pairs, there is large arbitrariness while determining the envelope function of time histories in seismic hazard analysis. In this paper, we describe a method to control the envelope functions of the time histories by introducing the most-likely combinations of magnitude and distance of the scenario earthquakes based on a probabilistic method, revise the software of the ellipse model for seismic hazard analysis, and give a computation example.展开更多
Several studies showed that the breast cancer incidence rates are higher in high-income (developed) countries, due to the link of breast cancer with several risk factors and the presence of systematic screening polici...Several studies showed that the breast cancer incidence rates are higher in high-income (developed) countries, due to the link of breast cancer with several risk factors and the presence of systematic screening policies. Some of the authors suggest that lower breast cancer incidence rates in low-income (developing) countries probably reflect international variation in hormonal factors and accessibility to early detection facilities. Recent studies showed that the breast cancer increased rapidly among women in Pakistan (a developing country) and it became the first malignancy among females of Pakistan. Although, the incidence rates may contain important evidence for understanding and control of the disease;however in Pakistan, the breast cancer incidence data have never been available in the last five decades since independence;rather, only hospital-based data are available. In this study, we intend to apply Functional Time Series (FTS) models to the breast cancer incidence rates of United State (developed country), and to see the difference between various components (age and time) of Functional Time Series (FTS) models applied independently on the breast cancer incidence rates of Karachi (Pakistan) and US. Past studies have already suggested that the incidence of US breast cancer cases was expected to increase in the coming decades. A progressive increase in the number of new cases is already predetermined by the high birth rate that occurred during the middle part of the century, and it will lead to nearly a doubling in the number of cases in about 4 decades. We also obtain 15 years predictions of breast cancer incidence rates in United States and compare them with the forecasts of incidence curves for Karachi. Development of methods for cancer incidence trend forecasting can provide a sound and accurate foundation for planning a comprehensive national strategy for optimal partitioning of research resources between the need for development of new treatments and the need for new research directed toward primary preventive measures.展开更多
Static disorder plays a crucial role in the electronic dynamics and spec-troscopy of complex molecular sys-tems.Traditionally,obtaining ob-servables averaged over static disor-der requires thousands of realiza-tions v...Static disorder plays a crucial role in the electronic dynamics and spec-troscopy of complex molecular sys-tems.Traditionally,obtaining ob-servables averaged over static disor-der requires thousands of realiza-tions via direct sampling of the dis-order distribution,leading to high computational costs.In this work,we extend the auxiliary degree-of-freedom based matrix product state(MPS)method to handle system-bath correlated thermal equilibrium initial states,which can capture static disorder effects using a one-shot quantum dynamical simulation.We validate the effectiveness of the extended method by computing the dipole-dipole time correlation function of the Holstein model relevant to the emission spectrum of molecular aggregates.Our results show that the one-shot method is very accu-rate with only a moderate increase in MPS bond dimension,thereby significantly reducing computational cost.Moreover,it enables the generation of a much larger number of samples than the conventional direct sampling method at negligible additional cost,thus reducing sta-tistical errors.This method provides a broadly useful tool for calculating equilibrium time cor-relation functions in system-bath coupled models with static disorder.展开更多
Permanent dipole moments induced high-order harmonic generation(HHG)signals offer a potential approach to producing elliptically or even circularly polarized X-ray attosecond sources.Previous studies on this topic hav...Permanent dipole moments induced high-order harmonic generation(HHG)signals offer a potential approach to producing elliptically or even circularly polarized X-ray attosecond sources.Previous studies on this topic have mainly focused on diatomic molecules such as CO and HeH.Based on this scheme,significant HHG signals in the direction perpendicular to the molecular axis can be observed in both the high-energy and low-energy regions.However,we found that the high-order harmonics induced by the permanent dipole moments of polyatomic complex molecules involve more intricate physical processes.Using time-dependent density functional theory,we simulated the dynamics of HHG from NH2COOH and NH2COSH interacting with linearly polarized lasers.We found that the harmonic signals in the direction perpendicular to the N-C bond were significantly enhanced in the high-energy photon region.Our analysis indicates that this is due to the complex molecular configuration of NH_(2)COOH and NH_(2)COSH:while the NH_(2) group has C_(2v) symmetry,both COOH and COSH groups lack this symmetry.This structural characteristic results in permanent dipole moments being felt only when electrons return to either COSH or COOH groups,but not to NH_(2) group.Additionally,our results reveal a multi-plateau structure in HHG signal along laser polarization direction,a phenomenon arising from multi-electron and multiorbital effects during interaction between complex molecule and strong laser field.展开更多
基金Project supported by the Postdoctoral Fellowship Program of China Postdoctoral Science Foundation(Grant No.GZC20231050)the National Natural Science Foundation of China(Grant Nos.12175193 and 11905183)the 13th Five-year plan for Education Science Funding of Guangdong Province(Grant No.2021GXJK349)。
文摘For non-stationary complex dynamic systems,a standardized algorithm is developed to compute time correlation functions,addressing the limitations of traditional methods reliant on the stationary assumption.The proposed algorithm integrates two-point and multi-point time correlation functions into a unified framework.Further,it is verified by a practical application in complex financial systems,demonstrating its potential in various complex dynamic systems.
基金supported by National Key Research and Development Program of China(2024YFE0214000)National Natural Science Foundation of China(62173308)+3 种基金Natural Science Foundation of Zhejiang Province of China(LRG25F030002)Zhejiang Province Leading Geese Plan(2025C01056)Jinhua Science and Technology Project(2022-1-042)Natural Science Foundation of Jiangsu Province(BK20240009).
文摘Dear Editor,This letter addresses distributed optimization for resource allocation problems with time-varying objective functions and time-varying constraints.Inspired by the distributed average tracking(DAT)approach,a distributed control protocol is proposed for optimal resource allocation.The convergence to a time-varying optimal solution within a predefined time is proved.Two numerical examples are given to illustrate the effectiveness of the proposed approach.
基金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.
文摘A series of blue and blue‑green Ir(Ⅲ)complexes has been investigated theoretically to explore their electronic structures,photophysical properties,efficiency roll‑off effect,and thermal activation delayed fluorescence(TADF)properties.All calculations were performed using density functional theory(DFT)and time‑dependent density functional theory(TDDFT).Calculations for electronic structures,frontier molecular orbital characteristics(which determine the efficiency roll‑off effect of the complexes),and photophysical properties were conducted using the Gaussian 09 software package.The calculation of spin‑orbit coupling matrix elements<T|HSOC|S>,which determine the TADF properties of the complexes,was performed using the ORCA software package.The calculation results show that the auxiliary ligand tetraphenylimidodiphosphinate(tpip),a strong electron‑withdrawing group,can mitigate the efficiency roll‑off effect of the complex.Furthermore,TADF is observed in one of the designed complexes,(F_(3)Phppy)_(2)Ir(tpip),where F_(3)Phppy=2‑[4‑(2,4,6‑trifluorophenyl)phenyl]pyridine.
文摘An earthquake of Ms= 6, 9 occurred at the Gonghe, Qinghai Province, China on April 26, 1990. Three larger aftershocks took place at the same region, Ms= 5. 0 on May 7, 1990, Ms= 6. 0 on Jan. 3, 1994 and Ms= 5. 7on Feb. 16, 1994. The long-period recordings of the main shock from China Digital Seismograph Network (CDSN) are deconvolved for the source time functions by the correspondent0 recordings of the three aftershocks asempirical Green's functions (EGFs). No matter which aftershock is taken as EGF, the relative source time functions (RSTFs) Obtained are nearly identical. The RSTFs suggest the Ms= 6. 9 event consists of at least two subevents with approximately equal size whose occurrence times are about 30 s apart, the first one has a duration of 12 s and a rise time of about 5 s, and the second one has a duration of 17 s and a rise time of about & s. COmParing the RSTFs obtained from P- and SH-phases respectively, we notice that those from SH-phases are a slightly more complex than those from p-phases, implying other finer subevents exist during the process of the main shock. It is interesting that the results from the EGF deconvolution of long-Period way form data are in good agreement with the results from the moment tensor inversion and from the EGF deconvolution of broadband waveform data. Additionally, the two larger aftershocks are deconvolved for their RSTFs. The deconvolution results show that the processes of the Ms= 6. 0 event on Jan. 3, 1994 and the Ms= 5. 7 event on Feb. 16,1994 are quite simple, both RSTFs are single impulses.The RSTFs of the Ms= 6. 9 main shock obtained from different stations are noticed to be azimuthally dependent, whose shapes are a slightly different with different stations. However, the RSTFs of the two smaller aftershocks are not azimuthally dependent. The integrations of RSTFs over the processes are quite close to each other, i. e., the scalar seismic moments estimated from different stations are in good agreement. Finally the scalar seismic moments of the three aftershocks are compared. The relative scalar seismic moment Of the three aftershocks deduced from the relative scalar seismic moments of the Ms=6. 9 main shock are very close to those inverted directly from the EGF deconvolution. The relative scalar seismic moment of the Ms =6. 9 main shock calculated using the three aftershocks as EGF are 22 (the Ms= 6. 0 aftershock being EGF), 26 (the Ms= 5. 7 aftershock being EGF) and 66 (the Ms= 5. 5 aftershock being EGF), respectively. Deducingfrom those results, the relative scalar sesimic moments of the Ms= 6. 0 to the Ms= 5. 7 events, the Ms= 6. 0 tothe Ms= 5. 5 events and the Ms= 5. 7 to the Ms= 5. 5 events are 1. 18, 3. 00 and 2. 54, respectively. The correspondent relative scalar seismic moments calculated directly from the waveform recordings are 1. 15, 3. 43, and 3. 05.
基金supported by the Key Program of the National Natural Science Foundation of China (No. 50334060)
文摘The surface subsidence is a common environmental hazard in mined-out area. Based on careful analysis of the regularity of surface subsidence in mined-out area, we proposed a new time function based on Harris curve model in consideration of the shortage of current surface subsidence time functions. By analyzing the characteristics of the new time function, we found that it could meet the dynamic process, the velocity change process and the acceleration change process during surface subsidence. Then its rationality had been verified through project cases. The results show that the proposed time function model can give a good reflection of the regularity of surface subsidence in mined-out area and can accurately predict surface subsidence. And the prediction data of the model are a little greater than measured data on condition of proper measured data quantity, which is safety in the engineering. This model provides a new method for the analysis of surface subsidence in mined-out area and reference for future prediction, and it is valuable to engineering application.
基金support from the National Key Basic Research Program of China (2016YFC0400207)the National Natural Science Foundation of China (51439006, 91425302)the 111 Program of Introducing Talents of Discipline to Universities (B14002)
文摘Plant capacity for water storage leads to time lags between basal stem sap flow and transpiration in various woody plants. Internal water storage depends on the sizes of woody plants. However, the changes and its influencing factors in time lags of basal stem flow during the development of herbaceous plants including crops remain unclear. A field experiment was conducted in an arid region of Northwest China to examine the time lag characteristics of sap flow in seed-maize and to calibrate the transpiration modeling. Cross-correlation analysis was used to estimate the time lags between stem sap flow and meteorological driving factors including solar radiation（R_s） and vapor pressure deficit of the air（VPD_（air））. Results indicate that the changes in seed-maize stem sap flow consistently lagged behind the changes in R_s and preceded the changes in VPD_（air） both on hourly and daily scales, suggesting that light-mediated stomatal closures drove sap flow responses. The time lag in the maize＇s sap flow differed significantly during different growth stages and the difference was potentially due to developmental changes in capacitance tissue and/or xylem during ontogenesis. The time lags between stem sap flow and R_s in both female plants and male plants corresponded to plant use of stored water and were independent of total plant water use. Time lags of sap flow were always longer in male plants than in female plants. Theoretically, dry soil may decrease the speed by which sap flow adjusts ahead of shifts in VPD_（air） in comparison with wet soil and also increase the speed by which sap flow adjusts to R_s. However, sap flow lags that were associated with R_s before irrigation and after irrigation in female plants did not shift. Time series analysis method provided better results for simulating seed-maize sap flow with advantages of allowing for fewer variables to be included. This approach would be helpful in improving the accuracy of estimation for canopy transpiration and conductance using meteorological measurements.
基金State Natural Science Foundation of China (49904004) and IPGP of France.Contribution No. 02FE2007, Institute of Geophysics, Ch
文摘A large earthquake (Mw=7.6) occurred in Jiji (Chi-Chi), Taiwan, China on September 20, 1999, and was followed by many moderate-size shocks in the following days. Two of the largest aftershocks with the magnitudes of Mw=6.1 and Mw=6.2, respectively, were used as empirical Green's functions (EGFs) to obtain the source time functions (STFs) of the main shock from long-period waveform data of the Global Digital Seismograph Network (GDSN) including IRIS, GEOSCOPE and CDSN. For the Mw=6.1 aftershock of September 22, there were 97 pairs of phases clear enough from 78 recordings of 26 stations; for the Mw=6.2 aftershock of September 25, there were 81 pairs of phases clear enough from 72 recordings of 24 stations. For each station, 2 types of STFs were retrieved, which are called P-STF and S-STF due to being from P and S phases, respectively. Totally, 178 STF individuals were obtained for source-process analysis of the main shock. It was noticed that, in general, STFs from most of the stations had similarities except that those in special azimuths looked different or odd due to the mechanism difference between the main shock and the aftershocks; and in detail, the shapes of the STFs varied with azimuth. Both of them reflected the stability and reliability of the retrieved STFs. The comprehensive analysis of those STFs suggested that this event consisted of two sub-events, the total duration time was about 26 s, and on the average, the second event was about 7 s later than the first one, and the moment-rate amplitude of the first event was about 15% larger than that of the second one.
基金sponsored by the Annual Earthquake Tracking Task,CEA(2017010214)
文摘Based on the existing continuous borehole strain observation,the multiquadric function fitting method was used to deal with time series data. The impact of difference kernel function parameters was discussed to obtain a valuable fitting result,from which the physical connotation of the original data and its possible applications were analyzed.Meanwhile,a brief comparison was made between the results of multiquadric function fitting and polynomial fitting.
基金supported by the National Natural Science Foundation of China(11201324)the Fok Ying Tuny Education Foundation(141114)the Sichuan Technology Program(2022ZYD0011,2022NFSC1852).
文摘For a general normed vector space,a special optimal value function called a maximal time function is considered.This covers the farthest distance function as a special case,and has a close relationship with the smallest enclosing ball problem.Some properties of the maximal time function are proven,including the convexity,the lower semicontinuity,and the exact characterizations of its subdifferential formulas.
文摘The Region-Time-Length Algorithm (RTL algorithm) is introduced and improved in the paper. Compared with the original definition, the influence of rupture length on RTL function is emphasized and the weights of epicentral distance function, time function, and rupture length function are ensured to be equal. The retrospective examinations of RTL algorithm in North China have indicated that the anomalies obtained by the improved RTL algorithm show the short or intermediate-short term precursory features in most cases. There are two types of RTL anomalous patterns before the main shock. For the I-type, the variation pattern of the VRTL, numerical values of the VRTL(x, y, z, t) function, is complete and most of them have shown the changing pattern of rising from 0 turning dropping or dropping from 0 turning rising. For the II-type, the variation pattern of VRTL is not complete, which increases or decreases quickly from 0 and there is no evident turning, the main shock generally occurs in the short period around the peak VRTL. The rising of VRTL indicates an increase of seismic activity relative to the background level, which means the enhancement of seismic activity, while the dropping of VRTL indicates the decrease of seismic activity relative to the background level, which represents the seismic quiescence to a certain extent. According to statistical examination results of RTL algorithm in North China, the methods to distinguish the intermediate and short-term anomalies and to estimate the occurrence time of the coming main shock are given in the paper. For both I and II-type RTL anomalies, the R-value, i.e., the forecasting score, is about 0.6 and 0.3 for the 3 months forecasting period and about 0.7 and 0.4 for the 6 months forecasting period. The preliminary discussion is also made for the influences of characteristic time-span t0, characteristic distance r0, and threshold magnitude M0 on computation of VRTL, as well as some other significant problems in application.
文摘The lumped time distribution functions were proposed, which can be used for describing the dynamicsystems with two or more than two states of the end of growing polymer chain during chain addition polymerization.Numerical analysis of the lumped time distribution functions was carried out. The method for calculating molecularweight distribution of polymer in the stable free radical polymerization and more general cases was developed basedon the lumped time distribution functions.
基金supported by the National Natural Science Foundation of China (Nos.41004020 and IS201102643)
文摘This paper applies a new formulation to do moment tensor inversion for earthquakes in the Kushiro area of Japan. Comparing with conventional moment tensor inversion method, the new one takes the effect of source time function into consideration. For the inversion, best solution is obtained by minimizing the difference between the observed seismograms and the synthetic ones. And the best-fitting focal depth is determined from the variance reduction. The results indicate that half duration of source time function is proportional to the magnitude of earthquakes. Large earthquakes have long half duration, whereas that of moderate-small earthquakes is comparatively shorter. The focal mechanisms of all three earthquakes are of thrust fault type, which is mainly ascribed to the collision of the North American plate with the Eurasia plate in the late Cretaceous or Paleogene.
基金supported by the National Natural Science Foundation of China (60974139)
文摘Finite time stability and stabilization are studied for hy-brid dynamic systems. By combining multiple Lyapunov function and finite time Lyapunov function, a sufficient condition of finite time stability is given for the system. Compared with the previ-ous works, our results have less conservativeness. Furthermore, based on the state partition of continuous and resetting parts of system, a hybrid feedback controller is constructed, which stabi-lizes the closed-loop systems in finite time. Finally, a numerical example is provided to demonstrate the effectiveness of the pro-posed method.
文摘Escape time algorithm is an effective theoretical algorithm of constructing fractal graphics. The key of this algorithm lies in the construction of escape time function. A new escape time function is presented based on the research of escape time algorithm. An accelerated escape time algorithm is carried out in this paper. The experiments have demonstrated that the new algorithm is not only as precise as the old, but also faster when it is used to construct Julia set.
文摘In this study,application of the spectral representation method for generation of endurance time excitation functions is introduced.Using this method,the intensifying acceleration time series is generated so that its acceleration response spectrum in any desired time duration is compatible with a time-scaled predefined acceleration response spectrum.For this purpose,simulated stationary acceleration time series is multiplied by the time dependent linear modulation function,then using a simple iterative scheme,it is forced to match a target acceleration response spectrum.It is shown that the generated samples have excellent conformity in low frequency,which is useful for nonlinear endurance time analysis.In the second part of this study,it is shown that this procedure can be extended to generate a set of spatially correlated endurance time excitation functions.This makes it possible to assess the performance of long structures under multi-support seismic excitation using endurance time analysis.
基金sponsored under the keyresearch project of social development of Zhejiang Province(2005C23075)
文摘It has been a period of time since the concept of scenario earthquake was proposed, but this concept has rarely been used in seismic safety evaluation in China since then. Meanwhile, because of the uncertainties of magnitudes-distances pairs, there is large arbitrariness while determining the envelope function of time histories in seismic hazard analysis. In this paper, we describe a method to control the envelope functions of the time histories by introducing the most-likely combinations of magnitude and distance of the scenario earthquakes based on a probabilistic method, revise the software of the ellipse model for seismic hazard analysis, and give a computation example.
文摘Several studies showed that the breast cancer incidence rates are higher in high-income (developed) countries, due to the link of breast cancer with several risk factors and the presence of systematic screening policies. Some of the authors suggest that lower breast cancer incidence rates in low-income (developing) countries probably reflect international variation in hormonal factors and accessibility to early detection facilities. Recent studies showed that the breast cancer increased rapidly among women in Pakistan (a developing country) and it became the first malignancy among females of Pakistan. Although, the incidence rates may contain important evidence for understanding and control of the disease;however in Pakistan, the breast cancer incidence data have never been available in the last five decades since independence;rather, only hospital-based data are available. In this study, we intend to apply Functional Time Series (FTS) models to the breast cancer incidence rates of United State (developed country), and to see the difference between various components (age and time) of Functional Time Series (FTS) models applied independently on the breast cancer incidence rates of Karachi (Pakistan) and US. Past studies have already suggested that the incidence of US breast cancer cases was expected to increase in the coming decades. A progressive increase in the number of new cases is already predetermined by the high birth rate that occurred during the middle part of the century, and it will lead to nearly a doubling in the number of cases in about 4 decades. We also obtain 15 years predictions of breast cancer incidence rates in United States and compare them with the forecasts of incidence curves for Karachi. Development of methods for cancer incidence trend forecasting can provide a sound and accurate foundation for planning a comprehensive national strategy for optimal partitioning of research resources between the need for development of new treatments and the need for new research directed toward primary preventive measures.
基金supported by the National Natural Science Foundation of China(No.22273005 and No.22422301)the Innovation Program for Quantum Science and Technology(No.2023ZD0300200)+1 种基金the National Security Academic Foundation(No.U2330201)the Fundamental Research Funds for the Central Universities.
文摘Static disorder plays a crucial role in the electronic dynamics and spec-troscopy of complex molecular sys-tems.Traditionally,obtaining ob-servables averaged over static disor-der requires thousands of realiza-tions via direct sampling of the dis-order distribution,leading to high computational costs.In this work,we extend the auxiliary degree-of-freedom based matrix product state(MPS)method to handle system-bath correlated thermal equilibrium initial states,which can capture static disorder effects using a one-shot quantum dynamical simulation.We validate the effectiveness of the extended method by computing the dipole-dipole time correlation function of the Holstein model relevant to the emission spectrum of molecular aggregates.Our results show that the one-shot method is very accu-rate with only a moderate increase in MPS bond dimension,thereby significantly reducing computational cost.Moreover,it enables the generation of a much larger number of samples than the conventional direct sampling method at negligible additional cost,thus reducing sta-tistical errors.This method provides a broadly useful tool for calculating equilibrium time cor-relation functions in system-bath coupled models with static disorder.
基金supported by the Fundamental Research Funds for the Central Universities(Grant Nos.GK202207012 and QCYRCXM-2022-241).
文摘Permanent dipole moments induced high-order harmonic generation(HHG)signals offer a potential approach to producing elliptically or even circularly polarized X-ray attosecond sources.Previous studies on this topic have mainly focused on diatomic molecules such as CO and HeH.Based on this scheme,significant HHG signals in the direction perpendicular to the molecular axis can be observed in both the high-energy and low-energy regions.However,we found that the high-order harmonics induced by the permanent dipole moments of polyatomic complex molecules involve more intricate physical processes.Using time-dependent density functional theory,we simulated the dynamics of HHG from NH2COOH and NH2COSH interacting with linearly polarized lasers.We found that the harmonic signals in the direction perpendicular to the N-C bond were significantly enhanced in the high-energy photon region.Our analysis indicates that this is due to the complex molecular configuration of NH_(2)COOH and NH_(2)COSH:while the NH_(2) group has C_(2v) symmetry,both COOH and COSH groups lack this symmetry.This structural characteristic results in permanent dipole moments being felt only when electrons return to either COSH or COOH groups,but not to NH_(2) group.Additionally,our results reveal a multi-plateau structure in HHG signal along laser polarization direction,a phenomenon arising from multi-electron and multiorbital effects during interaction between complex molecule and strong laser field.
