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.展开更多
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.展开更多
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 effect of source size and emission time on the proton-proton(p-p)momentum correlation function(Cpp(q))has been studied systematically.Assuming a spherical Gaussian source with space and time profile according to t...The effect of source size and emission time on the proton-proton(p-p)momentum correlation function(Cpp(q))has been studied systematically.Assuming a spherical Gaussian source with space and time profile according to the function S(r,t)~exp(-r2/2 r02-t/τ)in the correlation function calculation code(CRAB),the results indicate that one Cpp(q)distribution corresponds to a unique combination of source size r0 and emission timeτ.Considering the possible nuclear deformation from a spherical nucleus,an ellipsoidal Gaussian source characterized by the deformation parameter∈=ΔR/R has been simulated.There is almost no difference of Cpp(q)between the results of spherically and ellipsoidally shaped sources with small deformation.These results indicate that a unique source size r0 and emission time could be extracted from the p-p momentum correlation function,which is especially important for identifying the mechanism of twoproton emission from proton-rich nuclei.Furthermore,considering the possible existence of cluster structures within a nucleus,the double Gaussian source is assumed.The results show that the p-p momentum correlation function for a source with or without cluster structures has large systematical differences with the variance of r0 andτ.This may provide a possible method for experimentally observing the cluster structures in proton-rich nuclei.展开更多
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.展开更多
The distribution of continuous service time in call centers is investigated.A non-Maxwellian collision kernel combining two different value functions in the interaction rule are used to describe the evolution of conti...The distribution of continuous service time in call centers is investigated.A non-Maxwellian collision kernel combining two different value functions in the interaction rule are used to describe the evolution of continuous service time,respectively.Using the statistical mechanical and asymptotic limit methods,Fokker–Planck equations are derived from the corresponding Boltzmann-type equations with non-Maxwellian collision kernels.The steady-state solutions of the Fokker–Planck equation are obtained in exact form.Numerical experiments are provided to support our results under different parameters.展开更多
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.展开更多
Time synchronization(TS)is crucial for ensuring the secure and reliable functioning of the distribution power Internet of Things(IoT).Multi-clock source time synchronization(MTS)has significant advantages of high reli...Time synchronization(TS)is crucial for ensuring the secure and reliable functioning of the distribution power Internet of Things(IoT).Multi-clock source time synchronization(MTS)has significant advantages of high reliability and accuracy but still faces challenges such as optimization of the multi-clock source selection and the clock source weight calculation at different timescales,and the coupling of synchronization latency jitter and pulse phase difference.In this paper,the multi-timescale MTS model is conducted,and the reinforcement learning(RL)and analytic hierarchy process(AHP)-based multi-timescale MTS algorithm is designed to improve the weighted summation of synchronization latency jitter standard deviation and average pulse phase difference.Specifically,the multi-clock source selection is optimized based on Softmax in the large timescale,and the clock source weight calculation is optimized based on lower confidence bound-assisted AHP in the small timescale.Simulation shows that the proposed algorithm can effectively reduce time synchronization delay standard deviation and average pulse phase difference.展开更多
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.展开更多
Nonlinear variations in the coordinate time series of global navigation satellite system(GNSS) reference stations are strongly correlated with surface displacements caused by environmental loading effects,including at...Nonlinear variations in the coordinate time series of global navigation satellite system(GNSS) reference stations are strongly correlated with surface displacements caused by environmental loading effects,including atmospheric, hydrological, and nontidal ocean loading. Continuous improvements in the accuracy of surface mass loading products, performance of Earth models, and precise data-processing technologies have significantly advanced research on the effects of environmental loading on nonlinear variations in GNSS coordinate time series. However, owing to theoretical limitations, the lack of high spatiotemporal resolution surface mass observations, and the coupling of GNSS technology-related systematic errors, environmental loading and nonlinear GNSS reference station displacements remain inconsistent. The applicability and capability of these loading products across different regions also require further evaluation. This paper outlines methods for modeling environmental loading, surface mass loading products, and service organizations. In addition, it summarizes recent advances in applying environmental loading to address nonlinear variations in global and regional GNSS coordinate time series. Moreover, the scientific questions of existing studies are summarized, and insights into future research directions are provided. The complex nonlinear motion of reference stations is a major factor limiting the accuracy of the current terrestrial reference frame. Further refining the environmental load modeling method, establishing a surface mass distribution model with high spatiotemporal resolution and reliability, exploring other environmental load factors such as ice sheet and artificial mass-change effects, and developing an optimal data-processing model and strategy for reprocessing global reference station data consistently could contribute to the development of a millimeter-level nonlinear motion model for GNSS reference stations with actual physical significance and provide theoretical support for establishing a terrestrial reference frame with 1 mm accuracy by 2050.展开更多
The acquisition of neutron time spectrum data plays a pivotal role in the precise quantification of uranium via prompt fission neutron uranium logging(PFNUL).However,the impact of the detector dead-time effect remains...The acquisition of neutron time spectrum data plays a pivotal role in the precise quantification of uranium via prompt fission neutron uranium logging(PFNUL).However,the impact of the detector dead-time effect remains paramount in the accurate acquisition of the neutron time spectrum.Therefore,it is imperative for neutron logging instruments to establish a dead-time correction method that is not only uncomplicated but also practical and caters to various logging sites.This study has formulated an innovative equation for determining dead time and introduced a dead-time correction method for the neutron time spectrum,called the“dual flux method.”Using this approach,a logging instrument captures two neutron time spectra under disparate neutron fluxes.By carefully selecting specific“windows”on the neutron time spectrum,the dead time can be accurately ascertained.To substantiate its efficacy and discern the influencing factors,experiments were conducted utilizing a deuterium-tritium(D-T)neutron source,a Helium-3(3He)detector,and polyethylene shielding to collate and analyze the neutron time spectrum under varying neutron fluxes(at high voltages).The findings underscore that the“height”and“spacing”of the two windows are the most pivotal influencing factors.Notably,the“height”(fd)should surpass 2,and the“spacing”twd should exceed 200μs.The dead time of the 3 He detector determined in the experiment was 7.35μs.After the dead-time correction,the deviation of the decay coefficients from the theoretical values for the neutron time spectrum under varying neutron fluxes decreased from 12.4%to within 5%.Similarly,for the PFNUL instrument,the deviation in the decay coefficients decreased from 22.94 to 0.49%after correcting for the dead-time effect.These results demonstrate the exceptional efficacy of the proposed method in ensuring precise uranium quantification.The dual flux method was experimentally validated as a universal approach applicable to pulsed neutron logging instruments and holds immense significance for uranium exploration.展开更多
A source-time function plays an important role in generating finite-difference synthetic seismograms.In this regard,we consider the amplitude spectra of five different time functions.This investigation shows that when...A source-time function plays an important role in generating finite-difference synthetic seismograms.In this regard,we consider the amplitude spectra of five different time functions.This investigation shows that when the high-frequency signals of the time function carry more power than the main signal,the tailing phenomena would be associated with the main signal which is not desirable in producing finite-difference synthetic seismograms.It is illustrated that the time function denoted by is the optimal source-time function among those examined in this paper for generating finite-difference synthetic seismograms.This matter is also verified by making finite-difference synthetic seismograms over an elastic earth model consisting of a layer over a half space.展开更多
Background:Rehabilitation typically restores shoulder function after a common rotator cuff repair;however,it is unclear when to initiate physical therapy(PT)post-surgery.To understand effective start times,this study ...Background:Rehabilitation typically restores shoulder function after a common rotator cuff repair;however,it is unclear when to initiate physical therapy(PT)post-surgery.To understand effective start times,this study analyzed PT start times and their effect on shoulder function.The most effective shoulder outcomes were hypothesized to occur when PT started 0-4 weeks post-surgery.Methods:Data from 1002 patients who underwent rotator cuff repair were analyzed retrospectively from 52 outpatient clinics in the years of 2016-2021.The primary data of timings between surgery and the start of PT(0-2 weeks[0-14 days],2-4 weeks[15-28 days],4-6 weeks[29-42 days],and 6-14 weeks[43-100 days]),change in functional shoulder scores,number of PT visits utilized,and functional scores changed per visit were analyzed.Regression analyses of the impact of age,sex,payor source,and tear size were completed.Results:The change in functional scores and functional scores changed per visit are not affected by the initiation timing of PT postrotator cuff surgery,when controlling for baseline functional scores.This result is revealed despite initial functional scores being significantly different.Functional scores change per visit may have been affected by payor source(R^(2)=0.004,p≤0.010).PT start times may have been impacted by age(R^(2)=0.010,p=0.007),payor source(R^(2)=0.016,p=0.001),and tear size(R^(2)=0.007,p=0.026).Payor source may have influenced the number of PT visits(R^(2)=0.060,p<0.001).Conclusions:After rotator cuff surgery,patients should choose to complete rehabilitation to optimize shoulder function.However,the initiation timing of PT may not affect functional shoulder outcomes.展开更多
A predator-prey model with prey dispersal and Holling type-Ⅱ functional response is investigated.In this model,the time delay due to the gestation of the predator and stagestructure for the predator are considered.By...A predator-prey model with prey dispersal and Holling type-Ⅱ functional response is investigated.In this model,the time delay due to the gestation of the predator and stagestructure for the predator are considered.By analyzing the corresponding characteristic equations,the local stability of each of the nonnegative equilibria is discussed.The existence of Hopf bifurcations at the positive equilibrium is established.By using Lyapunov functionals and LaSalle’s invariance principle,sufficient conditions are obtained for the global stability of the positive equilibrium,the nonnegative boundary equilibrium and the trivial equilibrium of the model,respectively.Numerical simulations are carried out to illustrate the main results.展开更多
systematic verification and validation(V&V)of our previously proposed momentum source wave generation method is performed.Some settings of previous numerical wave tanks(NWTs)of regular and irregular waves have bee...systematic verification and validation(V&V)of our previously proposed momentum source wave generation method is performed.Some settings of previous numerical wave tanks(NWTs)of regular and irregular waves have been optimized.The H2-5 V&V method involving five mesh sizes with mesh refinement ratio being 1.225 is used to verify the NWT of regular waves,in which the wave height and mass conservation are mainly considered based on a Lv3(H s=0.75 m)and a Lv6(H s=5 m)regular wave.Additionally,eight different sea states are chosen to validate the wave height,mass conservation and wave frequency of regular waves.Regarding the NWT of irregular waves,five different sea states with significant wave heights ranging from 0.09 m to 12.5 m are selected to validate the statistical characteristics of irregular waves,including the profile of the wave spectrum,peak frequency and significant wave height.Results show that the verification errors for Lv3 and Lv6 regular wave on the most refined grid are−0.018 and−0.35 for wave height,respectively,and−0.14 and for−0.17 mass conservation,respectively.The uncertainty estimation analysis shows that the numerical error could be partially balanced out by the modelling error to achieve a smaller validation error by adjusting the mesh size elaborately.And the validation errors of the wave height,mass conservation and dominant frequency of regular waves under different sea states are no more than 7%,8% and 2%,respectively.For a Lv3(H_(s)=0.75 m)and a Lv6(H_(s)=5 m)regular wave,simulations are validated on the wave height in wave development section for safety factors FS≈1 and FS≈0.5-1,respectively.Regarding irregular waves,the validation errors of the significant wave height and peak frequency are both lower than 2%.展开更多
Unveiling the underlying physical mechanisms governing the fracture of brittle rocks is imperative for preventing rockbursts.The novelty of this study lies in the analysis of the dynamic response process of rock three...Unveiling the underlying physical mechanisms governing the fracture of brittle rocks is imperative for preventing rockbursts.The novelty of this study lies in the analysis of the dynamic response process of rock three-dimensional(3D)deformation under true triaxial stress,and the surge behavior of timedependent multifractal spectrum has been successfully used to warn of progressive failure inside the rock.Firstly,this study analyzed the dynamic adjustment trajectory of rock deformation,specifically lateral strain,within the framework of the Poisson effect.This analysis highlighted the intricate dependence of rock mechanical properties on the intermediate principal stress.Secondly,by defining the crack interval function(ICF),this study compared the disparities between the two crack growth stages(strengthening stage and weakening stage)under varying stress levels.It was found that the fracture activity of granite system has significant multifractal characteristics.Notably,the multifractal spectrum emerges as a valuable tool for characterizing the distinct fracture properties of rocks,encompassing both the crack scale and the associated energy.Finally,a quantitative criterion grounded in the multifractal parameters of the acoustic emission(AE)time series was formulated,and it indicates that the abrupt changes observed in the time-dependent fractal spectra can serve as precursor indicators for the progressive development of rockbursts.展开更多
Black carbon (BC) aerosols are considered key factors that contribute to rapid climate warming and ice melt in the Arctic region.However,compared with long-term observations from land-based stations,observational data...Black carbon (BC) aerosols are considered key factors that contribute to rapid climate warming and ice melt in the Arctic region.However,compared with long-term observations from land-based stations,observational data over the Arctic Ocean remain relatively scarce.Four Arctic scientific expeditions were conducted in the summer and early autumn of 2010,2012,2016,and 2018 via the Chinese research vessel Xue Long,during which the BC concentrations along the routes were measured via light absorption methods.In this work,the spatiotemporal distribution characteristics of BC over the Arctic Ocean were examined on the basis of these observations.The potential sources of BC along the various routes were analyzed via the weighted potential source contribution function and weighted concentrationweighted trajectory methods of the hybrid single-particle Lagrangian integrated trajectory model in conjunction with Arctic transport potential climate model simulations.The analysis results indicated that wildfires in the western Aleutian Islands,Siberia,and Far East regions were the primary contributors to the BC aerosol concentration observed along the Arctic expedition routes in summer,identifying these regions as major potential source areas.展开更多
Dear Editor,This letter investigates predefined-time optimization problems(OPs) of multi-agent systems(MASs), where the agent of MASs is subject to inequality constraints, and the team objective function accounts for ...Dear Editor,This letter investigates predefined-time optimization problems(OPs) of multi-agent systems(MASs), where the agent of MASs is subject to inequality constraints, and the team objective function accounts for impulse effects. Firstly, to address the inequality constraints,the penalty method is introduced. Then, a novel optimization strategy is developed, which only requires that the team objective function be strongly convex.展开更多
Ultrafast electron sources, which enable high spatiotemporal resolution in time-resolved electron microscopy and scanning probe microscopy, are receiving increased attention. The most widely used method for achieving ...Ultrafast electron sources, which enable high spatiotemporal resolution in time-resolved electron microscopy and scanning probe microscopy, are receiving increased attention. The most widely used method for achieving ultrafast electron sources involves irradiating metal tips by ultrashort laser pulses, causing electron beam emission via the photoelectric effect [including photon-driven(quantum) or field-driven(classical) emission]. However, the thermionic electrons emission process due to the heating effect of ultrashort lasers, particularly its dynamic aspects, has rarely been addressed in previous studies. In this paper, we improved the signal-to-noise ratio of a two-pulse correlation measurement on the tip electron emission by nearly two orders of magnitude using a delay time modulation method. This allowed us to obtain information on the temperature evolution of hot electrons and phonons in a non-equilibrium state, and to extract characteristic time scales for electron-phonon and phonon-phonon scattering. Our findings indicate that the thermionic electrons emission, unlike the instantaneous photoelectric effect, causes electron emission to lag behind the laser pulse by tens of picoseconds, thus significantly affecting the detection of ultrafast dynamics of samples. Furthermore, such a lagging effect was found to be sensitive to the local structure of the metal tip, offering new insights into the improved design of ultrafast electron sources.展开更多
With Biolog Eco microplate, metabolic characteristics and functional diver-sity of carbon source in microflora of ponds were researched based on recitculating ponds and control ponds in order to explore effects of eco...With Biolog Eco microplate, metabolic characteristics and functional diver-sity of carbon source in microflora of ponds were researched based on recitculating ponds and control ponds in order to explore effects of eco-adjustments on microflo-ra in ponds. The results indicate that total number of bacterium, microbial metabolism activity, and diversity index in P7, P8, P1 and P2 kept higher, fol owed by P3, P4, P5 and P6. The utilization rate of microbes on sugars achieved the highest （31.0%-48.7%）, fol owed by carboxylic acid （13.4%-18.0%）, amino acid （10.1%-20.5%）, polymers （9.4%-17.0%）, biopolymer （5.7%-9.7%） and phenol （4.95%-7.50%）. Principal component analysis divided microflora in different ponds, suggesting that microbial community has varied carbon source characteristics and nitrogen-containing compound and biopolymer metabolisms are most affected.展开更多
基金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 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.
文摘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 National Key R&D Program of China(No.2018YFA0404404)the National Natural Science Foundation of China(Nos.11925502,11935001,11961141003,11421505,11475244,and 11927901)+2 种基金the Shanghai Development Foundation for Science and Technology(No.19ZR1403100)the Strategic Priority Research Program of the CAS(No.XDB34030000)the Key Research Program of Frontier Sciences of the CAS(No.QYZDJ-SSW-SLH002).
文摘The effect of source size and emission time on the proton-proton(p-p)momentum correlation function(Cpp(q))has been studied systematically.Assuming a spherical Gaussian source with space and time profile according to the function S(r,t)~exp(-r2/2 r02-t/τ)in the correlation function calculation code(CRAB),the results indicate that one Cpp(q)distribution corresponds to a unique combination of source size r0 and emission timeτ.Considering the possible nuclear deformation from a spherical nucleus,an ellipsoidal Gaussian source characterized by the deformation parameter∈=ΔR/R has been simulated.There is almost no difference of Cpp(q)between the results of spherically and ellipsoidally shaped sources with small deformation.These results indicate that a unique source size r0 and emission time could be extracted from the p-p momentum correlation function,which is especially important for identifying the mechanism of twoproton emission from proton-rich nuclei.Furthermore,considering the possible existence of cluster structures within a nucleus,the double Gaussian source is assumed.The results show that the p-p momentum correlation function for a source with or without cluster structures has large systematical differences with the variance of r0 andτ.This may provide a possible method for experimentally observing the cluster structures in proton-rich nuclei.
基金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.
基金the Special Project of Yili Normal University(to improve comprehensive strength of disciplines)(Grant No.22XKZZ18)Yili Normal University Scientific Research Innovation Team Plan Project(Grant No.CXZK2021015)Yili Science and Technology Planning Project(Grant No.YZ2022B036).
文摘The distribution of continuous service time in call centers is investigated.A non-Maxwellian collision kernel combining two different value functions in the interaction rule are used to describe the evolution of continuous service time,respectively.Using the statistical mechanical and asymptotic limit methods,Fokker–Planck equations are derived from the corresponding Boltzmann-type equations with non-Maxwellian collision kernels.The steady-state solutions of the Fokker–Planck equation are obtained in exact form.Numerical experiments are provided to support our results under different parameters.
基金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.
基金supported by Science and Technology Project of China Southern Power Grid Company Limited under Grant Number 036000KK52200058(GDKJXM20202001).
文摘Time synchronization(TS)is crucial for ensuring the secure and reliable functioning of the distribution power Internet of Things(IoT).Multi-clock source time synchronization(MTS)has significant advantages of high reliability and accuracy but still faces challenges such as optimization of the multi-clock source selection and the clock source weight calculation at different timescales,and the coupling of synchronization latency jitter and pulse phase difference.In this paper,the multi-timescale MTS model is conducted,and the reinforcement learning(RL)and analytic hierarchy process(AHP)-based multi-timescale MTS algorithm is designed to improve the weighted summation of synchronization latency jitter standard deviation and average pulse phase difference.Specifically,the multi-clock source selection is optimized based on Softmax in the large timescale,and the clock source weight calculation is optimized based on lower confidence bound-assisted AHP in the small timescale.Simulation shows that the proposed algorithm can effectively reduce time synchronization delay standard deviation and average pulse phase difference.
基金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 Basic Science Center Project of the National Natural Science Foundation of China(42388102)the National Natural Science Foundation of China(42174030)+2 种基金the Special Fund of Hubei Luojia Laboratory(220100020)the Major Science and Technology Program for Hubei Province(2022AAA002)the Fundamental Research Funds for the Central Universities of China(2042022dx0001 and 2042023kfyq01)。
文摘Nonlinear variations in the coordinate time series of global navigation satellite system(GNSS) reference stations are strongly correlated with surface displacements caused by environmental loading effects,including atmospheric, hydrological, and nontidal ocean loading. Continuous improvements in the accuracy of surface mass loading products, performance of Earth models, and precise data-processing technologies have significantly advanced research on the effects of environmental loading on nonlinear variations in GNSS coordinate time series. However, owing to theoretical limitations, the lack of high spatiotemporal resolution surface mass observations, and the coupling of GNSS technology-related systematic errors, environmental loading and nonlinear GNSS reference station displacements remain inconsistent. The applicability and capability of these loading products across different regions also require further evaluation. This paper outlines methods for modeling environmental loading, surface mass loading products, and service organizations. In addition, it summarizes recent advances in applying environmental loading to address nonlinear variations in global and regional GNSS coordinate time series. Moreover, the scientific questions of existing studies are summarized, and insights into future research directions are provided. The complex nonlinear motion of reference stations is a major factor limiting the accuracy of the current terrestrial reference frame. Further refining the environmental load modeling method, establishing a surface mass distribution model with high spatiotemporal resolution and reliability, exploring other environmental load factors such as ice sheet and artificial mass-change effects, and developing an optimal data-processing model and strategy for reprocessing global reference station data consistently could contribute to the development of a millimeter-level nonlinear motion model for GNSS reference stations with actual physical significance and provide theoretical support for establishing a terrestrial reference frame with 1 mm accuracy by 2050.
基金supported by the National Natural Science Foundation of China(No.42374226)Jiangxi Provincial Natural Science Foundation(Nos.20232BAB201043 and 20232BCJ23006)+2 种基金Nuclear Energy Development Project(20201192-01)National Key Laboratory of Uranium Resource Exploration-Mining and Nuclear Remote Sensing(ECUT)(2024QZ-TD-09)Fundamental Science on Radioactive Geology and Exploration Technology Laboratory(2022RGET20).
文摘The acquisition of neutron time spectrum data plays a pivotal role in the precise quantification of uranium via prompt fission neutron uranium logging(PFNUL).However,the impact of the detector dead-time effect remains paramount in the accurate acquisition of the neutron time spectrum.Therefore,it is imperative for neutron logging instruments to establish a dead-time correction method that is not only uncomplicated but also practical and caters to various logging sites.This study has formulated an innovative equation for determining dead time and introduced a dead-time correction method for the neutron time spectrum,called the“dual flux method.”Using this approach,a logging instrument captures two neutron time spectra under disparate neutron fluxes.By carefully selecting specific“windows”on the neutron time spectrum,the dead time can be accurately ascertained.To substantiate its efficacy and discern the influencing factors,experiments were conducted utilizing a deuterium-tritium(D-T)neutron source,a Helium-3(3He)detector,and polyethylene shielding to collate and analyze the neutron time spectrum under varying neutron fluxes(at high voltages).The findings underscore that the“height”and“spacing”of the two windows are the most pivotal influencing factors.Notably,the“height”(fd)should surpass 2,and the“spacing”twd should exceed 200μs.The dead time of the 3 He detector determined in the experiment was 7.35μs.After the dead-time correction,the deviation of the decay coefficients from the theoretical values for the neutron time spectrum under varying neutron fluxes decreased from 12.4%to within 5%.Similarly,for the PFNUL instrument,the deviation in the decay coefficients decreased from 22.94 to 0.49%after correcting for the dead-time effect.These results demonstrate the exceptional efficacy of the proposed method in ensuring precise uranium quantification.The dual flux method was experimentally validated as a universal approach applicable to pulsed neutron logging instruments and holds immense significance for uranium exploration.
文摘A source-time function plays an important role in generating finite-difference synthetic seismograms.In this regard,we consider the amplitude spectra of five different time functions.This investigation shows that when the high-frequency signals of the time function carry more power than the main signal,the tailing phenomena would be associated with the main signal which is not desirable in producing finite-difference synthetic seismograms.It is illustrated that the time function denoted by is the optimal source-time function among those examined in this paper for generating finite-difference synthetic seismograms.This matter is also verified by making finite-difference synthetic seismograms over an elastic earth model consisting of a layer over a half space.
文摘Background:Rehabilitation typically restores shoulder function after a common rotator cuff repair;however,it is unclear when to initiate physical therapy(PT)post-surgery.To understand effective start times,this study analyzed PT start times and their effect on shoulder function.The most effective shoulder outcomes were hypothesized to occur when PT started 0-4 weeks post-surgery.Methods:Data from 1002 patients who underwent rotator cuff repair were analyzed retrospectively from 52 outpatient clinics in the years of 2016-2021.The primary data of timings between surgery and the start of PT(0-2 weeks[0-14 days],2-4 weeks[15-28 days],4-6 weeks[29-42 days],and 6-14 weeks[43-100 days]),change in functional shoulder scores,number of PT visits utilized,and functional scores changed per visit were analyzed.Regression analyses of the impact of age,sex,payor source,and tear size were completed.Results:The change in functional scores and functional scores changed per visit are not affected by the initiation timing of PT postrotator cuff surgery,when controlling for baseline functional scores.This result is revealed despite initial functional scores being significantly different.Functional scores change per visit may have been affected by payor source(R^(2)=0.004,p≤0.010).PT start times may have been impacted by age(R^(2)=0.010,p=0.007),payor source(R^(2)=0.016,p=0.001),and tear size(R^(2)=0.007,p=0.026).Payor source may have influenced the number of PT visits(R^(2)=0.060,p<0.001).Conclusions:After rotator cuff surgery,patients should choose to complete rehabilitation to optimize shoulder function.However,the initiation timing of PT may not affect functional shoulder outcomes.
基金Supported by the Social Science Foundation of Hebei Province(HB23TJO03)。
文摘A predator-prey model with prey dispersal and Holling type-Ⅱ functional response is investigated.In this model,the time delay due to the gestation of the predator and stagestructure for the predator are considered.By analyzing the corresponding characteristic equations,the local stability of each of the nonnegative equilibria is discussed.The existence of Hopf bifurcations at the positive equilibrium is established.By using Lyapunov functionals and LaSalle’s invariance principle,sufficient conditions are obtained for the global stability of the positive equilibrium,the nonnegative boundary equilibrium and the trivial equilibrium of the model,respectively.Numerical simulations are carried out to illustrate the main results.
基金supported by the National Key R&D Program of China(Grant No.2022YFB3303500).
文摘systematic verification and validation(V&V)of our previously proposed momentum source wave generation method is performed.Some settings of previous numerical wave tanks(NWTs)of regular and irregular waves have been optimized.The H2-5 V&V method involving five mesh sizes with mesh refinement ratio being 1.225 is used to verify the NWT of regular waves,in which the wave height and mass conservation are mainly considered based on a Lv3(H s=0.75 m)and a Lv6(H s=5 m)regular wave.Additionally,eight different sea states are chosen to validate the wave height,mass conservation and wave frequency of regular waves.Regarding the NWT of irregular waves,five different sea states with significant wave heights ranging from 0.09 m to 12.5 m are selected to validate the statistical characteristics of irregular waves,including the profile of the wave spectrum,peak frequency and significant wave height.Results show that the verification errors for Lv3 and Lv6 regular wave on the most refined grid are−0.018 and−0.35 for wave height,respectively,and−0.14 and for−0.17 mass conservation,respectively.The uncertainty estimation analysis shows that the numerical error could be partially balanced out by the modelling error to achieve a smaller validation error by adjusting the mesh size elaborately.And the validation errors of the wave height,mass conservation and dominant frequency of regular waves under different sea states are no more than 7%,8% and 2%,respectively.For a Lv3(H_(s)=0.75 m)and a Lv6(H_(s)=5 m)regular wave,simulations are validated on the wave height in wave development section for safety factors FS≈1 and FS≈0.5-1,respectively.Regarding irregular waves,the validation errors of the significant wave height and peak frequency are both lower than 2%.
基金funding support from the National Natural Science Foundation of China(Grant No.U2034207)the Natural Science Foundation of Hebei Province(Grant No.E2021210099).
文摘Unveiling the underlying physical mechanisms governing the fracture of brittle rocks is imperative for preventing rockbursts.The novelty of this study lies in the analysis of the dynamic response process of rock three-dimensional(3D)deformation under true triaxial stress,and the surge behavior of timedependent multifractal spectrum has been successfully used to warn of progressive failure inside the rock.Firstly,this study analyzed the dynamic adjustment trajectory of rock deformation,specifically lateral strain,within the framework of the Poisson effect.This analysis highlighted the intricate dependence of rock mechanical properties on the intermediate principal stress.Secondly,by defining the crack interval function(ICF),this study compared the disparities between the two crack growth stages(strengthening stage and weakening stage)under varying stress levels.It was found that the fracture activity of granite system has significant multifractal characteristics.Notably,the multifractal spectrum emerges as a valuable tool for characterizing the distinct fracture properties of rocks,encompassing both the crack scale and the associated energy.Finally,a quantitative criterion grounded in the multifractal parameters of the acoustic emission(AE)time series was formulated,and it indicates that the abrupt changes observed in the time-dependent fractal spectra can serve as precursor indicators for the progressive development of rockbursts.
基金supported by the National Natural Science Foundation of China (No.42201151)the Basic Research Fund of the Chinese Academy of Meteorological Sciences (Nos.2023Z004 and 2024Z007)。
文摘Black carbon (BC) aerosols are considered key factors that contribute to rapid climate warming and ice melt in the Arctic region.However,compared with long-term observations from land-based stations,observational data over the Arctic Ocean remain relatively scarce.Four Arctic scientific expeditions were conducted in the summer and early autumn of 2010,2012,2016,and 2018 via the Chinese research vessel Xue Long,during which the BC concentrations along the routes were measured via light absorption methods.In this work,the spatiotemporal distribution characteristics of BC over the Arctic Ocean were examined on the basis of these observations.The potential sources of BC along the various routes were analyzed via the weighted potential source contribution function and weighted concentrationweighted trajectory methods of the hybrid single-particle Lagrangian integrated trajectory model in conjunction with Arctic transport potential climate model simulations.The analysis results indicated that wildfires in the western Aleutian Islands,Siberia,and Far East regions were the primary contributors to the BC aerosol concentration observed along the Arctic expedition routes in summer,identifying these regions as major potential source areas.
基金supported in part by the National Natural Science Foundation of China(62276119)the Natural Science Foundation of Jiangsu Province(BK20241764)the Postgraduate Research & Practice Innovation Program of Jiangsu Province(KYCX22_2860)
文摘Dear Editor,This letter investigates predefined-time optimization problems(OPs) of multi-agent systems(MASs), where the agent of MASs is subject to inequality constraints, and the team objective function accounts for impulse effects. Firstly, to address the inequality constraints,the penalty method is introduced. Then, a novel optimization strategy is developed, which only requires that the team objective function be strongly convex.
基金supported by the National Key R&D Program under Grant No.2021YFA1400500the National Natural Science Foundation of China under Grant No.22273029+1 种基金the New Cornerstone Science Foundation through the New Cornerstone Investigator Program under Grant No.NCI202303 and the XPLORER PRIZEthe Beijing Outstanding Young Scientist Program under Grant No.JWZQ20240101002。
文摘Ultrafast electron sources, which enable high spatiotemporal resolution in time-resolved electron microscopy and scanning probe microscopy, are receiving increased attention. The most widely used method for achieving ultrafast electron sources involves irradiating metal tips by ultrashort laser pulses, causing electron beam emission via the photoelectric effect [including photon-driven(quantum) or field-driven(classical) emission]. However, the thermionic electrons emission process due to the heating effect of ultrashort lasers, particularly its dynamic aspects, has rarely been addressed in previous studies. In this paper, we improved the signal-to-noise ratio of a two-pulse correlation measurement on the tip electron emission by nearly two orders of magnitude using a delay time modulation method. This allowed us to obtain information on the temperature evolution of hot electrons and phonons in a non-equilibrium state, and to extract characteristic time scales for electron-phonon and phonon-phonon scattering. Our findings indicate that the thermionic electrons emission, unlike the instantaneous photoelectric effect, causes electron emission to lag behind the laser pulse by tens of picoseconds, thus significantly affecting the detection of ultrafast dynamics of samples. Furthermore, such a lagging effect was found to be sensitive to the local structure of the metal tip, offering new insights into the improved design of ultrafast electron sources.
基金Supported by National Modern Agricultural Technology System(CARS-46)NationalSci-tech Support Plan(2012BAD25B05,2012BAD25B01)National Department PublicBenefit Research Foundation(201203083)~~
文摘With Biolog Eco microplate, metabolic characteristics and functional diver-sity of carbon source in microflora of ponds were researched based on recitculating ponds and control ponds in order to explore effects of eco-adjustments on microflo-ra in ponds. The results indicate that total number of bacterium, microbial metabolism activity, and diversity index in P7, P8, P1 and P2 kept higher, fol owed by P3, P4, P5 and P6. The utilization rate of microbes on sugars achieved the highest （31.0%-48.7%）, fol owed by carboxylic acid （13.4%-18.0%）, amino acid （10.1%-20.5%）, polymers （9.4%-17.0%）, biopolymer （5.7%-9.7%） and phenol （4.95%-7.50%）. Principal component analysis divided microflora in different ponds, suggesting that microbial community has varied carbon source characteristics and nitrogen-containing compound and biopolymer metabolisms are most affected.
