Search

This paper is concerned with the low-rank approximation for large-scale nonsymmetric matrices. Inspired by the classical Nystrom method, which is a popular method to find the low-rank approximation for symmetric positive semidefinite matrices, we explore an extension of the Nystrom method to approximate nonsymmetric matrices. The proposed method is a generalized Nystrom method with column sketching and shows its advantages in accuracy and speed without sacri cing stability. And the numerical experiments will illustrate the robustness of our new methods in finding a desired low-rank approximation of nonsymmetric matrix.

To accurately predict the synthesis cross-sections of superheavy elements, identifying the optimal projectile-target combinations and the evaporation channels at specific collision energies, we have attempted to utilize high-quality microscopic nuclear masses (HQMNM) within the dinuclear system (DNS) model, which are obtained by fitting experimental data with the Skyrme energy density functional theory (DFT), as published in Phys. Lett. B 851 (2024) 138578. The atomic nuclear mass serves as a crucial input for the DNS model, as the Q-values and separation energies it generates directly influence the calculated fusion and survival probabilities. Our calculations have reproduced the experimental data for hot fusion and have been compared with results based on the finite-range droplet model (FRDM12) mass calculations. Compared to the FRDM12 mass results, we have found that the HQMNM provides a better fit to the experimental outcomes. For the specific reaction of \(^{48}\rm{Ca} + ^{243}\rm{Am} \rightarrow ^{291}\rm{Mc}^*\), we have conducted a detailed calculation of capture, fusion, and survival based on the HQMNM model and compared these with calculations based on other mass models. Based on these findings, we have systematically calculated available projectile target combinations for the synthesis of elements 119 and 120, and identified the optimal combinations. We provided the synthesis cross-sections, collision energies, and evaporation channels, offering a reference for conducting experiments on the synthesis of superheavy elements., Comment: 8 pages, 5 figures

To assess the potential of cold-fusion for synthesizing superheavy nuclei (SHN) with proton numbers 104-113, we systematically calculated 145 naturally occurring projectile-target combinations within the DNS model. Reactions predominantly show maximum cross-sections in the 1n to 2n channels, peaking near the Coulomb barrier with a sum of barrier and Q-value within 30 MeV. The maximum cross-section occurs below the Bass barrier, suggesting either the Bass model's limitation or significant deformation reducing the effective Coulomb barrier. Our calculations align well with experimental data, revealing that more neutron-rich projectiles slightly enhance fusion, though the effect is minor. For fixed targets (Pb, Bi), evaporation residue cross-sections decrease linearly with increasing projectile proton number, attributed to reduced fusion probability and lower fission barriers in heavier SHN. The touching potential $V_{\rm in}$ shows a linear trend with the product of projectile-target proton numbers, with neutron-rich systems exhibiting lower $V_{\rm in}$. Some reactions with $V_{\rm in} < V_{\rm S}$ may involve nucleon transfer before capture. Based on the DNS model, we identified optimal combinations and collision energies for synthesizing SHN with significant cross-sections. Collectively, our findings indicate that cold fusion is a promising avenue for creating proton-rich SHN around the drip line in the Z=104-113 region, offering distinct advantages over alternative mechanisms., Comment: 10 pages, 6 figures

Although quantum entanglement is an important resource, its characterization is quite challenging. The partial transposition is a common method to detect bipartite entanglement. In this paper, the authors study the partial-transpose(PT)-moments of two-qubit states,and completely describe the whole region, composed of the second and third PT-moments, for all two-qubit states. Furthermore, they determine the accurate region corresponding to all entangled two-qubit states. The states corresponding to those boundary points of the whole region, and to the border lines between separable and entangled states are analyzed. As an application, they characterize the entangled region of PT-moments for the two families of Werner states and Bell-diagonal states. The relations between entanglement and the pairs of PT-moments are revealed from these typical examples. They also numerically plot the whole region of possible PT-moments for all two-qubit X-states, and find that this region is almost the same as the whole region of PT-moments for all two-qubit states. Moreover, they extend their results to detect the entanglement of multiqubit states. By utilizing the PT-moment-based method to characterize the entanglement of the multiqubit states mixed by the GHZ and W states, they propose an operational way of verifying the genuine entanglement in such states., Comment: 31 pages, LaTeX, 9 figures

In this paper, we extend Chandrasekhar's method of calculating rotating black holes into $f(R)$ theory. We consider the Ricci scalar is a constant and derive the Kerr and Kerr-Ads metric by using the analytical mathematical method. Suppose that the spacetime is a 4-dimensional Riemannian manifold with a general stationary axisymmetric metric, we calculate Cartan's equation of structure and derive the Einstein tensor. In order to reduce the solving difficulty, we fix the gauge freedom to transform the metric into a more symmetric form. We solve the field equations in the two cases of the Ricci scalar $R=0$ and $R\neq 0$. In the case of $R=0$, the Ernst's equations are derived. We give the elementary solution of Ernst's equations and show the way to obtain more solutions including Kerr metric. In the case of $R\neq 0$, we reasonably assume that the solution to the equations consists of two parts: the first is Kerr part and the second is introduced by the Ricci scalar. Giving solution to the second part and combining the two parts, we obtain the Kerr-Ads metric. The calculations are carried out in a general $f(R)$ theory, indicating the Kerr and Kerr-Ads black holes exist widely in general $f(R)$ models. Furthermore, the whole solving process can be treated as a standard calculation procedure to obtain rotating black holes, which can be applied to other modified gravities.

Within the framework of the dinuclear system model, this study delves into the impact of various nuclear mass models on evaluating the fusion probability of superheavy nuclei. Nuclear mass models, as crucial inputs to the DNS model, exhibit slight variations in binding energy, quadrupole deformation, and extrapolation ability; these subtle differences can significantly influence the model's outcomes. Specifically, the study finds that nuclear mass plays a pivotal role in determining fusion probability, and Q-value. By numerically solving a set of master equations, the study examines how binding energies from different mass models affect the fusion probability of colliding nuclei, taking the example of $^{48}$Ca + $^{243}$Am $\rightarrow$ $^{291}$Mc. A careful analysis of the potential energy surface (PES) reveals that the inner fusion barriers lead to variations in fusion probabilities. Importantly, the study demonstrates that the synthesis cross sections of superheavy nuclei calculated using different nuclear mass models align well with experimental data, falling within an error range of one order of magnitude. This finding underscores the reliability of our model predictions. Looking ahead, the study utilizes five distinct nuclear mass models to predict the synthesis cross sections of superheavy elements 119 and 120, along with their associated uncertainties. These predictions offer valuable insights into the feasibility of synthesizing these elusive elements and pave the way for future experimental explorations.

We investigate the cosmological implications of a phantom dark energy model with bulk viscosity. We explore this model as a possible way to resolve the big rip singularity problem that plagues the phantom models. We use the latest type Ia supernova and Hubble parameter data to constrain the model parameters and find that the data favor a significant bulk viscosity over a non-constant potential term for the phantom field. We perform a dynamical analysis of the model and show that the only stable and physical attractor corresponds to a phantom-dominated era with a total equation of state that can be greater than $-1$ due to the viscosity. We also study the general effect of viscosity on the phantom field and the late time evolution of the universe. We apply the statefinder diagnostic to the model and find that it approaches a nearby fixed point asymptotically, indicating that the universe can escape the big rip singularity with the presence of bulk viscosity. We conclude that bulk viscosity can play an important role in affecting the late-time behavior as well as alleviating the singularity problem of the phantom universe., Comment: accepted by EPJC

We study the accretion of relativistic Vlasov gas onto a Kerr black hole, regarding the particles as distributed throughout all the space, other than just in the equatorial plane. We solve the relativistic Liouville equation in the full $3+1$ dimensional framework of Kerr geometry. For the flow that is stationary and axial symmetric, we prove that the distribution function is independent of the conjugate coordinates. For an explicit distribution that can approximate to Maxwell-J\"{u}ttner distribution, we further calculate the particle current density, the stress energy momentum tensor and the unit accretion rates of mass, energy and angular momentum. The analytic results at large distance are shown to be consistent with the limits of the numerical ones computed at finite distance. Especially, we show that the unit mass accretion rate agrees with the Schwarzschild result in the case of low temperature limit. Furthermore, we find from the numerical results that the three unit accretion rates vary with the angle in Kerr metric and the accretion of Vlasov gas would slow down the Kerr black hole. The closer to the equator, the faster it slows down the black hole., Comment: 28 pages, 11 figures, accepted for publication in Physical Review D

The bumblebee model is an extension of the Einstein-Maxwell theory that allows for the spontaneous breaking of the Lorentz symmetry of the spacetime. In this paper, we study the quasinormal modes of the spherical black holes in this model that are characterized by a global monopole. We analyze the two cases with a vanishing cosmological constant or a negative one (the anti-de Sitter case). We find that the black holes are stable under the perturbation of a massless scalar field. However, both the Lorentz symmetry breaking and the global monopole have notable impacts on the evolution of the perturbation. The Lorentz symmetry breaking may prolong or shorten the decay of the perturbation according to the sign of the breaking parameter. The global monopole, on the other hand, has different effects depending on whether a nonzero cosmological constant presences: it reduces the damping of the perturbations for the case with a vanishing cosmological constant, but has little influence for the anti-de Sitter case., Comment: accepted by Eur. Phys. J. C

Efficiently detecting entanglement based on measurable quantities is a basic problem for quantum information processing. Recently, the measurable quantities called partial-transpose (PT)-moments have been proposed to detect and characterize entanglement. In the recently published paper [L. Zhang \emph{et al.}, \href{https://doi.org/10.1002/andp.202200289}{Ann. Phys.(Berlin) \textbf{534}, 2200289 (2022)}], we have already identified the 2-dimensional (2D) region, comprised of the second and third PT-moments, corresponding to two-qubit entangled states, and described the whole region for all two-qubit states. In the present paper, we visualize the 3D region corresponding to all two-qubit states by further involving the fourth PT-moment (the last one for two-qubit states). The characterization of this 3D region can finally be achieved by optimizing some polynomials. Furthermore, we identify the dividing surface which separates the two parts of the whole 3D region corresponding to entangled and separable states respectively. Due to the measurability of PT-moments, we obtain a complete and operational criterion for the detection of two-qubit entanglement., Comment: 29 pages, LaTeX, 8 figures, 2 tables

We present a quantum algorithm for finding the minimum of a function based on multistep quantum computation and apply it for optimization problems with continuous variables, in which the variables of the problem are discretized to form the state space of the problem. Usually the cost for solving the problem increases dramatically with the size of the problem. In this algorithm, the dimension of the search space of the problem can be reduced exponentially step by step. We construct a sequence of Hamiltonians such that the search space of a Hamiltonian is nested in that of the previous one. By applying a multistep quantum computation process, the optimal vector is finally located in a small state space and can be determined efficiently. One of the most difficult problems in optimization is that a trial vector is trapped in a deep local minimum while the global minimum is missed, this problem can be alleviated in our algorithm and the runtime is proportional to the number of the steps of the algorithm, provided certain conditions are satisfied. We have tested the algorithm for some continuous test functions., Comment: 24 pages, 8 figures

The analysis of light variation of M87 can help us understand the disc evolution. In the past decade, M87 has experienced several short-term light variabilities related to flares. We also find there are year-scale X-ray variations in the core of M87. Their light variability properties are similar to clumpy-ADAF. By re-analyzing 56 $\it Chandra$ observations from 2007 to 2019, we distinguish the `non-flaring state' from `flaring state' in the light variability. After removing flaring state data, we identify 4 gas clumps in the nucleus and all of them can be well fitted by the clumpy-ADAF model. The average mass accretion rate is $\sim 0.16 \rm M_{\odot} yr^{-1}$. We analyze the photon index($\Gamma$)-flux(2-10keV) correlation between the non-flaring state and flaring state. For the non-flaring states, the flux is inversely proportional to the photon index. For the flaring states, we find no obvious correlation between the two parameters. In addition, we find that the flare always occurs at a high mass accretion rate, and after the luminosity of the flare reaches the peak, it will be accompanied by a sudden decrease in luminosity. Our results can be explained as that the energy released by magnetic reconnection destroys the structure of the accretion disc, thus the luminosity decreases rapidly and returns to normal levels thereafter.

The radial basis function (RBF) method is used for the numerical solution of the Poisson problem in high dimension. The approximate solution can be found by solving a large system of linear equations. Here we investigate the extent to which the RBF method can be accelerated using an efficient quantum algorithm for linear equations. We compare the theoretical performance of our quantum algorithm with that of a standard classical algorithm, the conjugate gradient method. We find that the quantum algorithm can achieve a polynomial speedup.

Using various latest cosmological datasets including Type-Ia supernovae, cosmic microwave background radiation, baryon acoustic oscillations, and estimations of the Hubble parameter, we test some dark energy models with parameterized equations of state and try to distinguish or select observation-preferred models. We obtain the best fitting results of the six models and calculate their values of the Akaike Information Criteria and Bayes Information Criterion. And we can distinguish these dark energy models from each other by using these two information criterions. However, the $\Lambda $CDM model remains the best fit model. Furthermore, we perform geometric diagnostics including statefinder and Om diagnostics to understand the geometric behaviour of the dark energy models. We find that the six DE models can be distinguished from each other and from $\Lambda$CDM, Chaplygin gas, quintessence models after the statefinder and Om diagnostics were performed. Finally, we consider the growth factor of the dark energy models with comparison to $\Lambda $CDM model. Still, we find the models can be distinguished from each other and from $\Lambda $CDM model through the growth factor approximation., Comment: 23 pages, 6 figures, to be published in CPL

We have calculated production cross sections of new superheavy elements with atomic number Z=119,120 in the fusion-evaporation reactions of $^{48}$Ca+$^{252}$Es, $^{48}$Ca+$^{257}$Fm, $^{49}$Sc+$^{252}$Es, $^{49}$Sc+$^{251}$Cf, $^{50}$Ti+$^{247}$Bk, $^{50}$Ti+$^{251}$Cf, $^{51}$V+$^{247}$Cm, $^{51}$V+$^{247}$Cf, $^{54}$Cr+$^{243}$Am, $^{54}$Cr+$^{247}$Cm, $^{56}$Mn+$^{244}$Pu, $^{56}$Mn+$^{243}$Am, $^{60}$Fe+$^{237}$Np, $^{60}$Fe+$^{244}$Pu, $^{61}$Co+$^{238}$U, $^{61}$Co+$^{237}$Np, $^{64}$Ni+$^{231}$Pa, $^{64}$Ni+$^{238}$U, $^{65}$Cu+$^{232}$Th, $^{65}$Cu+$^{231}$Pa, and $^{68}$Zn+$^{232}$Th within the dinuclear system model systematically. The inner fusion barriers have been extracted from the driving potential or potential energy surface which could be used to predict the relative fusion probability roughly. The influence of mass asymmetry of the colliding partners on the production of new superheavy elements (SHE) has been investigated systematically. It is found that fusion probability increase along with the increasing mass asymmetry of colliding systems. The Ti-induced reactions have the largest cross-sections of the new SHE. The dependence of production cross-sections of new superheavy elements on the isospin of projectile nuclei has been discussed. The new SHE of $^{289-295}$119 has been predicted as the synthesis cross sections around serval picobarns in the $^{46,48,50,52}$Ti-induced reactions. Production cross-section of the element of $^{295}$120 has been evaluated as large as 1 picobarn in the reactions $^{46}$Ti ($^{251}$Cf, 2n) $^{295}$120 at $E^*$ = 26 MeV. The optimal projectile-target combinations and beam energies for producing new SHE with atomic number Z=119-120 are proposed for the forthcoming experiments., Comment: 12 pages, 7 figures. arXiv admin note: text overlap with arXiv:2210.08941

In the framework of the dinuclear system model, the synthesis mechanism of the superheavy nuclides with atomic number $Z=112, 114, 115$ in the reactions of projectiles $^{40,48}$Ca bombarding on targets $^{238}$U, $^{242}$Pu, and $^{243}$Am at a wide incident energies (excitation energy from 0-100 MeV) have been investigated systematically. Based on the available experimental excitation functions, the dependence of calculated synthesis cross sections on collision orientations has been studied thoroughly. The TKEs of these collisions with the fixed collision orientation show its orientation dependence which can be used to predict the tendency of kinetic energy diffusion. The TKEs are dependent on incident energies which have been discussed. The method of Coulomb barrier distribution function has been applied in our calculations which could treat all of the collision orientations from the tip-tip to side-side approximately. The calculations of excitation functions of $^{48}$Ca + $^{238}$U, $^{48}$Ca + $^{242}$Pu, and $^{48}$Ca + $^{243}$Am have a nice agreement with the available experimental data. The isospin effect of projectiles on production cross sections of moscovium isotopes and the influence of entrance channel effect on the synthesis cross sections of superheavy nuclei have been discussed. The synthesis cross section of new moscovium isotopes $^{278-286}$Mc have been predicted as large as hundreds pb, in the fusion-evaporation reactions of $^{35,37}$Cl + $^{248}$Cf, $^{38,40}$Ar + $^{247}$Bk, $^{39,41}$K + $^{247}$Cm, $^{40,42,44,46}$Ca + $^{238}$Am, $^{45}$Sc + $^{242}$Pu, and $^{46,48,50}$Ti + $^{243}$Np, $^{51}$V + $^{238}$U at the excitation energy interval of 0-100 MeV., Comment: 15 pages, 8 figures, 1 table

We have calculated the multinucleon transfer reactions of $^{208}$Os, $^{208}$Pt, $^{208}$Hg, $^{208}$Pb,$^{208}$Po, $^{208}$Rn, $^{208}$Ra,$^{132,136}$Xe bombarding on $^{232}$Th and $^{248}$Cm at Coulomb barrier energies within the dinuclear system model, systematically. The results are in good agreement with the available experimental data. Coulomb effect and shell effect on production of actinides in these reactions have been investigated thoroughly. Potential energy surface and total kinetic energy mass distributions in the reactions $^{208}$Hg, $^{208}$Pb and$^{208}$Po colliding on $^{248}$Cm and $^{232}$Th are calculated and analyzed, respectively. It is found that PES and TKE spectra manifest the fragment formation mechanism in the multinucleon transfer reactions. The isospin effect and shell effect are shown in PES and TKE. Production cross-sections of multinucleon transfer products are highly dependent on the isobar projectiles with mass number $A=208$. The isobar projectiles with larger N/Z ratios are favorable for creating the neutron-rich target-like fragments. The isobar projectiles with larger charge number induced products tend to shift to proton-rich region. Coulomb potential coupled to shell effect is shown in production cross-sections of actinide isotopes. Based on the radioactive projectiles induced reactions, we have predicted massive new actinide isotopes around nuclear drip lines, even could access the superheavy nuclei region., Comment: 10 pages, 8 figures

We propose a new model for the viscosity of cosmic matters, which can be applied to different epochs of the universe. Using this model, we include the bulk viscosities as practical corrections to the perfect fluid models of the baryonic and dark matters since the material fluids in the real world may have viscosities due to thermodynamics. Such inclusion is put to the test within the framework of $f(T)$ gravity that is proved to be successful in describing the cosmic acceleration, where $T$ denotes the torsion scalar. We perform an observational fit to our model and constrain the cosmological and model parameters by using various latest cosmological datasets. Based on the fitting result, we discuss several cosmological implications including the dissipation of matters, the evolutionary history of the universe, $f(T)$ modification as an effective dark energy, and the Hubble tension problem. The corresponding findings are (i) The late time dissipation will make the density parameters of the matters vanish in the finite future. Moreover, the density ratio between the baryonic and dark matters will change over time. (ii) The radiation dominating era, matter dominating era and the accelerating era can be recovered and the model can successfully describe the known history of the universe. (iii) The $f(T)$ modification is the main drive of the acceleration expansion and currently mimics a phantom-like dark energy. But the universe will eventually enter a de Sitter expansion phase. (iv) The Hubble tension between local and global observations can be significantly alleviated in our model., Comment: accepted by Eur. Phys. J. C

We transform the problem of solving linear system of equations $A\mathbf{x}=\mathbf{b}$ to a problem of finding the right singular vector with singular value zero of an augmented matrix $C$, and present two quantum algorithms for solving this problem. The first algorithm solves the problem directly by applying the quantum eigenstate filtering algorithm with query complexity of $O\left( s\kappa \log \left( 1/\epsilon \right) \right) $ for a $s$-sparse matrix $C$, where $\kappa $ is the condition number of the matrix $A$, and $\epsilon $ is the desired precision. The second algorithm uses the quantum resonant transition approach, the query complexity scales as $O\left[s\kappa + \log\left( 1/\epsilon \right)/\log \log \left( 1/\epsilon \right) \right] $. Both algorithms meet the optimal query complexity in $\kappa $, and are simpler than previous algorithms., Comment: 13 pages, reformulate the algorithms, more materials are added

To enhance the scalability and performance of the traditional finite-difference time-domain (FDTD) methods, a three-dimensional summation-by-parts simultaneous approximation term (SBP-SAT) FDTD method is developed to solve complex electromagnetic problems. It is theoretically stable and can be further used for multiple mesh blocks with different mesh sizes. This paper mainly focuses on the fundamental theoretical aspects upon its three-dimensional implementation, the SAT for various boundary conditions, and the numerical dispersion properties and the comparison with the FDTD method. The proposed SBP-SAT FDTD method inherits all the merits of the FDTD method, which is matrix-free, easy to implement, and has the same level of accuracy with a negligible overhead of runtime (0.13\%) and memory usage (1.2\%). Four numerical examples are carried out to validate the effectiveness of the proposed method.

Total least squares (TLS) methods have been widely used in data fitting. Compared with the least squares method, for TLS problem we takes into account not only the observation errors, but also the errors in the measurement matrix. This is more realistic in practical applications. For the large-scale discrete ill-posed problem $Ax \approx b$, we introduce the quantum-inspired techniques to approximate the truncated total least squares (TTLS) solution. We analyze the accuracy of the quantum-inspired truncated total least squares algorithm and perform numerical experiments to demonstrate the efficiency of our method.

A summation-by-parts simultaneous approximation term (SBP-SAT) finite-difference time-domain (FDTD) subgridding method is proposed to model geometrically fine structures in this paper. Compared with our previous work, the proposed SBP-SAT FDTD method uses the staggered Yee's grid without adding or modifying any field components through field extrapolation on the boundaries to make the discrete operators satisfy the SBP property. The accuracy of extrapolation keeps consistency with that of the second-order finite-difference scheme near the boundaries. In addition, the SATs are used to weakly enforce the tangential boundary conditions between multiple mesh blocks with different mesh sizes. With carefully designed interpolation matrices and selected free parameters of the SATs, no dissipation occurs in the whole computational domain. Therefore, its long-time stability is theoretically guaranteed. Three numerical examples are carried out to validate its effectiveness. Results show that the proposed SBP-SAT FDTD subgridding method is stable, accurate, efficient, and easy to implement based on existing FDTD codes with only a few modifications., Comment: arXiv admin note: text overlap with arXiv:2110.09054

Within the framework of the dinuclear system model, the reaction mechanisms for synthesizing target-like isotopes from Bk to compound nuclei Lv are thoroughly investigated in complete and incomplete fusion reaction of $^{48}$Ca +$^{248}$Cm around Coulomb barrier energies. Production cross-section of $^{292,293}$Lv as a function of excitation energy in fusion-evaporation reactions and target-like isotopic yields in multinucleon transfer reactions are evaluated, in which a statistical approach is used to describe the decay process of excited nuclei. The available experimental data can be reproduced well with the model reasonably. The products of all possible formed isotopes in the dynamical pre-equilibrium process for collision partners at incident energy $E_{\rm lab}$ = 5.5 MeV/nucleon are exported, systematically. It is found that the quasi-fission fragments are dominant in the yields. The optimal pathway from the target to compound nuclei shows up along the valley of potential surface energy. The effective impact parameter of two colliding partners leading to compound nuclei is selected from head-on collision to semi-central collision with $L$ = 52 $\rm \hbar$. The timescale boundary between complete fusion and multinucleon transfer reactions is about 5.7$\times 10^{-21}$ s with effective impact parameters. Synthesis cross-section of unknown neutron-rich actinides from Bk to Rf have been predicted around several nanobarns.