skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on March 1, 2026

Title: An age-structured syphilis model, II: optimal control and numerical simulation
The World Health Organization estimated that 8 million adults between 15 and 49 years old acquired syphilis globally in 2022. China CDC reported that there were 530 116 cases of syphilis in mainland China in 2023. Since syphilis is a sexually transmitted disease and age structure of the host population plays a crucial role, in this series of two papers we develop an age-structured model with four infection stages (primary, secondary, latent and tertiary) to study the transmission dynamics of syphilis. In part I (Wuet al. 2025Proc. R. Soc. A481: 20240218 (doi:10.1098/rspa.2024.0218)), we investigated the well-posedness of the model and studied stability of the steady states. In part II, first, we consider the optimal control of the age-structured model. Second, utilizing the Markov Chain Monte Carlo method, we calibrate the reported syphilis data in China by using a demographic model. Finally, we apply the relevant simulation results to numerically simulate the age-structured model. Our results indicate that (i) for the syphilis demographic model, the basic reproduction number R 0 2.4876 with CI (95%) (1.6823, 3.1434); (ii) tertiary stage infection is more severe in the elderly population; (iii) reducing the number of secondary and latent stage syphilis individuals can effectively reduce the total number of infected populations.  more » « less
Award ID(s):
2424605 2052648
PAR ID:
10618150
Author(s) / Creator(s):
;
Publisher / Repository:
Ruan, Shigui
Date Published:
Journal Name:
Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences
Volume:
481
Issue:
2310
ISSN:
1471-2946
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract We consider the problem of covering multiple submodular constraints. Given a finite ground setN, a weight function$$w: N \rightarrow \mathbb {R}_+$$ w : N R + ,rmonotone submodular functions$$f_1,f_2,\ldots ,f_r$$ f 1 , f 2 , , f r overNand requirements$$k_1,k_2,\ldots ,k_r$$ k 1 , k 2 , , k r the goal is to find a minimum weight subset$$S \subseteq N$$ S N such that$$f_i(S) \ge k_i$$ f i ( S ) k i for$$1 \le i \le r$$ 1 i r . We refer to this problem asMulti-Submod-Coverand it was recently considered by Har-Peled and Jones (Few cuts meet many point sets. CoRR.arxiv:abs1808.03260Har-Peled and Jones 2018) who were motivated by an application in geometry. Even with$$r=1$$ r = 1 Multi-Submod-Covergeneralizes the well-known Submodular Set Cover problem (Submod-SC), and it can also be easily reduced toSubmod-SC. A simple greedy algorithm gives an$$O(\log (kr))$$ O ( log ( k r ) ) approximation where$$k = \sum _i k_i$$ k = i k i and this ratio cannot be improved in the general case. In this paper, motivated by several concrete applications, we consider two ways to improve upon the approximation given by the greedy algorithm. First, we give a bicriteria approximation algorithm forMulti-Submod-Coverthat covers each constraint to within a factor of$$(1-1/e-\varepsilon )$$ ( 1 - 1 / e - ε ) while incurring an approximation of$$O(\frac{1}{\epsilon }\log r)$$ O ( 1 ϵ log r ) in the cost. Second, we consider the special case when each$$f_i$$ f i is a obtained from a truncated coverage function and obtain an algorithm that generalizes previous work on partial set cover (Partial-SC), covering integer programs (CIPs) and multiple vertex cover constraints Bera et al. (Theoret Comput Sci 555:2–8 Bera et al. 2014). Both these algorithms are based on mathematical programming relaxations that avoid the limitations of the greedy algorithm. We demonstrate the implications of our algorithms and related ideas to several applications ranging from geometric covering problems to clustering with outliers. Our work highlights the utility of the high-level model and the lens of submodularity in addressing this class of covering problems. 
    more » « less
  2. Abstract We conduct 3D magnetohydrodynamic simulations of decaying turbulence in the context of the solar wind. To account for the spherical expansion of the solar wind, we implement the expanding box model. The initial turbulence comprises uncorrelated counterpropagating Alfvén waves and exhibits an isotropic power spectrum. Our findings reveal the consistent generation of negative residual energy whenever nonlinear interactions are present, independent of the normalized cross helicityσcand compressibility. The spherical expansion facilitates this process. The resulting residual energy is primarily distributed in the perpendicular direction, withS2(b) − S2(u) ∝ lor equivalently E r k 2 . HereS2(b) andS2(u) are second-order structure functions of magnetic field and velocity respectively. In most runs,S2(b) develops a scaling relation S 2 ( b ) l 1 / 2 ( E b k 3 / 2 ). In contrast,S2(u) is consistently shallower thanS2(b), which aligns with in situ observations of the solar wind. We observe that the higher-order statistics of the turbulence, which act as a proxy for intermittency, depend on the initialσcand are strongly affected by the expansion effect. Generally, the intermittency is more pronounced when the expansion effect is present. Finally, we find that in our simulations, although the negative residual energy and intermittency grow simultaneously as the turbulence evolves, the causal relation between them seems to be weak, possibly because they are generated on different scales. 
    more » « less
  3. Abstract Let$$p_{1},\ldots ,p_{n}$$ p 1 , , p n be a set of points in the unit square and let$$T_{1},\ldots ,T_{n}$$ T 1 , , T n be a set of$$\delta $$ δ -tubes such that$$T_{j}$$ T j passes through$$p_{j}$$ p j . We prove a lower bound for the number of incidences between the points and tubes under a natural regularity condition (similar to Frostman regularity). As a consequence, we show that in any configuration of points$$p_{1},\ldots , p_{n} \in [0,1]^{2}$$ p 1 , , p n [ 0 , 1 ] 2 along with a line$$\ell _{j}$$ j through each point$$p_{j}$$ p j , there exist$$j\neq k$$ j k for which$$d(p_{j}, \ell _{k}) \lesssim n^{-2/3+o(1)}$$ d ( p j , k ) n 2 / 3 + o ( 1 ) . It follows from the latter result that any set of$$n$$ n points in the unit square contains three points forming a triangle of area at most$$n^{-7/6+o(1)}$$ n 7 / 6 + o ( 1 ) . This new upper bound for Heilbronn’s triangle problem attains the high-low limit established in our previous work arXiv:2305.18253. 
    more » « less
  4. Abstract We studyℓnorms ofℓ2-normalized eigenfunctions of quantum cat maps. For maps with short quantum periods (constructed by Bonechi and de Biévre in F Bonechi and S De Bièvre (2000,Communications in Mathematical Physics,211, 659–686)) we show that there exists a sequence of eigenfunctionsuwith u log N 1 / 2 . For general eigenfunctions we show the upper bound u log N 1 / 2 . Here the semiclassical parameter is h = 2 π N 1 . Our upper bound is analogous to the one proved by Bérard in P Bérard (1977,Mathematische Zeitschrift,155, 249-276) for compact Riemannian manifolds without conjugate points. 
    more » « less
  5. Abstract A steady-state, semi-analytical model of energetic particle acceleration in radio-jet shear flows due to cosmic-ray viscosity obtained by Webb et al. is generalized to take into account more general cosmic-ray boundary spectra. This involves solving a mixed Dirichlet–Von Neumann boundary value problem at the edge of the jet. The energetic particle distribution functionf0(r,p) at cylindrical radiusrfrom the jet axis (assumed to lie along thez-axis) is given by convolving the particle momentum spectrum f 0 ( , p ) with the Green’s function G ( r , p ; p ) , which describes the monoenergetic spectrum solution in which f 0 δ ( p p ) asr→ ∞ . Previous work by Webb et al. studied only the Green’s function solution for G ( r , p ; p ) . In this paper, we explore for the first time, solutions for more general and realistic forms for f 0 ( , p ) . The flow velocityu=u(r)ezis along the axis of the jet (thez-axis).uis independent ofz, andu(r) is a monotonic decreasing function ofr. The scattering time τ ( r , p ) = τ 0 ( p / p 0 ) α in the shear flow region 0 <r<r2, and τ ( r , p ) = τ 0 ( p / p 0 ) α ( r / r 2 ) s , wheres> 0 in the regionr>r2is outside the jet. Other original aspects of the analysis are (i) the use of cosmic ray flow lines in (r,p) space to clarify the particle spatial transport and momentum changes and (ii) the determination of the probability distribution ψ p ( r , p ; p ) that particles observed at (r,p) originated fromr→ ∞ with momentum p . The acceleration of ultrahigh-energy cosmic rays in active galactic nuclei jet sources is discussed. Leaky box models for electron acceleration are described. 
    more » « less