In this paper, we calculate the absolute tensor square of the Dirichlet L-functions and show that it is expressed as an Euler product over pairs of primes. The method is to construct an equation to link primes to a se...In this paper, we calculate the absolute tensor square of the Dirichlet L-functions and show that it is expressed as an Euler product over pairs of primes. The method is to construct an equation to link primes to a series which has the factors of the absolute tensor product of the Dirichlet L-functions. This study is a generalization of Akatsuka’s theorem on the Riemann zeta function, and gives a proof of Kurokawa’s prediction proposed in 1992.展开更多
The subthalamic nucleus(STN)is considered the best target for deep brain stimulation treatments of Parkinson’s disease(PD).It is difficult to localize the STN due to its small size and deep location.Multichannel micr...The subthalamic nucleus(STN)is considered the best target for deep brain stimulation treatments of Parkinson’s disease(PD).It is difficult to localize the STN due to its small size and deep location.Multichannel microelectrode arrays(MEAs)can rapidly and precisely locate the STN,which is important for precise stimulation.In this paper,16-channel MEAs modified with multiwalled carbon nanotube/poly(3,4-ethylenedioxythiophene):poly(styrene sulfonate)(MWCNT/PEDOT:PSS)nanocomposites were designed and fabricated,and the accurate and rapid identification of the STN in PD rats was performed using detection sites distributed at different brain depths.These results showed that nuclei in 6-hydroxydopamine hydrobromide(6-OHDA)-lesioned brains discharged more intensely than those in unlesioned brains.In addition,the MEA simultaneously acquired neural signals from both the STN and the upper or lower boundary nuclei of the STN.Moreover,higher values of spike firing rate,spike amplitude,local field potential(LFP)power,and beta oscillations were detected in the STN of the 6-OHDA-lesioned brain,and may therefore be biomarkers of STN localization.Compared with the STNs of unlesioned brains,the power spectral density of spikes and LFPs synchronously decreased in the delta band and increased in the beta band of 6-OHDA-lesioned brains.This may be a cause of sleep and motor disorders associated with PD.Overall,this work describes a new cellular-level localization and detection method and provides a tool for future studies of deep brain nuclei.展开更多
For any given positive integer n ≥ 1, the Euler function φ(n) is defined to be the number of positive integers not exceeding n which are relatively prime to n. w(n) is defined to be the number of different prime...For any given positive integer n ≥ 1, the Euler function φ(n) is defined to be the number of positive integers not exceeding n which are relatively prime to n. w(n) is defined to be the number of different prime divisors of n. Some kind of equations involving Euler's function is studied in the paper.展开更多
Using Euler’s first-order explicit(EE)method and the peridynamic differential operator(PDDO)to discretize the time and internal crystal-size derivatives,respectively,the Euler’s first-order explicit method–peridyna...Using Euler’s first-order explicit(EE)method and the peridynamic differential operator(PDDO)to discretize the time and internal crystal-size derivatives,respectively,the Euler’s first-order explicit method–peridynamic differential operator(EE–PDDO)was obtained for solving the one-dimensional population balance equation in crystallization.Four different conditions during crystallization were studied:size-independent growth,sizedependent growth in a batch process,nucleation and size-independent growth,and nucleation and size-dependent growth in a continuous process.The high accuracy of the EE–PDDO method was confirmed by comparing it with the numerical results obtained using the second-order upwind and HR-van methods.The method is characterized by non-oscillation and high accuracy,especially in the discontinuous and sharp crystal size distribution.The stability of the EE–PDDO method,choice of weight function in the PDDO method,and optimal time step are also discussed.展开更多
In this paper,by choosing some appropriate test functions,we prove the Weyl’s lemma for triharmonic functions based on the new type of mean value formulas.
A generalized form of the error function, Gp(x)=pΓ(1/p)∫0xe−tpdt, which is directly associated with the gamma function, is evaluated for arbitrary real values of p>1and 0x≤+∞by employing a fast-converging power...A generalized form of the error function, Gp(x)=pΓ(1/p)∫0xe−tpdt, which is directly associated with the gamma function, is evaluated for arbitrary real values of p>1and 0x≤+∞by employing a fast-converging power series expansion developed in resolving the so-called Grandi’s paradox. Comparisons with accurate tabulated values for well-known cases such as the error function are presented using the expansions truncated at various orders.展开更多
This paper introduces a hybrid approach combining Green’s function Monte Carlo(GFMC)method with projected entangled pair state(PEPS)ansatz.This hybrid method regards PEPS as a trial state and a guiding wave function ...This paper introduces a hybrid approach combining Green’s function Monte Carlo(GFMC)method with projected entangled pair state(PEPS)ansatz.This hybrid method regards PEPS as a trial state and a guiding wave function in GFMC.By leveraging PEPS’s proficiency in capturing quantum state entanglement and GFMC’s efficient parallel architecture,the hybrid method is well-suited for the accurate and efficient treatment of frustrated quantum spin systems.As a benchmark,we applied this approach to study the frustrated J_(1)–J_(2) Heisenberg model on a square lattice with periodic boundary conditions(PBCs).Compared with other numerical methods,our approach integrating PEPS and GFMC shows competitive accuracy in the performance of ground-state energy.This paper provides systematic and comprehensive discussion of the approach of our previous work[Phys.Rev.B 109235133(2024)].展开更多
The solution of Poisson’s Equation plays an important role in many areas, including modeling high-intensity and high-brightness beams in particle accelerators. For the computational domain with a large aspect ratio, ...The solution of Poisson’s Equation plays an important role in many areas, including modeling high-intensity and high-brightness beams in particle accelerators. For the computational domain with a large aspect ratio, the integrated Green’s function method has been adopted to solve the 3D Poisson equation subject to open boundary conditions. In this paper, we report on the efficient implementation of this method, which can save more than a factor of 50 computing time compared with the direct brute force implementation and its improvement under certain extreme conditions.展开更多
Let φ(n) denote the Euler-totient function, we study the distribution of solutions of φ(n) ≤ x in arithmetic progressions, where n ≡ l(mod q) and an asymptotic formula was obtained by Perron formula.
Understanding the neural underpinning of human gait and balance is one of the most pertinent challenges for 21st-century translational neuroscience due to the profound impact that falls and mobility disturbances have ...Understanding the neural underpinning of human gait and balance is one of the most pertinent challenges for 21st-century translational neuroscience due to the profound impact that falls and mobility disturbances have on our aging population.Posture and gait control does not happen automatically,as previously believed,but rather requires continuous involvement of central nervous mechanisms.To effectively exert control over the body,the brain must integrate multiple streams of sensory information,including visual,vestibular,and somatosensory signals.The mechanisms which underpin the integration of these multisensory signals are the principal topic of the present work.Existing multisensory integration theories focus on how failure of cognitive processes thought to be involved in multisensory integration leads to falls in older adults.Insufficient emphasis,however,has been placed on specific contributions of individual sensory modalities to multisensory integration processes and cross-modal interactions that occur between the sensory modalities in relation to gait and balance.In the present work,we review the contributions of somatosensory,visual,and vestibular modalities,along with their multisensory intersections to gait and balance in older adults and patients with Parkinson’s disease.We also review evidence of vestibular contributions to multisensory temporal binding windows,previously shown to be highly pertinent to fall risk in older adults.Lastly,we relate multisensory vestibular mechanisms to potential neural substrates,both at the level of neurobiology(concerning positron emission tomography imaging)and at the level of electrophysiology(concerning electroencephalography).We hope that this integrative review,drawing influence across multiple subdisciplines of neuroscience,paves the way for novel research directions and therapeutic neuromodulatory approaches,to improve the lives of older adults and patients with neurodegenerative diseases.展开更多
This paper is concerned with the following fourth-order three-point boundary value problem , where , we discuss the existence of positive solutions to the above problem by applying to the fixed point theory in cones a...This paper is concerned with the following fourth-order three-point boundary value problem , where , we discuss the existence of positive solutions to the above problem by applying to the fixed point theory in cones and iterative technique.展开更多
For a set S of real numbers, we introduce the concept of S-almost automorphic functions valued in a Banach space. It generalizes in particular the space of Z-almost automorphic functions. Considering the space of S-al...For a set S of real numbers, we introduce the concept of S-almost automorphic functions valued in a Banach space. It generalizes in particular the space of Z-almost automorphic functions. Considering the space of S-almost automorphic functions, we give sufficient conditions of the existence and uniqueness of almost automorphic solutions of a differential equation with a piecewise constant argument of generalized type. This is done using the Banach fixed point theorem.展开更多
Long-term levodopa administration can lead to the development of levodopa-induced dyskinesia.Gamma oscillations are a widely recognized hallmark of abnormal neural electrical activity in levodopa-induced dyskinesia.Cu...Long-term levodopa administration can lead to the development of levodopa-induced dyskinesia.Gamma oscillations are a widely recognized hallmark of abnormal neural electrical activity in levodopa-induced dyskinesia.Currently,studies have reported increased oscillation power in cases of levodopa-induced dyskinesia.However,little is known about how the other electrophysiological parameters of gamma oscillations are altered in levodopa-induced dyskinesia.Furthermore,the role of the dopamine D3 receptor,which is implicated in levodopa-induced dyskinesia,in movement disorder-related changes in neural oscillations is unclear.We found that the cortico-striatal functional connectivity of beta oscillations was enhanced in a model of Parkinson’s disease.Furthermore,levodopa application enhanced cortical gamma oscillations in cortico-striatal projections and cortical gamma aperiodic components,as well as bidirectional primary motor cortex(M1)↔dorsolateral striatum gamma flow.Administration of PD128907(a selective dopamine D3 receptor agonist)induced dyskinesia and excessive gamma oscillations with a bidirectional M1↔dorsolateral striatum flow.However,administration of PG01037(a selective dopamine D3 receptor antagonist)attenuated dyskinesia,suppressed gamma oscillations and cortical gamma aperiodic components,and decreased gamma causality in the M1→dorsolateral striatum direction.These findings suggest that the dopamine D3 receptor plays a role in dyskinesia-related oscillatory activity,and that it has potential as a therapeutic target for levodopa-induced dyskinesia.展开更多
Nowadays,presynaptic dopaminergic positron emission tomography,which assesses deficiencies in dopamine synthesis,storage,and transport,is widely utilized for early diagnosis and differential diagnosis of parkinsonism....Nowadays,presynaptic dopaminergic positron emission tomography,which assesses deficiencies in dopamine synthesis,storage,and transport,is widely utilized for early diagnosis and differential diagnosis of parkinsonism.This review provides a comprehensive summary of the latest developments in the application of presynaptic dopaminergic positron emission tomography imaging in disorders that manifest parkinsonism.We conducted a thorough literature search using reputable databases such as PubMed and Web of Science.Selection criteria involved identifying peer-reviewed articles published within the last 5 years,with emphasis on their relevance to clinical applications.The findings from these studies highlight that presynaptic dopaminergic positron emission tomography has demonstrated potential not only in diagnosing and differentiating various Parkinsonian conditions but also in assessing disease severity and predicting prognosis.Moreover,when employed in conjunction with other imaging modalities and advanced analytical methods,presynaptic dopaminergic positron emission tomography has been validated as a reliable in vivo biomarker.This validation extends to screening and exploring potential neuropathological mechanisms associated with dopaminergic depletion.In summary,the insights gained from interpreting these studies are crucial for enhancing the effectiveness of preclinical investigations and clinical trials,ultimately advancing toward the goals of neuroregeneration in parkinsonian disorders.展开更多
The Owen’s T function is presented in four new ways, one of them as a series similar to the Euler’s arctangent series divided by 2π, which is its majorant series. All possibilities enable numerically stable ...The Owen’s T function is presented in four new ways, one of them as a series similar to the Euler’s arctangent series divided by 2π, which is its majorant series. All possibilities enable numerically stable and fast convergent computation of the bivariate normal integral with simple recursion. When tested computation on a random sample of one million parameter triplets with uniformly distributed components and using double precision arithmetic, the maximum absolute error was 3.45 × 10<sup>-</sup><sup>16</sup>. In additional testing, focusing on cases with correlation coefficients close to one in absolute value, when the computation may be very sensitive to small rounding errors, the accuracy was retained. In rare potentially critical cases, a simple adjustment to the computation procedure was performed—one potentially critical computation was replaced with two equivalent non-critical ones. All new series are suitable for vector and high-precision computation, assuming they are supplemented with appropriate efficient and accurate computation of the arctangent and standard normal cumulative distribution functions. They are implemented by the R package Phi2rho, available on CRAN. Its functions allow vector arguments and are ready to work with the Rmpfr package, which enables the use of arbitrary precision instead of double precision numbers. A special test with up to 1024-bit precision computation is also presented.展开更多
We present a method to unify the calculation of Green's functions for an electromagnetic(EM) transmitting source embedded in a homogeneous stratified medium.A virtual interface parallel to layer interfaces is intro...We present a method to unify the calculation of Green's functions for an electromagnetic(EM) transmitting source embedded in a homogeneous stratified medium.A virtual interface parallel to layer interfaces is introduced through the source location.The potentials for Green's function are derived by decomposing the partial wave solutions to Helmholtz's equations into upward and downward within boundaries.The amplitudes of the potentials in each stratum are obtained recursively from the initial amplitudes at the source level.The initial amplitudes are derived by coupling with the transmitting sources and following the discontinuity of the tangential electric and magnetic fields at the source interface.Only the initial terms are related to the transmitting sources and thus need to be modified for different transmitters,whereas the kernel connected with the stratified media stays unchanged.Hence,the present method can be easily applied to EM transmitting sources with little modification.The application of the proposed method to the marine controlled-source electromagnetic method(MCSEM) demonstrates its simplicity and flexibility.展开更多
文摘In this paper, we calculate the absolute tensor square of the Dirichlet L-functions and show that it is expressed as an Euler product over pairs of primes. The method is to construct an equation to link primes to a series which has the factors of the absolute tensor product of the Dirichlet L-functions. This study is a generalization of Akatsuka’s theorem on the Riemann zeta function, and gives a proof of Kurokawa’s prediction proposed in 1992.
基金funded by the National Natural Science Foundation of China(Nos.L2224042,T2293731,62121003,61960206012,61973292,62171434,61975206,and 61971400)the Frontier Interdisciplinary Project of the Chinese Academy of Sciences(No.XK2022XXC003)+2 种基金the National Key Research and Development Program of China(Nos.2022YFC2402501 and 2022YFB3205602)the Major Program of Scientific and Technical Innovation 2030(No.2021ZD02016030)the Scientific Instrument Developing Project of he Chinese Academy of Sciences(No.GJJSTD20210004).
文摘The subthalamic nucleus(STN)is considered the best target for deep brain stimulation treatments of Parkinson’s disease(PD).It is difficult to localize the STN due to its small size and deep location.Multichannel microelectrode arrays(MEAs)can rapidly and precisely locate the STN,which is important for precise stimulation.In this paper,16-channel MEAs modified with multiwalled carbon nanotube/poly(3,4-ethylenedioxythiophene):poly(styrene sulfonate)(MWCNT/PEDOT:PSS)nanocomposites were designed and fabricated,and the accurate and rapid identification of the STN in PD rats was performed using detection sites distributed at different brain depths.These results showed that nuclei in 6-hydroxydopamine hydrobromide(6-OHDA)-lesioned brains discharged more intensely than those in unlesioned brains.In addition,the MEA simultaneously acquired neural signals from both the STN and the upper or lower boundary nuclei of the STN.Moreover,higher values of spike firing rate,spike amplitude,local field potential(LFP)power,and beta oscillations were detected in the STN of the 6-OHDA-lesioned brain,and may therefore be biomarkers of STN localization.Compared with the STNs of unlesioned brains,the power spectral density of spikes and LFPs synchronously decreased in the delta band and increased in the beta band of 6-OHDA-lesioned brains.This may be a cause of sleep and motor disorders associated with PD.Overall,this work describes a new cellular-level localization and detection method and provides a tool for future studies of deep brain nuclei.
基金Foundation item: Supported by the National Natural Science Foundation of China(10671056)
文摘For any given positive integer n ≥ 1, the Euler function φ(n) is defined to be the number of positive integers not exceeding n which are relatively prime to n. w(n) is defined to be the number of different prime divisors of n. Some kind of equations involving Euler's function is studied in the paper.
文摘Using Euler’s first-order explicit(EE)method and the peridynamic differential operator(PDDO)to discretize the time and internal crystal-size derivatives,respectively,the Euler’s first-order explicit method–peridynamic differential operator(EE–PDDO)was obtained for solving the one-dimensional population balance equation in crystallization.Four different conditions during crystallization were studied:size-independent growth,sizedependent growth in a batch process,nucleation and size-independent growth,and nucleation and size-dependent growth in a continuous process.The high accuracy of the EE–PDDO method was confirmed by comparing it with the numerical results obtained using the second-order upwind and HR-van methods.The method is characterized by non-oscillation and high accuracy,especially in the discontinuous and sharp crystal size distribution.The stability of the EE–PDDO method,choice of weight function in the PDDO method,and optimal time step are also discussed.
基金Supported by National Natural Science Foundation of China(Grant Nos.11801006 and 12071489).
文摘In this paper,by choosing some appropriate test functions,we prove the Weyl’s lemma for triharmonic functions based on the new type of mean value formulas.
文摘A generalized form of the error function, Gp(x)=pΓ(1/p)∫0xe−tpdt, which is directly associated with the gamma function, is evaluated for arbitrary real values of p>1and 0x≤+∞by employing a fast-converging power series expansion developed in resolving the so-called Grandi’s paradox. Comparisons with accurate tabulated values for well-known cases such as the error function are presented using the expansions truncated at various orders.
基金Project supported by the National Natural Science Foundation of China(Grant No.11934020)the Innovation Program for Quantum Science and Technology(Grant No.2021ZD0302402).
文摘This paper introduces a hybrid approach combining Green’s function Monte Carlo(GFMC)method with projected entangled pair state(PEPS)ansatz.This hybrid method regards PEPS as a trial state and a guiding wave function in GFMC.By leveraging PEPS’s proficiency in capturing quantum state entanglement and GFMC’s efficient parallel architecture,the hybrid method is well-suited for the accurate and efficient treatment of frustrated quantum spin systems.As a benchmark,we applied this approach to study the frustrated J_(1)–J_(2) Heisenberg model on a square lattice with periodic boundary conditions(PBCs).Compared with other numerical methods,our approach integrating PEPS and GFMC shows competitive accuracy in the performance of ground-state energy.This paper provides systematic and comprehensive discussion of the approach of our previous work[Phys.Rev.B 109235133(2024)].
文摘The solution of Poisson’s Equation plays an important role in many areas, including modeling high-intensity and high-brightness beams in particle accelerators. For the computational domain with a large aspect ratio, the integrated Green’s function method has been adopted to solve the 3D Poisson equation subject to open boundary conditions. In this paper, we report on the efficient implementation of this method, which can save more than a factor of 50 computing time compared with the direct brute force implementation and its improvement under certain extreme conditions.
基金Supported by the National Natural Science Foundation of China(11271249) Supported by the Scientific and Technological Research Program of Chongqing Municipal Education Commission(1601213) Supported by the Scientific Research Program of Yangtze Normal University(2012XJYBO31)
文摘Let φ(n) denote the Euler-totient function, we study the distribution of solutions of φ(n) ≤ x in arithmetic progressions, where n ≡ l(mod q) and an asymptotic formula was obtained by Perron formula.
文摘Understanding the neural underpinning of human gait and balance is one of the most pertinent challenges for 21st-century translational neuroscience due to the profound impact that falls and mobility disturbances have on our aging population.Posture and gait control does not happen automatically,as previously believed,but rather requires continuous involvement of central nervous mechanisms.To effectively exert control over the body,the brain must integrate multiple streams of sensory information,including visual,vestibular,and somatosensory signals.The mechanisms which underpin the integration of these multisensory signals are the principal topic of the present work.Existing multisensory integration theories focus on how failure of cognitive processes thought to be involved in multisensory integration leads to falls in older adults.Insufficient emphasis,however,has been placed on specific contributions of individual sensory modalities to multisensory integration processes and cross-modal interactions that occur between the sensory modalities in relation to gait and balance.In the present work,we review the contributions of somatosensory,visual,and vestibular modalities,along with their multisensory intersections to gait and balance in older adults and patients with Parkinson’s disease.We also review evidence of vestibular contributions to multisensory temporal binding windows,previously shown to be highly pertinent to fall risk in older adults.Lastly,we relate multisensory vestibular mechanisms to potential neural substrates,both at the level of neurobiology(concerning positron emission tomography imaging)and at the level of electrophysiology(concerning electroencephalography).We hope that this integrative review,drawing influence across multiple subdisciplines of neuroscience,paves the way for novel research directions and therapeutic neuromodulatory approaches,to improve the lives of older adults and patients with neurodegenerative diseases.
文摘This paper is concerned with the following fourth-order three-point boundary value problem , where , we discuss the existence of positive solutions to the above problem by applying to the fixed point theory in cones and iterative technique.
文摘For a set S of real numbers, we introduce the concept of S-almost automorphic functions valued in a Banach space. It generalizes in particular the space of Z-almost automorphic functions. Considering the space of S-almost automorphic functions, we give sufficient conditions of the existence and uniqueness of almost automorphic solutions of a differential equation with a piecewise constant argument of generalized type. This is done using the Banach fixed point theorem.
基金supported by the National Natural Science Foundation of China,No.82071254(to WZ).
文摘Long-term levodopa administration can lead to the development of levodopa-induced dyskinesia.Gamma oscillations are a widely recognized hallmark of abnormal neural electrical activity in levodopa-induced dyskinesia.Currently,studies have reported increased oscillation power in cases of levodopa-induced dyskinesia.However,little is known about how the other electrophysiological parameters of gamma oscillations are altered in levodopa-induced dyskinesia.Furthermore,the role of the dopamine D3 receptor,which is implicated in levodopa-induced dyskinesia,in movement disorder-related changes in neural oscillations is unclear.We found that the cortico-striatal functional connectivity of beta oscillations was enhanced in a model of Parkinson’s disease.Furthermore,levodopa application enhanced cortical gamma oscillations in cortico-striatal projections and cortical gamma aperiodic components,as well as bidirectional primary motor cortex(M1)↔dorsolateral striatum gamma flow.Administration of PD128907(a selective dopamine D3 receptor agonist)induced dyskinesia and excessive gamma oscillations with a bidirectional M1↔dorsolateral striatum flow.However,administration of PG01037(a selective dopamine D3 receptor antagonist)attenuated dyskinesia,suppressed gamma oscillations and cortical gamma aperiodic components,and decreased gamma causality in the M1→dorsolateral striatum direction.These findings suggest that the dopamine D3 receptor plays a role in dyskinesia-related oscillatory activity,and that it has potential as a therapeutic target for levodopa-induced dyskinesia.
基金supported by the Research Project of the Shanghai Health Commission,No.2020YJZX0111(to CZ)the National Natural Science Foundation of China,Nos.82021002(to CZ),82272039(to CZ),82171252(to FL)+1 种基金a grant from the National Health Commission of People’s Republic of China(PRC),No.Pro20211231084249000238(to JW)Medical Innovation Research Project of Shanghai Science and Technology Commission,No.21Y11903300(to JG).
文摘Nowadays,presynaptic dopaminergic positron emission tomography,which assesses deficiencies in dopamine synthesis,storage,and transport,is widely utilized for early diagnosis and differential diagnosis of parkinsonism.This review provides a comprehensive summary of the latest developments in the application of presynaptic dopaminergic positron emission tomography imaging in disorders that manifest parkinsonism.We conducted a thorough literature search using reputable databases such as PubMed and Web of Science.Selection criteria involved identifying peer-reviewed articles published within the last 5 years,with emphasis on their relevance to clinical applications.The findings from these studies highlight that presynaptic dopaminergic positron emission tomography has demonstrated potential not only in diagnosing and differentiating various Parkinsonian conditions but also in assessing disease severity and predicting prognosis.Moreover,when employed in conjunction with other imaging modalities and advanced analytical methods,presynaptic dopaminergic positron emission tomography has been validated as a reliable in vivo biomarker.This validation extends to screening and exploring potential neuropathological mechanisms associated with dopaminergic depletion.In summary,the insights gained from interpreting these studies are crucial for enhancing the effectiveness of preclinical investigations and clinical trials,ultimately advancing toward the goals of neuroregeneration in parkinsonian disorders.
文摘The Owen’s T function is presented in four new ways, one of them as a series similar to the Euler’s arctangent series divided by 2π, which is its majorant series. All possibilities enable numerically stable and fast convergent computation of the bivariate normal integral with simple recursion. When tested computation on a random sample of one million parameter triplets with uniformly distributed components and using double precision arithmetic, the maximum absolute error was 3.45 × 10<sup>-</sup><sup>16</sup>. In additional testing, focusing on cases with correlation coefficients close to one in absolute value, when the computation may be very sensitive to small rounding errors, the accuracy was retained. In rare potentially critical cases, a simple adjustment to the computation procedure was performed—one potentially critical computation was replaced with two equivalent non-critical ones. All new series are suitable for vector and high-precision computation, assuming they are supplemented with appropriate efficient and accurate computation of the arctangent and standard normal cumulative distribution functions. They are implemented by the R package Phi2rho, available on CRAN. Its functions allow vector arguments and are ready to work with the Rmpfr package, which enables the use of arbitrary precision instead of double precision numbers. A special test with up to 1024-bit precision computation is also presented.
基金supported by CNSF(Granted No.40874050)Chinese High Technology Project(Granted No.2011YQ05006010)
文摘We present a method to unify the calculation of Green's functions for an electromagnetic(EM) transmitting source embedded in a homogeneous stratified medium.A virtual interface parallel to layer interfaces is introduced through the source location.The potentials for Green's function are derived by decomposing the partial wave solutions to Helmholtz's equations into upward and downward within boundaries.The amplitudes of the potentials in each stratum are obtained recursively from the initial amplitudes at the source level.The initial amplitudes are derived by coupling with the transmitting sources and following the discontinuity of the tangential electric and magnetic fields at the source interface.Only the initial terms are related to the transmitting sources and thus need to be modified for different transmitters,whereas the kernel connected with the stratified media stays unchanged.Hence,the present method can be easily applied to EM transmitting sources with little modification.The application of the proposed method to the marine controlled-source electromagnetic method(MCSEM) demonstrates its simplicity and flexibility.