Search

Abstract Optical solitons have long been of interest both from a fundamental perspective and because of their application potential. Both cubic (Kerr) and quadratic nonlinearities can lead to soliton formation, but quadratic solitons can practically benefit from stronger nonlinearity and achieve substantial wavelength conversion. However, despite their rich physics, quadratic cavity solitons have been used only for broadband frequency comb generation, especially in the mid-infrared. Here, we show that the formation dynamics of mid-infrared quadratic cavity solitons, specifically temporal simultons in optical parametric oscillators, can be effectively leveraged to enhance molecular sensing. We demonstrate significant sensitivity enhancement while circumventing constraints of traditional cavity enhancement mechanisms. We perform experiments sensing CO2 using cavity simultons around 4 μm and achieve an enhancement of 6000. Additionally, we demonstrate large sensitivity at high concentrations of CO2, beyond what can be achieved using an equivalent high-finesse linear cavity by orders of magnitude. Our results highlight a path for utilizing quadratic cavity nonlinear dynamics and solitons for molecular sensing beyond what can be achieved using linear methods.

ABSTRACTThe increasing incidence of diseases caused by Coxsackievirus A6 (CV-A6) and the presence of various mutants in the population present significant public health challenges. Given the concurrent development of multiple vaccines in China, it is challenging to objectively and accurately evaluate the level of neutralizing antibody response to different vaccines. The choice of the detection strain is a crucial factor that influences the detection of neutralizing antibodies. In this study, the National Institutes for Food and Drug Control collected a prototype strain (Gdula), one subgenotype D1, as well as 13 CV-A6 candidate vaccine strains and candidate detection strains (subgenotype D3) from various institutions and manufacturers involved in research and development. We evaluated cross-neutralization activity using plasma from naturally infected adults (n = 30) and serum from rats immunized with the aforementioned CV-A6 strains. Although there were differences between the geometric mean titer (GMT) ranges of human plasma and murine sera, the overall trends were similar. A significant effect of each strain on the neutralizing antibody test (MAX/MIN 48.0 ∼16410.3) was observed. Among all strains, neutralization of the S112 strain by 15 different sera resulted in higher neutralizing antibody titers (GMTS112 = 132.0) and more consistent responses across different genotypic immune sera (MAX/MIN = 48.0). Therefore, S112 may serve as a detection strain for NtAb testing in various vaccines, minimizing bias and making it suitable for evaluating the immunogenicity of the CV-A6 vaccine.

SRY‐box transcription factor 6 (SOX6) is a member of the SOX gene family and inhibits the proliferation of cervical cancer cells by inducing cell cycle arrest. However, the final cell fate and significance of these cell‐cycle‐arrested cervical cancer cells induced by SOX6 remains unclear. Here, we report that SOX6 inhibits the proliferation of cervical cancer cells by inducing cellular senescence, which is mainly mediated by promoting transforming growth factor beta 2 (TGFB2) gene expression and subsequently activating the TGFβ2–Smad2/3–p53–p21WAF1/CIP1–Rb pathway. SOX6 promotes TGFB2 gene expression through the MAP4K4–MAPK (JNK/ERK/p38)–ATF2 and WT1–ATF2 pathways, which is dependent on its high‐mobility group (HMG) domain. In addition, the SOX6‐induced senescent cervical cancer cells are resistant to cisplatin treatment. ABT‐263 (navitoclax) and ABT‐199 (venetoclax), two classic senolytics, can specifically eliminate the SOX6‐induced senescent cervical cancer cells, and thus significantly improve the chemosensitivity of cisplatin‐resistant cervical cancer cells. This study uncovers that the MAP4K4/WT1–ATF2–TGFβ2 axis mediates SOX6‐induced cellular senescence, which is a promising therapeutic target in improving the chemosensitivity of cervical cancer.

Abstract Optical frequency comb is an enabling technology for a multitude of applications from metrology to ranging and communications. The tremendous progress in sources of optical frequency combs has mostly been centered around the near-infrared spectral region, while many applications demand sources in the visible and mid-infrared, which have so far been challenging to achieve, especially in nanophotonics. Here, we report widely tunable frequency comb generation using optical parametric oscillators in lithium niobate nanophotonics. We demonstrate sub-picosecond frequency combs tunable beyond an octave extending from 1.5 up to 3.3 μm with femtojoule-level thresholds on a single chip. We utilize the up-conversion of the infrared combs to generate visible frequency combs reaching 620 nm on the same chip. The ultra-broadband tunability and visible-to-mid-infrared spectral coverage of our source highlight a practical and universal path for the realization of efficient frequency comb sources in nanophotonics, overcoming their spectral sparsity.

IntroductionHIV molecular network based on genetic distance (GD) has been extensively utilized. However, the GD threshold for the non-B subtype differs from that of subtype B. This study aimed to optimize the GD threshold for inferring the CRF01_AE molecular network.MethodsNext-generation sequencing data of partial CRF01_AE pol sequences were obtained for 59 samples from 12 transmission pairs enrolled from a high-risk cohort during 2009 and 2014. The paired GD was calculated using the Tamura-Nei 93 model to infer a GD threshold range for HIV molecular networks.Results2,019 CRF01_AE pol sequences and information on recent HIV infection (RHI) from newly diagnosed individuals in Shenyang from 2016 to 2019 were collected to construct molecular networks to assess the ability of the inferred GD thresholds to predict recent transmission events. When HIV transmission occurs within a span of 1-4 years, the mean paired GD between the sequences of the donor and recipient within the same transmission pair were as follow: 0.008, 0.011, 0.013, and 0.023 substitutions/site. Using these four GD thresholds, it was found that 98.9%, 96.0%, 88.2%, and 40.4% of all randomly paired GD values from 12 transmission pairs were correctly identified as originating from the same transmission pairs. In the real world, as the GD threshold increased from 0.001 to 0.02 substitutions/site, the proportion of RHI within the molecular network gradually increased from 16.6% to 92.3%. Meanwhile, the proportion of links with RHI gradually decreased from 87.0% to 48.2%. The two curves intersected at a GD of 0.008 substitutions/site.DiscussionA suitable range of GD thresholds, 0.008-0.013 substitutions/site, was identified to infer the CRF01_AE molecular transmission network and identify HIV transmission events that occurred within the past three years. This finding provides valuable data for selecting an appropriate GD thresholds in constructing molecular networks for non-B subtypes.

Abstract With the continuous in-depth study of the interaction mechanism between viruses and hosts, the virus has become a promising tool in cancer treatment. In fact, many oncolytic viruses with selectivity and effectiveness have been used in cancer therapy. Human enterovirus is one of the most convenient sources to generate oncolytic viruses, however, the high seroprevalence of some enteroviruses limits its application which urges to exploit more oncolytic enteroviruses. In this study, coxsackievirus B5/Faulkner (CV-B5/F) was screened for its potential oncolytic effect against non-small cell lung cancers (NSCLCs) through inducing apoptosis and autophagy. For refractory NSCLCs, DNA-dependent protein kinase (DNA-PK) or ataxia telangiectasia mutated protein (ATM) inhibitors can synergize with CV-B5/F to promote refractory cell death. Here, we showed that viral infection triggered endoplasmic reticulum (ER) stress-related pro-apoptosis and autophagy signals, whereas repair for double-stranded DNA breaks (DSBs) contributed to cell survival which can be antagonized by inhibitor-induced cell death, manifesting exacerbated DSBs, apoptosis, and autophagy. Mechanistically, PERK pathway was activated by the combination of CV-B5/F and inhibitor, and the irreversible ER stress-induced exacerbated cell death. Furthermore, the degradation of activated STING by ERphagy promoted viral replication. Meanwhile, no treatment-related deaths due to CV-B5/F and/or inhibitors occurred. Conclusively, our study identifies an oncolytic CV-B5/F and the synergistic effects of inhibitors of DNA-PK or ATM, which is a potential therapy for NSCLCs.

The detection and surveillance of ship targets in coastal waters is not only a crucial technology for the advancement of ship intelligence, but also holds great significance for the safety and economic development of coastal areas. However, due to poor visibility in foggy conditions, the effectiveness of ship detection in coastal waters during foggy weather is limited. In this paper, we propose an improved version of YOLOv8s, termed YOLOv8s-Fog, which provides a multi-target detection network specifically designed for nearshore scenes in foggy weather. This improvement involves adding coordinate attention to the neck of YOLOv8 and replacing the convolution in C2f with deformable convolution. Additionally, to expand the dataset, we construct and synthesize a collection of ship target images captured in coastal waters on days with varying degrees of fog, using the atmospheric scattering model and monocular depth estimation. We compare the improved model with the standard YOLOv8s model, as well as several other object detection models. The results demonstrate superior performance achieved by the improved model, achieving an average accuracy of 74.4% (mAP@0.5), which is 1.2% higher than that achieved by the standard YOLOv8s model.

The authors introduce and demonstrate cross-comb spectroscopy in the mid-infrared as a variant of dual-comb spectroscopy. It provides enhanced performance and allows mid-infrared spectral information to be obtained by near-infrared detection.

Neutralizing antibodies (NtAbs) against severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) are indicators of vaccine efficacy that enable immunity surveillance. However, the rapid mutation of SARS-CoV-2 variants prevents the timely establishment of standards required for effective XBB vaccine evaluation. Therefore, we prepared four candidate standards (No. 11, No. 44, No. 22, and No. 33) using plasma, purified immunoglobulin, and a broad-spectrum neutralizing monoclonal antibody. Collaborative calibration was conducted across nine Chinese laboratories using neutralization methods against 11 strains containing the XBB and BA.2.86 sublineages. This study demonstrated the reduced neutralization potency of the first International Standard antibodies to SARS-CoV-2 variants of concern against XBB variants. No. 44 displayed broad-spectrum neutralizing activity against XBB sublineages, effectively reduced interlaboratory variability for nearly all XBB variants, and effectively minimized the geometric mean titer (GMT) difference between the live and pseudotyped virus. No. 22 showed a broader spectrum and higher neutralizing activity against all strains but failed to reduce interlaboratory variability. Thus, No. 44 was approved as a National Standard for NtAbs against XBB variants, providing a unified NtAb measurement standard for XBB variants for the first time. Moreover, No. 22 was approved as a national reference reagent for NtAbs against SARS-CoV-2, offering a broad-spectrum activity reference for current and potentially emerging variants.

Abstract To investigate the protective efficacy and mechanism of ZF2001 (a protein subunit vaccine with conditional approval in China) to SARS-CoV-2 Delta variant-induced severe pneumonia, the lethal challenge model of K18-hACE2 transgenic mice was used in this study. An inactivated-virus vaccine at the research and development stage (abbreviated as RDINA) was compared to ZF2001. We found that ZF2001 and RDINA could provide the protective effect against Delta variant-induced severe cases, as measured by the improved survival rates, the reduced virus loads, the alleviated lung histopathology and the high neutralizing antibody geomean titers, compared to aluminum adjuvant group. To prevent and control Omicron or other variant epidemics, further improvements in vaccine design and compatibilities with the novel adjuvant are required to achieve better immunogenicity.

Neutralizing antibody (NtAb) levels are key indicators in the development and evaluation of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) vaccines. Establishing a unified and reliable WHO International Standard (IS) for NtAb is crucial for the calibration and harmonization of NtAb detection assays. National and other WHO secondary standards are key links in the transfer of IS to working standards but are often overlooked. The Chinese National Standard (NS) and WHO IS were developed by China and WHO in September and December 2020, respectively, the application of which prompted and coordinated sero-detection of vaccine and therapy globally. Currently, a second-generation Chinese NS is urgently required owing to the depletion of stocks and need for calibration to the WHO IS. The Chinese National Institutes for Food and Drug Control (NIFDC) developed two candidate NSs (samples 33 and 66–99) traced to the IS according to the WHO manual for the establishment of national secondary standards through a collaborative study of nine experienced labs. Either NS candidate can reduce the systematic error among different laboratories and the difference between the live virus neutralization (Neut) and pseudovirus neutralization (PsN) methods, ensuring the accuracy and comparability of NtAb test results among multiple labs and methods, especially for samples 66–99. At present, samples 66–99 have been approved as the second-generation NS, which is the first NS calibrated tracing to the IS with 580 (460–740) International Units (IU)/mL and 580 (520–640) IU/mL by Neut and PsN, respectively. The use of standards improves the reliability and comparability of NtAb detection, ensuring the continuity of the use of the IS unitage, which effectively promotes the development and application of SARS-CoV-2 vaccines in China.

Abstract There are nearly 40% of cervical cancer patients showing poor response to neoadjuvant chemotherapy that can be induced by autophagy, however, the underlying mechanism has not yet been fully clarified. We previously found that Sex-determining region of Y-related high-mobility-group box 6 (SOX6), a tumor suppressor gene or oncogene in several cancers, could induce autophagy in cervical cancer. Accordingly, this study aims to investigate the mechanism of SOX6-induced autophagy and its potential significance in the platinum-based chemotherapy of cervical cancer. Firstly, we found that SOX6 could promote autophagy in cervical cancer cells depending on its HMG domain. Mitogen-activated protein kinase kinase kinase kinase-4 (MAP4K4) gene was identified as the direct target gene of SOX6, which was transcriptionally upregulated by binding the HMG domain of SOX6 protein to its double-binding sites within MAP4K4 gene promoter. MAP4K4 mediated the SOX6-induced autophagy through inhibiting PI3K-Akt-mTOR pathway and activating MAPK/ERK pathway. Further, the sensitivity of cervical cancer cells to cisplatin chemotherapy could be reduced by the SOX6-induced autophagy in vitro and in vivo, while such a phenomenon could be turned over by autophagy-specific inhibitor and MAP4K4 inhibitor, respectively. Moreover, cisplatin itself could promote the expression of endogenous SOX6 and subsequently the MAP4K4-mediated autophagy in cervical cancer cells, which might in turn reduce the sensitivity of these cells to cisplatin treatment. These findings uncovered the underlying mechanism and potential significance of SOX6-induced autophagy, and shed new light on the usage of MAP4K4 inhibitor or autophagy-specific inhibitor for sensitizing cervical cancer cells to the platinum-based chemotherapy.

Mumps, a disease caused by the mumps virus (MuV), has been spread widely across the world, especially among children and adolescents. Recent frequent local mumps outbreaks were reported worldwide, which may be caused by the decline in the neutralization ability of the existing attenuated live mumps vaccines against circulating MuV strains which were different from the genotype A or B vaccine strains. There is an urgent need to understand the genotypes of MuV strains currently circulated globally and in China. The gene sequences of MuV strains circulated globally were collected and phylogenetic trees were constructed using different strategies. The results showed that the MuV strains previously circulated globally were predominantly genotype G, while genotype F was predominantly circulated in China, followed by genotype G. The molecular evolution of genotype F MuV strains circulated in China is at a low genetic mutation rate, and the analysis of population dynamics pattern indicates that the incidence of genotype F mumps in China showed a rebound trend. These findings provide a basis for the selection or design of vaccine strains, and the decision of the evaluation strains for immunogenicity and protective efficacy, which laid the foundation for the research and development, as well as the application of next-generation MuV vaccines.

As a sedimentary mineral, most sandstone type uranium deposits are formed in petroliferous basins. Therefore, we can fully tap the residual economic value of historical logging and 3D seismic data measured for oil and gas to search for sandstone type uranium deposits. However, a large amount of acoustic logging data are missing in the target stratum of the uranium reservoir in that it is not the main stratum of oil and gas. A reconstructed method of acoustic logging data based on clustering analysis and with the low-frequency compensation of deterministic inversion is proposed to solve this problem. Secondly, we can use these logging data with seismic data to obtain the 3D inversion data volume representing the sand body of the uranium reservoir based on seismic lithological inversion. Then, we can also delimit the 3D spatial range of sandstone type uranium deposits in petroliferous basins based on the calibration of uranium anomaly and sub-body detection. Finally, a 3D field data example is given to test and analyze the effectiveness of the above research schemes.

The precision of phase estimation with interferometers can be greatly enhanced using non-classical quantum states, and the SU(1,1) interferometer is an elegant scheme, which generates two-mode squeezed state internally and also amplifies the signal. It has been shown in Anderson et al (2017 Phys. Rev. A 95 063843) that the photon-number measurement can achieve the Heisenberg limit, but only for estimating a small phase shift. We relax the constraint on the range of phase by considering two adaptive protocols: one also uses the photon-number measurement with a specially tuned sequence of feedback phase; the other implements the yet-to-be-realized optimal measurement but without fine tuning.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) protein subunit vaccine is one of the mainstream technology platforms for the development of COVID-19 vaccines, and most R&D units use the receptor-binding domain (RBD) or spike (S) protein as the main target antigen. The complexity of vaccine design, sequence, and expression systems makes it urgent to establish common antigen assays to facilitate vaccine development. In this study, we report the development of a double-antibody sandwich enzyme-linked immunosorbent assay (ELISA) to determine the antigen content of SARS-CoV-2 protein subunit vaccines based on the United States Pharmacopeia and ICH (international conference on harmonization) Q14 and Q2 (R2) requirements. A monoclonal antibody (mAb), 20D8, was identified as the detection antibody based on its high RBD binding activity (EC50 = 8.4 ng/mL), broad-spectrum anti-variant neutralizing activity (EC50: 2.7–9.8 ng/mL for pseudovirus and EC50: 9.6–127 ng/mL for authentic virus), good in vivo protection, and a recognized linear RBD epitope (369–379 aa). A porcine anti-RBD polyclonal antibody was selected as the coating antibody. Assay performance met the requirements of the analytical target profile with an accuracy and precision of ≥90% and adequate specificity. Within the specification range of 70–143%, the method capability index was >0.96; the misjudgment probability was

This study reconstructed molecular networks of human immunodeficiency virus (HIV) transmission history in an area affected by an epidemic of multiple HIV-1 subtypes and assessed the efficacy of strengthened early antiretroviral therapy (ART) and regular interventions in preventing HIV spread. We collected demographic and clinical data of 2221 treatment-naïve HIV-1–infected patients in a long-term cohort in Shenyang, Northeast China, between 2008 and 2016. HIV pol gene sequencing was performed and molecular networks of CRF01_AE, CRF07_BC, and subtype B were inferred using HIV-TRACE with separate optimized genetic distance threshold. We identified 168 clusters containing ≥ 2 cases among CRF01_AE-, CRF07_BC-, and subtype B-infected cases, including 13 large clusters (≥ 10 cases). Individuals in large clusters were characterized by younger age, homosexual behavior, more recent infection, higher CD4 counts, and delayed/no ART (P < 0.001). The dynamics of large clusters were estimated by proportional detection rate (PDR), cluster growth predictor, and effective reproductive number (Re). Most large clusters showed decreased or stable during the study period, indicating that expansion was slowing. The proportion of newly diagnosed cases in large clusters declined from 30 to 8% between 2008 and 2016, coinciding with an increase in early ART within 6 months after diagnosis from 24 to 79%, supporting the effectiveness of strengthened early ART and continuous regular interventions. In conclusion, molecular network analyses can thus be useful for evaluating the efficacy of interventions in epidemics with a complex HIV profile.

Fast Marching Method (FMM) boasts high calculation efficiency and strong adaptability and stability while being applied to seismic traveltime. However, when it is applied to the largescale model, the calculation precision of FMM is insufficient. FMM has poor calculation precision near the source, which is an essential reason for the low accuracy of the whole algorithm. This paper puts forward a joint traveltime calculation method to address the problem. Wavefront Construction (WFC) with a relatively high calculation accuracy rather than FMM is adopted for calculation of the grid nodes’ traveltime near the source. After that, FMM is used to calculate the seismic traveltime in the remaining area. Joint traveltime calculation method greatly improves the calculation accuracy of the source’s surrounding area and the calculation accuracy of FMM. According to the new method, FMM is still adopted for the calculation of most grid nodes in the model, so the high calculation efficiency of FMM is maintained. Multiple numerical models are utilized to verify the above conclusions in the paper.

Kirchhoff beam migration (KBM) is a ray-based seismic imaging method, which can handle multi-arrivals caused by model complexity. Apart from its high imaging precision, it also retains the merits of Kirchhoff migration, such as efficiency, stability, and flexibility. However, two issues should be taken into consideration when this method is expanded to the complicated surface conditions: first, the computational accuracy deficiency of the original local plane-wave decomposition method cannot suit for low signalto-noise ratio seismic data; second, as the rays traveling, the beam width increases rapidly, which cannot guarantee the computational accuracy of the corresponding grid points' attribute information. In addition, the insufficient coverage of the beam in the shallow part of the model might affect the imaging quality of this region. Kirchhoff dynamic focused beam migration based on compressed sensing is proposed to resolve these two problems. For the first problem, the local plane-wave decomposition method based on compressed sensing is introduced into KBM to enhance its computational accuracy. To solve the second problem, this paper adopts the dynamic focused beam to replace the original simplified Gaussian beam in the migration method, control the divergence of beam, and increase the coverage of beam in the shallow part of the model. Both Marmousi model and Canadian Foothills model are employed in this paper to test the new migration imaging method.

Kirchhoff beam migration (KBM) is a seismic-signal-based imaging method that considers both computational efficiency and computational precision. This method requires the decomposition of the seismic records, at each center-window, to plane waves in different directions. KBM achieves this step via the conventional linear Radon transform (LRT), which has obvious drawbacks. First, many noises and spatial aliasing exist in the transform result; and second, the resolution of the result does not meet a pre-defined threshold. These drawbacks could affect the imaging ability of KBM. To solve these problems, this paper introduced a high-resolution LRT method based on compressive sensing to KBM to improve the quality of the obtained plane waves and imaging results. This paper verifies the seismic-signal imaging method via multiple numerical models.

This paper aims at the problem that in deep water the migration quality is poor, and starts with the influence that velocity model accuracy has on migration, studying influence that variable seawater velocity makes on migration effect. At first, variable seawater velocity influenced by temperature, pressure and salinity is defined to replace the true seawater velocity. Then variable seawater velocity’s influence on interface migration location, layer sickness and migration energy focusing degree are analyzed in theory. And finally a deep water layered medium model containing variable seawater velocity, a syncline wedge shape model and a complex seafloor velocity model are constructed. By changing the seawater velocity of each model and comparing migration results of constant seawater-velocity model and variable seawater-velocity model, we can draw the conclusion: Under the condition of deep water, variable seawater-velocity’s impact on the quality of seismic migration is significant, which not only can change the location of geologic body migration result, but also can influence the resolution of geologic interface in the migration section and maybe can cause migration illusion. Investigación de los efectos de la variación en la velocidad del agua marina sobre las imágenes de migración en la geología de aguas profundas Resumen Este artículo se enfoca en el problema de la baja calidad de la migración en aguas profundas. Se analiza la influencia que tiene el modelo de precisión de velocidad en la migración y se estudia el impacto que la variación de velocidad del agua marina tiene en el efecto de movimiento. En primera instancia, se define la variación de la velocidad del agua marina afectada por la temperatura, la presión y la salinidad para reemplazar la velocidad del agua marina actual. Luego se analiza la teoría de la influencia de la velocidad del agua marina sobre la interfaz de la ubicación de migración, el grosor de la capa y la energía de movimiento de acuerdo con la inclinación. Además, finalmente, se construyó un modelo medio por capas de aguas profundas que contiene las variaciones de velocidad del agua marina, un modelo de forma de cuña sinclinal, y un modelo de complejo de velocidades del lecho marino. Al cambiar la velocidad del agua marina en cada modelo y comparar los resultados de migración del modelo de velocidad de agua marina constante con los del modelo modelo de velocidad variable se encontró la conclusión. Bajo las condiciones de agua profunda, el impacto de la variación de la velocidad del agua es significativo en la migración sísmica, lo que no solo puede cambiar el resultado de migración en la ubicación del cuerpo geológico sino que además puede influir en la resolución de la interconexión geológica en la sección de movimiento y podría causar la ilusión de migración.

In seismology, the local slope is essential for various seismic forward and inversion methods, including seismic migration imaging, seismic tomography, structure prediction, and other core seismic data processing methods. However, high-frequency random noise is unavoidable in real seismic data. Therefore, a robust estimation of event-slope attributes becomes important under high-frequency random noise. This research explores a robust and efficient method for estimating local slopes. The relationship between the frequency response of the Hilbert transform and Fourier transform is analyzed. The derivative operator constructed by the Hilbert transform can effectively weaken the effect on the energy enhancement of the high-frequency random noise. Minimizing the quadratic function of the proportionality factors was used to obtain the proportionality factors of the noise correction. Finally, the Hilbert transform and noise correction derived a new linear plane-wave destruction filter operator. The linear operator can effectively estimate the seismic events' local slope with high-frequency random noise. The calculation results of the numerical examples show that the linear plane-wave destruction filter estimation local slope method proposed in this study has high calculation accuracy and efficiency. [ABSTRACT FROM AUTHOR]