Suboptimal subspace construction for log-determinant approximation

Variance reduction is a crucial idea for Monte Carlo simulation and the stochastic Lanczos quadrature method is a dedicated method to approximate the trace of a matrix function. Inspired by their advantages, we combine these two techniques to approximate the log-determinant of large-scale symmetric positive definite matrices. Key questions to be answered for such a method are how to construct or choose an appropriate projection subspace and derive guaranteed theoretical analysis. This paper applies some probabilistic approaches including the projection-cost-preserving sketch and matrix concentration inequalities to construct a suboptimal subspace. Furthermore, we provide some insights on choosing design parameters in the underlying algorithm by deriving corresponding approximation error and probabilistic error estimations. Numerical experiments demonstrate our method's effectiveness and illustrate the quality of the derived error bounds

An analysis on stochastic Lanczos quadrature with asymmetric quadrature nodes

The stochastic Lanczos quadrature method has garnered significant attention recently. Upon examination of the error analyses given by Ubaru, Chen and Saad and Cortinovis and Kressner, certain notable inconsistencies arise. It turns out that the former's results are valid for cases with symmetric quadrature nodes and may not be adequate for many practical cases such as estimating log determinant of matrices. This paper analyzes probabilistic error bound of the stochastic Lanczos quadrature method for cases with asymmetric quadrature nodes. Besides, an optimized error allocation technique is employed to minimize the overall number of matrix vector multiplications required by the stochastic Lanczos quadrature method.Comment: 20 pages, 3 figure

CSP and Takeout Genes Modulate the Switch between Attraction and Repulsion during Behavioral Phase Change in the Migratory Locust

Behavioral plasticity is the most striking trait in locust phase transition. However, the genetic basis for behavioral plasticity in locusts is largely unknown. To unravel the molecular mechanisms underlying the behavioral phase change in the migratory locust Locusta migratoria, the gene expression patterns over the time courses of solitarization and gregarization were compared by oligonucleotide microarray analysis. Data analysis revealed that several gene categories relevant to peripheral olfactory perception are strongly regulated in a total of 1,444 differentially expressed genes during both time courses. Among these candidate genes, several CSP (chemosensory protein) genes and one takeout gene, LmigTO1, showed higher expression in gregarious and solitarious locusts, respectively, and displayed opposite expression trends during solitarization and gregarization. qRT-PCR experiments revealed that most CSP members and LmigTO1 exhibited antenna-rich expressions. RNA interference combined with olfactory behavioral experiments confirmed that the CSP gene family and one takeout gene, LmigTO1, are involved in the shift from repulsion to attraction between individuals during gregarization and in the reverse transition during solitarization. These findings suggest that the response to locust-emitted olfactory cues regulated by CSP and takeout genes is involved in the behavioral phase change in the migratory locust and provide a previously undescribed molecular mechanism linked to the formation of locust aggregations

Multi-scenario renewable energy absorption capacity assessment method based on the attention-enhanced time convolutional network

As the penetration rate of renewable energy in modern power grids continues to increase, the assessment of renewable energy absorption capacity plays an increasingly important role in the planning and operation of power and energy systems. However, traditional methods for assessing renewable energy absorption capacity rely on complex mathematical modeling, resulting in low assessment efficiency. Assessment in a single scenario determined by the source-load curve is difficult because it fails to reflect the random fluctuation characteristics of the source-load, resulting in inaccurate assessment results. To address and solve the above challenges, this paper proposes a multi-scenario renewable energy absorption capacity assessment method based on an attention-enhanced time convolutional network (ATCN). First, a source-load scene set is generated based on a generative adversarial network (GAN) to accurately characterize the uncertainty on both sides of the source and load. Then, the dependence of historical time series information in multiple scenarios is fully mined using the attention mechanism and temporal convolution network (TCN). Finally, simulation and experimental verification are carried out using a provincial power grid located in southwest China. The results show that the method proposed in this article has higher evaluation accuracy and speed than the traditional model

Search for the decay J/ψ→γ+invisibleJ/\psi\to\gamma + \rm {invisible}

We search for $J/\psi$ radiative decays into a weakly interacting neutral particle, namely an invisible particle, using the $J/\psi$ produced through the process $\psi(3686)\to\pi^+\pi^-J/\psi$ in a data sample of $(448.1\pm2.9)\times 10^6$ $\psi(3686)$ decays collected by the BESIII detector at BEPCII. No significant signal is observed. Using a modified frequentist method, upper limits on the branching fractions are set under different assumptions of invisible particle masses up to 1.2 $\mathrm{\ Ge\kern -0.1em V}/c^2$. The upper limit corresponding to an invisible particle with zero mass is 7.0$\times 10^{-7}$ at the 90\% confidence level

Observation of ηc→ωω\eta_c\to\omega\omega in J/ψ→γωωJ/\psi\to\gamma\omega\omega

Using a sample of $(1310.6\pm7.0)\times10^6$ $J/\psi$ events recorded with the BESIII detector at the symmetric electron positron collider BEPCII, we report the observation of the decay of the $(1^1 S_0)$ charmonium state $\eta_c$ into a pair of $\omega$ mesons in the process $J/\psi\to\gamma\omega\omega$. The branching fraction is measured for the first time to be $\mathcal{B}(\eta_c\to\omega\omega)= (2.88\pm0.10\pm0.46\pm0.68)\times10^{-3}$, where the first uncertainty is statistical, the second systematic and the third is from the uncertainty of $\mathcal{B}(J/\psi\to\gamma\eta_c)$. The mass and width of the $\eta_c$ are determined as $M=(2985.9\pm0.7\pm2.1)\,$MeV/$c^2$ and $\Gamma=(33.8\pm1.6\pm4.1)\,$MeV.Comment: 13 pages, 6 figure

Observation and study of the decay J/ψ→ϕηη′J/\psi\rightarrow\phi\eta\eta'

We report the observation and study of the decay $J/\psi\rightarrow\phi\eta\eta'$ using $1.3\times{10^9}$ $J/\psi$ events collected with the BESIII detector. Its branching fraction, including all possible intermediate states, is measured to be $(2.32\pm0.06\pm0.16)\times{10^{-4}}$. We also report evidence for a structure, denoted as $X$, in the $\phi\eta'$ mass spectrum in the $2.0-2.1$ GeV/$c^2$ region. Using two decay modes of the $\eta'$ meson ($\gamma\pi^+\pi^-$ and $\eta\pi^+\pi^-$), a simultaneous fit to the $\phi\eta'$ mass spectra is performed. Assuming the quantum numbers of the $X$ to be $J^P = 1^-$, its significance is found to be 4.4$\sigma$, with a mass and width of $(2002.1 \pm 27.5 \pm 21.4)$ MeV/$c^2$ and $(129 \pm 17 \pm 9)$ MeV, respectively, and a product branching fraction $\mathcal{B}(J/\psi\rightarrow\eta{}X)\times{}\mathcal{B}(X\rightarrow\phi\eta')=(9.8 \pm 1.2 \pm 1.7)\times10^{-5}$. Alternatively, assuming $J^P = 1^+$, the significance is 3.8$\sigma$, with a mass and width of $(2062.8 \pm 13.1 \pm 7.2)$ MeV/$c^2$ and $(177 \pm 36 \pm 35)$ MeV, respectively, and a product branching fraction $\mathcal{B}(J/\psi\rightarrow\eta{}X)\times{}\mathcal{B}(X\rightarrow\phi\eta')=(9.6 \pm 1.4 \pm 2.0)\times10^{-5}$. The angular distribution of $J/\psi\rightarrow\eta{}X$ is studied and the two $J^P$ assumptions of the $X$ cannot be clearly distinguished due to the limited statistics. In all measurements the first uncertainties are statistical and the second systematic.Comment: 10 pages, 6 figures and 4 table

Evidence of a resonant structure in the e+e−→π+D0D∗−e^+e^-\to \pi^+D^0D^{*-} cross section between 4.05 and 4.60 GeV

The cross section of the process $e^+e^-\to \pi^+D^0D^{*-}$ for center-of-mass energies from 4.05 to 4.60~GeV is measured precisely using data samples collected with the BESIII detector operating at the BEPCII storage ring. Two enhancements are clearly visible in the cross section around 4.23 and 4.40~GeV. Using several models to describe the dressed cross section yields stable parameters for the first enhancement, which has a mass of 4228.6 \pm 4.1 \pm 6.3 \un{MeV}/c^2 and a width of 77.0 \pm 6.8 \pm 6.3 \un{MeV}, where the first uncertainties are statistical and the second ones are systematic. Our resonant mass is consistent with previous observations of the $Y(4220)$ state and the theoretical prediction of a $D\bar{D}_1(2420)$ molecule. This result is the first observation of $Y(4220)$ associated with an open-charm final state. Fits with three resonance functions with additional $Y(4260)$, $Y(4320)$, $Y(4360)$, $\psi(4415)$, or a new resonance, do not show significant contributions from either of these resonances. The second enhancement is not from a single known resonance. It could contain contributions from $\psi(4415)$ and other resonances, and a detailed amplitude analysis is required to better understand this enhancement

Observation of Ds+→pnˉD^+_s\rightarrow p\bar{n} and confirmation of its large branching fraction

The baryonic decay $D^+_s\rightarrow p\bar{n}$ is observed, and the corresponding branching fraction is measured to be $(1.21\pm0.10\pm0.05)\times10^{-3}$, where the first uncertainty is statistical and second systematic. The data sample used in this analysis was collected with the BESIII detector operating at the BEPCII $e^+e^-$ double-ring collider with a center-of-mass energy of 4.178~GeV and an integrated luminosity of 3.19~fb$^{-1}$. The result confirms the previous measurement by the CLEO Collaboration and is of greatly improved precision, which may deepen our understanding of the dynamical enhancement of the W-annihilation topology in the charmed meson decays