Aberrant intrinsic functional brain topology in methamphetamine-dependent individuals after six-months of abstinence

Xiang Li; Jinyu Cong; Kunmeng Liu; Pingping Wang; Min Sun; Benzheng Wei; Xiang Li; Jinyu Cong; Kunmeng Liu; Pingping Wang; Min Sun; Benzheng Wei

doi:10.3934/mbe.2023867

Mathematical Biosciences and Engineering

2023, Volume 20, Issue 11: 19565-19583. doi: 10.3934/mbe.2023867

Previous Article Next Article

Research article Special Issues

Aberrant intrinsic functional brain topology in methamphetamine-dependent individuals after six-months of abstinence

1.
First Clinical Medical College, Shandong University of Traditional Chinese Medicine, Jinan 250355, China
2.
Center for Medical Artificial Intelligence, Shandong University of Traditional Chinese Medicine, Qingdao 266112, China
3.
Qingdao Academy of Chinese Medical Sciences, Shandong University of Traditional Chinese Medicine, Qingdao 266112, China
4.
Shandong Detoxification Monitoring and Treatment Institute, Zibo 255311, China

Academic Editor: Shangce Gao

Received: 10 August 2023 Revised: 05 October 2023 Accepted: 16 October 2023 Published: 25 October 2023

Our aim was to explore the aberrant intrinsic functional topology in methamphetamine-dependent individuals after six months of abstinence using resting-state functional magnetic imaging (rs-fMRI). Eleven methamphetamines (MA) abstainers who have abstained for six months and eleven healthy controls (HC) were recruited for rs-fMRI examination. The graph theory and functional connectivity (FC) analysis were employed to investigate the aberrant intrinsic functional brain topology between the two groups at multiple levels. Compared with the HC group, the characteristic shortest path length ( ${L}_{p}$ ) showed a significant decrease at the global level, while the global efficiency ( ${E}_{glob}$ ) and local efficiency ( ${E}_{loc}$ ) showed an increase considerably. After FDR correction, we found significant group differences in nodal degree and nodal efficiency at the regional level in the ventral attentional network (VAN), dorsal attentional network (DAN), somatosensory network (SMN), visual network (VN) and default mode network (DMN). In addition, the NBS method presented the aberrations in edge-based FC, including frontoparietal network (FPN), subcortical network (SCN), VAN, DAN, SMN, VN and DMN. Moreover, the FC of large-scale functional brain networks revealed a decrease within the VN and SCN and between the networks. These findings suggest that some functions, e.g., visual processing skills, object recognition and memory, may not fully recover after six months of withdrawal. This leads to the possibility of relapse behavior when confronted with MA-related cues, which may contribute to explaining the relapse mechanism. We also provide an imaging basis for revealing the neural mechanism of MA-dependency after six months of abstinence.

Keywords:

Citation: Xiang Li, Jinyu Cong, Kunmeng Liu, Pingping Wang, Min Sun, Benzheng Wei. Aberrant intrinsic functional brain topology in methamphetamine-dependent individuals after six-months of abstinence[J]. Mathematical Biosciences and Engineering, 2023, 20(11): 19565-19583. doi: 10.3934/mbe.2023867

Related Papers:

[1]	Annalisa Quaini, Sunčica Čanić, David Paniagua . Numerical characterization of hemodynamics conditions near aortic valve after implantation of left ventricular assist device. Mathematical Biosciences and Engineering, 2011, 8(3): 785-806. doi: 10.3934/mbe.2011.8.785
[2]	Lorenzo Civilla, Agnese Sbrollini, Laura Burattini, Micaela Morettini . An integrated lumped-parameter model of the cardiovascular system for the simulation of acute ischemic stroke: description of instantaneous changes in hemodynamics. Mathematical Biosciences and Engineering, 2021, 18(4): 3993-4010. doi: 10.3934/mbe.2021200
[3]	Li Cai, Yu Hao, Pengfei Ma, Guangyu Zhu, Xiaoyu Luo, Hao Gao . Fluid-structure interaction simulation of calcified aortic valve stenosis. Mathematical Biosciences and Engineering, 2022, 19(12): 13172-13192. doi: 10.3934/mbe.2022616
[4]	Scott R. Pope, Laura M. Ellwein, Cheryl L. Zapata, Vera Novak, C. T. Kelley, Mette S. Olufsen . Estimation and identification of parameters in a lumped cerebrovascular model. Mathematical Biosciences and Engineering, 2009, 6(1): 93-115. doi: 10.3934/mbe.2009.6.93
[5]	Xinyu Zhang, Yu Meng, Mei Jiang, Lin Yang, Kuixing Zhang, Cuiting Lian, Ziwei Li . Machine learning-based evaluation of application value of pulse wave parameter model in the diagnosis of hypertensive disorder in pregnancy. Mathematical Biosciences and Engineering, 2023, 20(5): 8308-8319. doi: 10.3934/mbe.2023363
[6]	Yuqian Mei, Debao Guan, Xinyu Tong, Qian Liu, Mingcheng Hu, Guangxin Chen, Caijuan Li . Association of cerebral infarction with vertebral arterial fenestration using non-Newtonian hemodynamic evaluation. Mathematical Biosciences and Engineering, 2022, 19(7): 7076-7090. doi: 10.3934/mbe.2022334
[7]	Li Cai, Qian Zhong, Juan Xu, Yuan Huang, Hao Gao . A lumped parameter model for evaluating coronary artery blood supply capacity. Mathematical Biosciences and Engineering, 2024, 21(4): 5838-5862. doi: 10.3934/mbe.2024258
[8]	Lijian Xu, Lekang Yin, Youjun Liu, Fuyou Liang . A computational study on the influence of aortic valve disease on hemodynamics in dilated aorta. Mathematical Biosciences and Engineering, 2020, 17(1): 606-626. doi: 10.3934/mbe.2020031
[9]	Giovanna Guidoboni, Alon Harris, Lucia Carichino, Yoel Arieli, Brent A. Siesky . Effect of intraocular pressure on the hemodynamics of the central retinal artery: A mathematical model. Mathematical Biosciences and Engineering, 2014, 11(3): 523-546. doi: 10.3934/mbe.2014.11.523
[10]	Rachel Clipp, Brooke Steele . An evaluation of dynamic outlet boundary conditions in a 1D fluid dynamics model. Mathematical Biosciences and Engineering, 2012, 9(1): 61-74. doi: 10.3934/mbe.2012.9.61

Abstract

1. Introduction

The cardiovascular system is a complex dynamic system, which contains rich dynamics information to predict the physiological or pathological state. The traditional method for studying cardiovascular hemodynamics is to perform surgery on animals or humans, which is usually limited by ethics and experimental conditions ^[1,2]. As the method and tool in vitro, the numerical model and mock circulatory model can provide technical support for clinical noninvasive diagnosis ^[3]. The numerical modeling is the most widespread research technique due to its proven capability, efficiency and minimized cost in comparison to in-vitro studies in mock flow loop, or in-vivo studies in animals and humans ^[4].

The cardiovascular lumped parameter (CLP) model is a simplified numerical method to study hemodynamics. Researchers have done a lot of work in the CLP modeling in the past ten years. Some lumped parameter models for a single circulation loop were modeled ^[5,6,7]. For example, Figliola et al. ^[5] modeled a lumped parameter model including only right heart and pulmonary circulation. Simaan et al. ^[7] modeled a fifth-order lumped parameter circuit model for left heart and systemic circulation. These simplified Windkessel models can be used to simulate the hemodynamics of a single circulation loop. Some complete CLP models have been developed ^{[8,9,10,11,12,13]}. For example, a CLP model including systemic circulation, pulmonary circulation, and the heart with four chambers was proposed ^[8]. To simulate the hemodynamics of patients with heart failure, Gu et al. ^[9] established the CLP for patients with heart failure, and verified the effectiveness of physiological mechanism by the model. Abdi ^[10] used a lumped method to institute a cardiovascular model with 42 chambers, including arteries, veins, capillaries, and heart chambers. This model was used to evaluate the effect of the abnormal heartbeat on cardiovascular system performance. Although these models are capable of simulating the hemodynamics of the intact cardiovascular system, they are not combined with ventricular assist devices.

Since their first surgical implantation in the 1960s ^[14], ventricular assist devices play an active role in the treatment of heart failure. The new generation of LVADs use axial or centrifugal rotary pumps ^[15]. Compared with the first generation of pulsating pumps, the rotary pumps have the advantages of small size, high reliability, and durability. It is helpful for clinical management to comprehend the coupling relationship between the native heart and the rotary pump. Coupling modeling of the cardiovascular system and ventricular assist device (VAD) is a feasible way to explain the hemodynamic relationship between a native heart and a rotary pump. Shi et al. ^[4,14] proposed a coupled model for the cardiovascular system and VAD to simulate the hemodynamics under left ventricular (LV) failure with VAD support. CARDIOSIM^© is a developed simulator based on the CLP model by Lazzari ^{[15,16,17,18]}. It not only can be used to simulate the physiological and pathological hemodynamics, but also can be combined with various types of ventricular assist devices (pulsating pumps, rotary pumps, etc.). CARDIOSIM^© has also been applied in medical education ^[19]. Lim et al. ^[1] developed a coupled model between the cardiovascular system and implantable rotary blood pump, and it includes the left heart, right heart, systemic circulation, pulmonary circulation, and implantable rotary blood pump. Di et al. ^[20,21,22] established a complete CLP model, which could be connected with pulsating pumps or rotary pumps. In the case of cardiac dysplasia or heart failure, it is helpful to understand the effects of single or two VADs on hemodynamics. However, the autonomic regulation of the cardiovascular system is not considered in the above models.

Baroreflex is a kind of short-term autonomic nerve regulation, which plays a great role in short-term pressure control ^[23]. The work of Ursino ^[24] and Ottesen ^[25] indicates that the baroreflex system can maintain the stability of the hemodynamic environment by regulating neural effectors. Consequently, it is more autonomous and accurate for simulating hemodynamics that a combination of the CLP model and the baroreflex model to form a closed-loop control system. Certainly, baroreflex regulation was also considered in some coupled models for the cardiovascular system and VADs ^{[26,27,28,29,30,31,32]}. In these models, ventricular characteristics or systemic vascular resistance are usually altered to simulate the hemodynamics of the cardiovascular system with impaired heart function. Bozkurt et al. ^[27] reduced left ventricular systolic elastance to 0.5 mmHg/mL from 2.5 mmHg/mL to simulate a dilated cardiomyopathy (DCM) condition in the cardiovascular system model. Cox et al. ^[28] changed left ventricular contractility, initial left ventricular volume and systemic vascular resistance to simulate end-stage DCM. Lim et al. ^[32] decreased the end-systolic elastance value to 40% of normal to implement a heart-failure model. Heart failure is a complex syndrome, it is not only related to the cardiovascular system but also the baroreflex system. These models only consider ventricular characteristics or systemic vascular resistance to simulate heart failure, which is different from clinical observations. However, other hemodynamic parameters in the cardiovascular system, such as systemic vascular compliance and pulmonary vascular resistance, as well as some parameters in the baroreflex system, should also be considered. Otherwise, some changes in hemodynamic parameters during LVAD support and the prediction results of the LVAD operating condition will be affected. The application scope of the above model is limited. Therefore, it is necessary to develop a comprehensive numerical model and consider the mechanism of changes in the cardiovascular system and the baroreflex system to improve the simulation and prediction effect.

In this study, we aim to establish a CLP model coupled with the baroreflex model and the LVAD model to form a comprehensive numerical model, which provides a platform for the hemodynamics simulation and the noninvasive evaluation of LVAD operating condition.

This paper is organized in the following manner: The modeling of the coupling cardiovascular system is introduced in Section 2. Section 3 presents some results including the hemodynamics under different conditions and the hemodynamic variables at different pump speeds. Section 4 discusses the operating condition of HeartMate Ⅲ and its effect on hemodynamics. In addition, the limitations of this study are also described in Section 4. Lastly, conclusions are drawn on the outcomes of this study in Section 5.

2. Materials and method

The comprehensive numerical model is illustrated in Figure 1. The numerical model was made up of three parts: a developed CLP model for describing the human cardiovascular circulation system, a baroreflex model for short-term neuroregulation, and an LVAD model for assisting diseased heart. The rotary pump was connected with the diseased heart in parallel. Thus, the pump inlet cannulae was connected to the left ventricular apex and the outlet cannulae was connected to the aorta.

Figure 1. Block diagram of numerical model.

$R$ ,

$L$ and

$C$ denote resistance, inertance and compliance,

$P$ and

$Q$ denote pressure and flow rate,

$D$ denotes ideal diodes,

$H$ is pressure difference between the inlet and outlet of the rotary pump,

$t$ and

$T$ are time variable and heart period,

${E_{\max }}$ denotes ventricular end-systolic elastance, subscripts

${\rm{MV}}$ ,

${\rm{AV}}$ ,

${\rm{TV}}$ and

${\rm{PV}}$ are mitral, aortic, tricuspid and pulmonary valves,

$lv$ ,

$la$ ,

$rv$ and

$ra$ denote left ventricle, left atrium, right ventricle and right atrium,

$ao$ and

$pa$ denote aorta and pulmonary artery,

$sar$ and

$par$ denote systemic artery and pulmonary arteriole,

$svr$ and

$pvr$ denote systemic vascular resistance and pulmonary vascular resistance,

$sv$ and

$pv$ denote systemic and pulmonary veins,

$irp$ is implantable rotary pump.

DownLoad: Full-Size Img PowerPoint

2.1. CLP Model

The lumped parameter method is to concentrate on the human cardiovascular circulation system into several chambers with organ and tissue functions. Each chamber is composed of Windkessel elements representing hemodynamic parameters. This method can not only reduce the number of parameters to improve computational efficiency, but also obtain global hemodynamics. In this paper, the cardiovascular circulation system was divided into four parts: left heart, right heart, systemic circulation, and pulmonary circulation. The heart is a four-chamber pump with time-varying elasticity. The direction of blood flow is controlled by the four valves of the heart. Each chamber of heart was modeled as a Windkessel model.

2.1.1. Ventricles

Suga et al. ^[33] show that the pressure-volume relationship of the ventricle can be expressed by a time-varying ventricle elastance. The ventricle elastance characterizes the elastic of ventricular myocardium and describes the dynamics of ventricular systole and diastole. Numerically, the ventricle elastance is generally defined as the ratio of pressure and volume in the ventricle. In the research field of the cardiovascular model, the ventricle elastance has been recognized and used ^{[1,4,7,11,13,18,21]}.

The time-varying elastances of the left and right ventricles:

${E_{lv}}(t) = \frac{1}{{{C_{lv}}(t)}} = \frac{{{P_{lv}}(t)}}{{{V_{lv}}(t) - {V_{lv, 0}}}}$

(1)

${E_{rv}}(t) = \frac{1}{{{C_{rv}}\;(t)}} = \frac{{{P_{rv}}\;(t)}}{{{V_{rv}}\;(t) - {V_{rv, 0}}}}$

(2)

where ${E_{lv}}(t)$ and ${E_{rv}}(t)$ are the time-varying elastances of left and right ventricles, ${C_{lv}}(t)$ and ${C_{rv}}(t)$ are the time-varying compliances of left and right ventricles, which are the reciprocals of ${E_{lv}}(t)$ and ${E_{rv}}(t)$ , ${P_{lv}}(t)$ and ${P_{rv}}(t)$ are pressures of left and right ventricles, ${V_{lv}}(t)$ and ${V_{rv}}(t)$ are volumes of left and right ventricles, ${V_{lv, 0}}$ and ${V_{rv, 0}}$ are theoretical volumes in the left and right ventricles at zero pressure.

However, when the pressure and volume of the ventricle are not be provided, the time-varying elastance cannot be determined. Accordingly, the calculation formulas commonly used in engineering are as follows ^[7]:

${E_{lv}}(t) = ({E_{\max , lv}} - {E_{\min , lv}}){E_{n, lv}}({t_n}) + {E_{\min , lv}}$

(3)

${E_{rv}}(t) = ({E_{\max , rv}} - {E_{\min , rv}}){E_{n, rv}}({t_n}) + {E_{\min , rv}}$

(4)

where ${E_{n, lv}}({t_n})$ and ${E_{n, rv}}({t_n})$ are the normalized elastic functions. The actions of the left and right ventricles are synchronous, therefore, the same normalized elastic function is used for the left ventricle and right ventricle.

${E_{n,lv}}({t_n}) = {E_{n,rv}}({t_n}) = 1.55 \times [\frac{{{{(\frac{{{t_n}}}{{0.7}})}^{1.9}}}}{{1 + {{(\frac{{{t_n}}}{{0.7}}\;)}^{1.9}}}}\;] \times [\frac{1}{{1 + {{(\frac{{{t_n}}}{{1.17}})}^{21.9}}\;\;}}]$

(5)

${t_n} = \frac{t}{{0.2 + 0.15T}}$

(6)

$T = 60/{\rm{HR}}$

(7)

where $t$ and $T$ are time variable and heart period, ${\rm{HR}}$ is heart rate, ${E_{\max, lv}}$ and ${E_{\min, lv}}$ are left ventricular end-systolic elastance and left ventricular end-diastolic elastance, ${E_{\max, rv}}$ and ${E_{\min, rv}}$ are right ventricular end-systolic elastance and right ventricular end-diastolic elastance.

2.1.2. Atria

The atrium also contracts and relaxes periodically similar to the ventricle, but its contractility is much smaller than that of the ventricle. The main function of the atrium is to store venous blood. Similarly, the time-varying elastance was adopted to describe the pressure-volume relationship of the atrium. The states of the left and right atria are modeled by the time-varying compliances ${C_{la}}(t)$ and ${C_{ra}}(t)$ . ${E_{la}}(t)$ and ${E_{ra}}(t)$ are the time-varying elastances of left and right atria. ${E_{la}}(t)$ and ${E_{ra}}(t)$ can be expressed by Eqs (8) and (9). The meanings of other parameters have resemblance with those for ventricles.

${E_{la}}(t) = \frac{1}{{{C_{la}}(t)}} = \frac{{{P_{la}}(t)}}{{{V_{la}}(t) - {V_{la, 0}}}}$

(8)

${E_{ra}}(t) = \frac{1}{{{C_{ra}}(t)}} = \frac{{{P_{ra}}(t)}}{{{V_{ra}}(t) - {V_{ra, 0}}}}$

(9)

Although the activities of the atrium and ventricle are asynchronous, the formulas for calculating atrial elastances are similar to those for the ventricles. The expressions of ${E_{la}}(t)$ and ${E_{ra}}(t)$ are as follows ^[14]:

${E_{la}}(t) = \frac{{({E_{\max ,la}}\;\; - {E_{\min ,la}}\;\;)}}{2}{e_{la}}(t) + {E_{\min ,la}}$

(10)

${E_{ra}}(t) = \frac{{({E_{\max ,ra}}\;\; - {E_{\min ,ra}}\;\;)}}{2}{e_{ra}}(t) + {E_{\min ,ra}}$

(11)

where ${E_{\max, la}}$ and ${E_{\min, la}}$ are left atrial end-systolic elastance and left atrial end-diastolic elastance, ${E_{\max, ra}}$ and ${E_{\min, ra}}$ are right atrial end-systolic elastance and right atrial end-diastolic elastance, ${e_{la}}(t)$ and ${e_{ra}}(t)$ are the activation functions of the left and right atria.

The synchronous property of left and right atria gives rise to the same activation functions for left and right atria.

${e_{la}}(t) = {e_{ra}}(t) = \left\{ \begin{gathered} 0, \quad \quad \quad \quad \quad \, \, \, \, \;\;\, 0 \leqslant t \lt {T_{asb}} \\ 1 - \cos (\frac{{t - {T_{asb}}}}{{{T_{ase}}{\rm{ - }}{T_{asb}}}}2\pi ), {T_{asb}} \leqslant t \lt {T_{ase}} \\ 0, \quad \quad \quad \quad \quad \, \, \, \, \;\;\, {T_{ase}} \leqslant t \lt T \\ \end{gathered} \right.$

(12)

where ${T_{asb}}$ is the beginning time of atrium contraction and ${T_{asb}} = 0.9T$ , ${T_{ase}}$ is the end time of atrium contraction and ${T_{ase}} = 0.09T + {T_{asb}}$ .

2.1.3. Valves

The heart valves prevent blood from flowing back, it comprises the mitral and tricuspid valves between the atria and ventricles, as well as the aortic and pulmonary valves between the ventricles and arteries. The opening and closing of the valve are similar to the short circuit and open circuit of an ideal diode. Consequently, the circuit model for the valve can be represented by resistor in series ideal diode. Take the aortic valve in for example. The circuit model of the aortic valve is composed of aortic valve resistance ${R_{{\rm{AV}}}}$ and ideal diode ${D_{{\rm{AV}}}}$ , and the opening or closing of the valve is controlled by the difference between the left ventricular pressure ${P_{lv}}$ and the aortic pressure ${P_{ao}}$ . If ${P_{lv}}{\rm{ - }}{P_{ao}} > 0$ , the aortic valve is opened and the flow rate through the valve is ${{({P_{lv}} - {P_{ao}})} / {{R_{{\rm{AV}}}}}}$ , otherwise the aortic valve is closed and the flow rate is 0.

The states of the valve are determined by the transvalvular pressure difference. The function of transvalvular pressure difference $\gamma (\xi)$ is introduced:

$\gamma (\xi ) = \left\{ \begin{gathered} \xi , \quad {\rm{if}} \quad\xi \gt 0 \\ 0, \quad {\rm{if}} \quad \xi \leqslant 0 \\ \end{gathered} \right.$

(13)

The valve flow rate ${Q_V}$ is as follow:

${Q_V} = \frac{1}{{{R_V}}}\gamma (\xi )$

(14)

where $\xi$ is the transvalvular pressure difference, ${Q_V}$ is the flow rate through the valve and the $V$ is the abbreviation of the valve.

2.1.4. Systemic circulation and pulmonary circulation

The systemic circulation is separated into aorta, artery, arteriole, capillary, and vein segments. The aorta was regarded as a vascular chamber with resistance, inertance, and compliance. The artery segment was considered as a vascular chamber with resistance and compliance. The primary role of the arteriole and capillary is to provide resistance to blood flow. Arteriole and capillary were modeled as a vascular chamber with resistance. The vascular resistance is provided by the artery, arteriole, and capillary. All levels of veins were modeled as a vascular chamber. Collecting and storing blood are primary roles played by veins that are possessed of greater compliance. The model for pulmonary circulation is similar to that for systemic circulation.

There is an analogy between fluid network and electric network. An electrical analogy is typically used where pressure is analogous to voltage, flow rate is analogous to electric current, resistance is analogous to resistor, inertance is analogous to inductor, and compliance is analogous to capacitor ^[12,34]. On account of an electrical network is applied to the blood circulation system, this paper used the Kirchhoff's law to expound the coupling relationship between the branches of the blood circulation system. Using the notation given in Figure 1, the equations of the instantaneous pressure and flow rate for systemic circulation are given as follows.

The flow rate balance equations are:

${C_{ao}} \cdot \frac{{d{P_{ao}}(t)}}{{dt}} = \frac{{\gamma ({P_{lv}}(t) - {P_{ao}}(t))}}{{{R_{{\rm{AV}}}}}} - {Q_{ao}}(t) + {Q_{irp}}$

(15)

${C_{sar}} \cdot \frac{{d{P_{sar}}(t)}}{{dt}} = {Q_{ao}}(t) - \frac{{{P_{sar}}(t) - {P_{sv}}(t)}}{{{R_{svr}}(t)}}$

(16)

${C_{sv}} \cdot \frac{{d{P_{sv}}(t)}}{{dt}} = \frac{{{P_{sar}}(t) - {P_{sv}}(t)}}{{{R_{svr}}(t)}} - \frac{{{P_{sv}}(t) - {P_{ra}}(t)}}{{{R_{sv}}}}$

(17)

The pressure balance equation is:

${L_{ao}} \cdot \frac{{d{Q_{ao}}(t)}}{{dt}} = {P_{ao}}(t) - {R_{ao}} \cdot {Q_{ao}}(t) - {P_{sar}}(t)$

(18)

2.2. Baroreflex model

According to Ursino et al. ^[24], the baroreflex model was conducted in the order of afferent neural pathway, efferent neural pathways and effectors. In this study, the baroreflex model was modified to control the corresponding effectors in the cardiovascular system: systemic vascular resistance, left ventricular end-systolic elastance, right ventricular end-systolic elastance, and heart period. The open-loop control system for baroreflex regulation is shown in Figure 2.

Figure 2. Block diagram describing the baroreflex model in complex domain.

DownLoad: Full-Size Img PowerPoint

2.2.1. Afferent neural pathway

Ursino et al. ^[24] described the relationship between arterial pressure and the activity of the afferent nerves by considering two main aspects arising from experiments. First, the relationship between arterial pressure and afferent nerve activity in static conditions exhibits a sigmoidal shape with upper saturation and lower threshold. Second, the arterial baroreceptors are sensitive not only to the arterial pressure but also to its rate of change. These two properties have been reproduced by representing the arterial baroreceptors plus the afferent neural pathway as a first-order linear differential equation, with a static gain and a rate-dependent gain, connected in series with a sigmoidal static function.

${\tau _p} \cdot \frac{{d\tilde p(t)}}{{dt}} + \tilde p(t) = {\tau _z} \cdot \frac{{d{p_{ao}}(t)}}{{dt}} + {p_{ao}}(t)$

(19)

${f_{as}}(t) = \frac{{{f_{\min }}\; + {f_{\max }}\; \cdot {e^{[\tilde p\;(t) - {p_n}\;]/{k_a}}}}}{{1 + {e^{[\tilde p(t) - {p_{n\;}}]/{k_a}}}}}$

(20)

where ${\tau _p}$ and ${\tau _z}$ are the time constants, aortic pressure ${p_{ao}}$ is used as the input pressure to arterial baroreceptors, $\tilde p$ is the output variable of a first-order linear differentiator, ${f_{as}}$ is the frequency of spikes in the afferent fibers, ${f_{\min }}$ and ${f_{\max }}$ are the lower and upper saturation of the frequency discharge, ${p_n}$ is the pressure value at the central point of the sigmoidal functional, ${k_a}$ is a parameter with the dimension of pressure, related to the slope of the static function at the central point, and ${1 /{{k_a}}}$ denotes the baroreceptor sensitivity.

2.2.2. Efferent neural pathways

The efferent pathway includes the efferent sympathetic pathway and efferent vagal pathway. Increasing the frequency of spikes in the afferent fibers results in a decrease in the frequency of the sympathetic fibers and in an increase in vagal activity. These relationships have been reproduced by using a monotonically decreasing monoexponential static curve for the sympathetic efferent activity and a monotonically increasing exponential curve, with upper saturation, for the vagal activity.

${f_{es}}(t) = {f_{es, \infty }} + ({f_{es, 0}} - {f_{es, \infty }}) \cdot {e^{ - {k_{es}} \cdot {f_{as}}(t)}}$

(21)

${f_{ev}}(t) = \frac{{{f_{ev,0}} + {f_{ev,\infty }}\; \cdot \;{e^{[{f_{as}}\;(t) - {f_{as,0}}\;]/{k_{ev}}}}}}{{1 + {e^{[{f_{as}}(t) - {f_{as,0}}\;]/{k_{ev}}}}}}$

(22)

where ${f_{es}}$ is the frequency of spikes in the efferent sympathetic nerves, ${f_{es, \infty }}$ , ${f_{es, 0}}$ and ${k_{es}}$ are constants with ${f_{es, 0}} > {f_{es, \infty }}$ , ${f_{ev}}$ is the frequency of spikes in the efferent vagal fibers, ${f_{ev, \infty }}$ , ${f_{ev, 0}}$ and ${k_{ev}}$ are constant parameters with ${f_{ev, 0}} > {f_{ev, \infty }}$ , and ${f_{as, 0}}$ is the central value in Eq 20.

2.2.3. Effectors

All the effectors in the cardiovascular model change in response to stimulation of sympathetic nerves. The response of a generic effector to sympathetic stimulation includes a pure delay, a monotonic logarithmic static function, and a linear first-order dynamics with a real time constant. The use of logarithmic relationships is justified because the effector static response exhibits a rapid slope at low sympathetic frequencies and then progressively flattens. Of course, the static function is monotonically increasing when the systemic vascular resistance and the end-systolic elastances (left and right ventricles) are considered but monotonically decreasing when heart period is considered. Finally, the first-order dynamics simulates the time required for the effectors to progressively complete its action.

${\sigma _\theta }(t) = \left\{ \begin{gathered} {G_\theta } \cdot \ln [{f_{es}}(t - {D_\theta }) - {f_{es, \min }} + 1], \;if\;{f_{es}} \geqslant {f_{es, \min }} \\ 0, \quad \quad \quad \quad \quad \quad \quad \quad \;\, \, \, \, \, \quad \;\;\, \, if\;{f_{es}} \lt {f_{es, \min }} \\ \end{gathered} \right.$

(23)

$\frac{{d\Delta \theta (t)}}{{dt}} = \frac{1}{{{\tau _\theta }}} \cdot [ - \Delta \theta (t) + {\sigma _\theta }(t)]$

(24)

$\theta (t) = \Delta \theta (t) + {\theta _0}$

(25)

where $\theta$ denotes the generic controlled parameter ( ${R_{svr}}$ , ${E_{\max, lv}}$ , ${E_{\max, rv}}$ , or $T$ ), $\Delta \theta$ is the parameter change caused by sympathetic stimulation, ${\sigma _\theta }$ is the output of the static characteristic, ${\tau _\theta }$ and ${D_\theta }$ are the time constant and pure latency, ${f_{es, \min }}$ is a threshold for sympathetic stimulation, ${G_\theta }$ is a constant gain factor, and ${\theta _0}$ denotes generic controlled parameter in the absence of cardiac innervation.

Heart period is affected not only by the sympathetic but also by the vagal activity. However, the static relationship is assumed to be linear, i.e., heart period increases proportionally to vagal stimulation. The final heart period level is then computed by summing the positive changes induced by the vagal stimulation, the negative changes induced by the sympathetic stimulation, and a constant level representing heart period in the absence of cardiac innervation. Hence, the two mechanisms interact linearly in the heart period control. Of course, this is a crude simplification of the real sympathetic-parasympathetic interaction. However, it provides an acceptable approximation of heart rate changes under a variety of physiological conditions.

${\sigma _{{T_s}}}(t) = \left\{ \begin{gathered} {G_{{T_s}}} \cdot \ln [{f_{es}}(t - {D_{{T_s}}}) - {f_{es, \min }} + 1], \;if\;{f_{es}} \geqslant {f_{es, \min }} \\ 0, \quad \quad \quad \quad \quad \quad \quad \quad \, \;\;\, \, \, \, \;\;\;\;\;\, \, \, if\;{f_{es}} \lt {f_{es, \min }} \\ \end{gathered} \right.$

(26)

$\frac{{d\Delta {T_s}(t)}}{{dt}} = \frac{1}{{{\tau _{{T_s}}}}} \cdot [ - \Delta {T_s}(t) + {\sigma _{{T_s}}}(t)]$

(27)

${\sigma _{{T_v}}}(t) = {G_{{T_v}}} \cdot {f_{ev}}(t - {D_{{T_v}}})$

(28)

$\frac{{d\Delta {T_v}(t)}}{{dt}} = \frac{1}{{{\tau _{{T_v}}}}} \cdot [ - \Delta {T_v}(t) + {\sigma _{{T_v}}}(t)]$

(29)

$T = {T_0} + \Delta {T_s} + \Delta {T_v}$

(30)

where $\Delta {T_s}$ and $\Delta {T_v}$ denote sympathetic and vagal stimulation, ${\sigma _{{T_s}}}$ and ${\sigma _{{T_v}}}$ are the sympathetic and vagal activities, ${\tau _{{T_s}}}$ and ${D_{{T_s}}}$ are the time constant and pure latency for sympathetic, ${\tau _{{T_v}}}$ and ${D_{{T_v}}}$ are the time constant and pure latency for vagus, ${G_{{T_s}}}$ and ${G_{{T_v}}}$ are constant gain factors, and ${T_0}$ denotes heart period in the absence of cardiac innervation.

2.3. LVAD model

A continuous flow LVAD model was connected to the left ventricular and aortic of the CLP model in . It added one state variable to the original model, namely the pressure difference between the pump outlet and pump inlet $H$ as a function of pump flow rate and speed ^[4]:

$H = {\beta _a} \cdot {Q_{irp}}^2(t) + {\beta _b} \cdot {Q_{irp}}(t) \cdot \omega (t) + {\beta _c} \cdot {\omega ^2}(t)$

(31)

where ${Q_{irp}}$ is the pump flow rate (mL/s), $\omega$ is pump speed (revolutions per minute, rpm), and ${\beta _a}$ , ${\beta _b}$ , ${\beta _c}$ are coefficients associated with a rotary pump. In this study, ${\beta _a}$ = −1.8×10⁻³ mmHg·s²·mL⁻², ${\beta _b}$ = −1.2×10⁻⁵ mmHg·s·mL⁻¹·rpm⁻¹ and ${\beta _c}$ = 7.3×10⁻⁶ mmHg·rpm⁻² are derived from a centrifugal pump HeartMate III. Also, corresponding parameters for the resistance and inertance of the outflow cannulae are ${R_i}$ = 0.0677 mmHg·s·mL⁻¹, and ${L_i}$ = 0.0127 mmHg·s²·mL⁻¹. HeartMate III is directly connected to the left ventricular apex and it does not possess inflow cannulae. Therefore, in this study, the parameters for the resistance and inertance of the inflow cannulae are set to 0.

A nonlinear time-varying resistor ${R_k}$ is inserted between the inflow cannulae and the left ventricle to simulate the potential suction event ^[35].

${R_k}(t) = \left\{ \begin{gathered} 0\;, \quad \quad \quad \quad \;\, \, \;\;\;if\;{P_{lv}}(t) \gt \bar P \\ \alpha ({P_{lv}}(t) - \bar P)\;, \quad if\;{P_{lv}}(t) \leqslant \bar P \\ \end{gathered} \right.$

(32)

where $\alpha$ is weight parameter and its value is taken as −3.5 s/mL, $\bar P$ is threshold pressure, and its value is taken as 1 mmHg.

2.4. Simulation protocol

In this study, most parameters under normal conditions in the CLP model and the baroreflex model have values taken from the published literature. The parameters in the CLP model were assigned the values shown in Table A.1 based on the parameter selection in ^[7,14,31]. The parameters in the baroreflex model were set as shown in Table A.2 according to Ref. ^[24]. Table A.1 and Table A.2 were listed in Appendix A.

Compared with normal conditions, heart failure can cause a series of changes, including changes in the cardiovascular system and the baroreflex system. Dilated cardiomyopathy (DCM) is the most common type of nonischemic cardiomyopathy and a major cause of congestive heart failure. The hemodynamic characteristics of patients with DCM usually show some changes, including decreased ventricular end-systolic elastance, increased heart rate, increased systemic vascular resistance, decreased systemic vascular compliance, and increased pulmonary vascular resistance ^[9]. Meanwhile, in the baroreflex system, the baroreflex sensitivity and baroreflex gain for heart failure decreased ^[36,37]. These qualitative and quantitative changes can be mapped into the heart failure model. In the heart failure model, initial changes in left ventricular end-systolic elastance, heart rate, and systemic vascular resistance were simulated by decreasing ${E_{\max, lv0}}$ , ${T_0}$ , and increasing ${R_{svr0}}$ . Besides, increases in ${C_{sar}}$ , ${R_{pvr}}$ , and ${R_{pv}}$ represented increased systemic arterial compliance and pulmonary circulation resistance. A decrease in ${1 / {{k_a}}}$ indicated a decrease in baroreflex sensitivity. The gain constants ( ${G_{{R_{svr}}}}$ , ${G_{{E_{\max, \, lv}}}}$ , ${G_{{E_{\max, \, rv}}}}$ , ${G_{{T_s}}}$ , and ${G_{{T_v}}}$ ) were reduced to simulate the change in the baroreflex gain. The parameters were adjusted to simulate the pathological condition of LV failure, and it was set as shown in Table A.3. Table A.3 was listed in Appendix A.

In this study, a series of hemodynamic responses were conducted under normal condition (no LV failure or LVAD support), LV failure condition (no LVAD support) and LV failure condition with different levels of LVAD support (2200 rpm, 3700 rpm, 4300 rpm). In addition, some variables were compared under different conditions, including pressure, flow rate, volume, heart rate, etc. The area of ventricular pressure-volume (P-V) loop can usually be used to measure ventricular stroke work ^[21]. Consequently, the ventricular P-V loop can provide information about the unloading of the ventricle. The P-V loops under different conditions were also compared in this study.

With different levels of pump support, the LVAD operating condition can be evaluated by changes in pump flow rate and physiological variables. Firstly, the HeartMate III pump speed increased linearly from 0 rpm to 5500 rpm within 60 s. The pump speed range, without pump regurgitation and ventricular suction, could be preliminarily determined by the information of the instantaneous pump flow rate during the linear increase in pump speed. The preliminary speed range was obtained by an instantaneous increase in speed, and it may have errors compared to the stable result with constant speed support. Thus, it is necessary to analyze the hemodynamic variables under constant speed support. Then, the constant pump speeds with 100 rpm interval were selected in the preliminary pump speed range, and the stable values of hemodynamic variables corresponding to constant speeds were calculated.

According to the numerical model described above, a program was developed in MATLAB 2010b (The MathWorks Inc., USA) to perform different simulation experiments. The ode45 solver was used to solve time-varying differential equations. In the system of differential equations, the variable values obtained in the last ode45 calculation was used as the initial value in the next ode45 calculation, and the calculation was repeated until the set time was reached. The time step was set to 0.001 second. All simulation times were set to 60 seconds. The system reaches a stable periodic solution after about 15 seconds. The stable result from the 58th to 60th seconds was chosen in all simulated cases for analysis.

3. Results

3.1. Hemodynamic response under normal condition

Figure 3 depicts the changes in pressure, flow, and volume under normal condition.

Figure 3. Pressures in left ventricle (

${P_{lv}}$ ), left atrium (

${P_{la}}$ ), aorta (

${P_{ao}}$ ), systemic vein (

${P_{sv}}$ ), right ventricle (

${P_{rv}}$ ), right atrium (

${P_{ra}}$ ), pulmonary artery (

${P_{pa}}$ ), and pulmonary vein (

${P_{pv}}$ ), total flow rates through aorta (

${Q_{ao}}$ ) and pulmonary artery (

${Q_{pa}}$ ), and volumes in left ventricle (

${V_{lv}}$ ), left atrium (

${V_{la}}$ ), right ventricle (

${V_{rv}}$ ) and right atrium (

${V_{ra}}$ ) under normal condition.

DownLoad: Full-Size Img PowerPoint

The heart rate was about 68 bpm (beat per minute, bpm).

The left ventricular systolic pressure was 120 mmHg, the aortic pressure changed between 73 and 118 mmHg, and the systemic venous pressure was about 7 mmHg. The right ventricular pressure changed between 3 and 27 mmHg, pulmonary artery changed between 10 and 27 mmHg, and pulmonary venous pressure was about 7 mmHg. The pressures of the left and right atria fluctuated slightly, ranging from 4 to 7 mmHg and 3 to 7 mmHg, respectively.

The left or right ventricle had a stroke volume of about 81 mL and a cardiac output of about 5.5 L/min.

The peak flow rate of the aorta or pulmonary artery was close to 700 mL/s.

These data and waveforms are all consistent with hemodynamic information in the Textbook of Medical Physiology ^[38,39].

3.2. Hemodynamic response under LV failure condition

The hemodynamic response under the LV failure condition is given in Figure 4.

Figure 4. Pressures in left ventricle (

${P_{lv}}$ ), left atrium (

${P_{la}}$ ), aorta (

${P_{ao}}$ ), systemic vein (

${P_{sv}}$ ), right ventricle (

${P_{rv}}$ ), right atrium (

${P_{ra}}$ ), pulmonary artery (

${P_{pa}}$ ), and pulmonary vein (

${P_{pv}}$ ), total flow rates through aorta (

${Q_{ao}}$ ) and pulmonary artery (

${Q_{pa}}$ ), and volumes in left ventricle (

${V_{lv}}$ ), left atrium (

${V_{la}}$ ), right ventricle (

${V_{rv}}$ ) and right atrium (

${V_{ra}}$ ) under LV failure condition.

DownLoad: Full-Size Img PowerPoint

The heart rate was adjusted to 91 bpm. LV failure caused an increase in heart rate.

The aortic pressure range was reduced from the normal value of 73–118 mmHg to the current 53−77 mmHg, the value of aortic pulse pressure (the difference between the maximum systolic pressure and minimum diastolic pressure) was reduced from the normal value of 45 mmHg to the current value of 24 mmHg, and the left ventricular systolic pressure was reduced from the normal value of 120 mmHg to the current value of 78 mmHg. The pressures of the left atrium and pulmonary vein were elevated to about 10 mmHg, almost increased 43% from the normal value of about 7 mmHg. When the LV failure develops further and the pressure of the pulmonary vein is elevated to above 18 mmHg, pulmonary edema will manifest itself and serious consequence of pulmonary failure will follow ^[39].

The peak flow rate of the aorta or pulmonary artery was reduced about 60% of the normal healthy values.

It is observed that there are prominent changes to all four-chamber volumes in the heart. The current LV volume varied between 179 and 211 mL, greatly elevated from the normal range of 57−138 mL. The stroke volume decreased from the normal value of 81 mL to the current value of 32 mL. The ejection fraction (the ratio of stroke volume to end-diastolic volume) is usually not less than 50%, but the current ejection fraction of the left ventricle was about 15%. The current value of cardiac output in the left ventricle (LVCO) was 2.9 L/min, which was reduced 47% from the normal value of 5.5 L/min. The left atrial volume was also increased from the normal range of 28−54 mL to the current range of 43−65 mL. Accompanying the dilation in the left heart, there was a great volume reduction in the right ventricle and atrium. The right ventricle volume was reduced from the normal range of 39−120 mL to the current range of 26−58 mL. The right atrial volume was also reduced from the normal range of 22−47 mL to the current range of 12−26 mL.

3.3. Hemodynamic response under LVAD support

Figure 5 shows the changes in pressure, flow rate, and volume under HeartMate III support at low (2200 rpm), medium (3700 rpm), and high (4300 rpm) speeds.

Figure 5. Pressures in left ventricle (

${P_{lv}}$ ), left atrium (

${P_{la}}$ ), aorta (

${P_{ao}}$ ), systemic vein (

${P_{sv}}$ ), right ventricle (

${P_{rv}}$ ), right atrium (

${P_{ra}}$ ), pulmonary artery (

${P_{pa}}$ ), and pulmonary vein (

${P_{pv}}$ ), total flow rates through aorta (

${Q_{ao}}$ ) and pulmonary artery (

${Q_{pa}}$ ), and volumes in left ventricle (

${V_{lv}}$ ), left atrium (

${V_{la}}$ ), right ventricle (

${V_{rv}}$ ) and right atrium (

${V_{ra}}$ ) under HeartMate III support at 2200 rpm, 3700 rpm or 4300 rpm.

DownLoad: Full-Size Img PowerPoint

Under LV failure condition with pump support at 4300 rpm, left ventricular pressure approached 0 mmHg and the pressure and flow rate in aorta almost lost its pulsation. Further increased in the support level could cause a collapse in the left ventricle.

At 2200 rpm pump operating speed, the flow rates in the rotary pump and aorta had the intermittent regurgitation (the negative flow rate in numerical model represents regurgitation) and the aortic peak flow rate was larger than that of rotary pump over a heart period. It suggested that the total aortic flow rate was provided by both the left ventricle and the rotary pump.

At 3700 rpm or 4300 rpm, the total aortic flow rate and the pump flow rate completely coincided and the aortic pressure was greater than the left ventricular pressure over a heart period. It indicated that the current aortic valve was closed and the total aortic flow rate was completely provided by the rotary pump.

Figure 6 gives the left and right ventricular P-V loops under normal condition, LV failure condition, and different levels of pump support.

Figure 6. The (a) left ventricular P-V loop and (b) right ventricular P-V loop under normal condition (Normal), LV failure condition (LVF) and LVAD support at 2200 rpm, 3700 rpm or 4300 rpm.

DownLoad: Full-Size Img PowerPoint

In contrast to the P-V loop in normal condition, the left ventricular P-V loop shifted to the right due to left ventricular end-diastolic volume dilation in the LV failure condition (Figure 6a). At the same time, the area of the left ventricular P-V loop decreased because of the decreased left ventricular end-systolic pressure and stroke volume. In addition, the right ventricular P-V loop shifted to left and its area decreased in LV failure condition relative to that in normal condition (Figure 6b). Compared with the P-V loop in LV failure condition, the left ventricular P-V loop shifted to the right and its area (stroke work) increased due to increased left ventricular end-diastolic volume and end-systolic pressure at 2200 rpm operating speed. Moreover, at 3700 rpm or 4300 rpm pump operating speed, the left ventricular P-V loop shifted to left and its area decreased, even near to 0 (Figure 6a). With the increase of pump support level, the right ventricular P-V loop continuously shifted to the right and its area (stroke work) was also increasing relative to that in LV failure condition (Figure 6b). And the area of the right ventricular P-V loop at 3700 rpm was relatively closer to that under normal condition.

Ranges of pressures and volumes in left and right heart, arterial and venous pressures, cardiac output, ventricular end-systolic elastances, systemic vascular resistance, and heart rate are shown in Table 1 under different conditions. Under LVAD support, the cardiac output in the left heart was provided by both left ventricle and rotary pump.

Table 1. Simulation values of hemodynamic variables under normal condition, LV failure condition and LVAD support at 2200 rpm, 3700 rpm or 4300 rpm.

	Normal	LVF	2200 rpm	3700 rpm	4300 rpm
${P_{lv}}$ ( ${\rm{mmHg}}$ )	4−120	8−78	8−79	4−26	1−3
${P_{rv}}$ ( ${\rm{mmHg}}$ )	3−27	2−20	2−21	3−29	5−37
${P_{la}}$ ( ${\rm{mmHg}}$ )	4−8	6−11	7−13	5−8	2−4
${P_{ra}}$ ( ${\rm{mmHg}}$ )	3−7	1−3	2−4	3−7	5−11
${P_{ao}}$ ( ${\rm{mmHg}}$ )	73−118	53−77	41−90	95−108	119−120
${P_{pa}}$ ( ${\rm{mmHg}}$ )	10−26	13−20	13−20	12−28	10−36
${\rm{MAP}}$ ( ${\rm{mmHg}}$ )	97	64.2	63.9	100	119.5
${\rm{MPAP}}$ ( ${\rm{mmHg}}$ )	18	17	17	20	23
${P_{sv}}$ ( ${\rm{mmHg}}$ )	7	3	4	7	11
${P_{pv}}$ ( ${\rm{mmHg}}$ )	7	10	10	8	5
${V_{lv}}$ ( ${\rm{mL}}$ )	57−138	179−211	190−270	94−153	28−58
${V_{rv}}$ ( ${\rm{mL}}$ )	39−120	26−58	28−66	40−123	56−193
${V_{la}}$ ( ${\rm{mL}}$ )	28−54	43−65	50−73	32−58	16−30
${V_{ra}}$ ( ${\rm{mL}}$ )	22−47	12−26	13−27	23−50	34-75
${\rm{LVCO}}$ ( ${\rm{L/min}}$ )	5.5	2.9	3.1	6.0	8.7
${\rm{RVCO}}$ ( ${\rm{L/min}}$ )	5.5	2.9	3.1	6.0	8.7
${E_{\max, lv}}$ ( ${\rm{mmHg}}/{\rm{mL}}$ )	2.19	0.44	0.38	0.25	0.14
${E_{\max, \, rv}}$ ( ${\rm{mmHg}}/{\rm{mL}}$ )	0.84	1.19	1.09	0.88	0.71
${R_{svr}}$ ( ${\rm{mmHg}} \cdot {\rm{s}} \cdot {\rm{m}}{{\rm{L}}^{{\rm{ - }}1}}$ )	0.95	1.23	1.12	0.89	0.71
${\rm{HR}}$ ( ${\rm{bpm}}$ )	68	91	83	72	64
* MAP and MPAP denote mean aortic pressure and mean pulmonary artery pressure, LVCO and RVCO denote cardiac output in left ventricle and right ventricle

| Show Table

DownLoad: CSV

Compared with the LV failure condition, the ranges of pressures and volumes in the left or right heart increased at 2200 rpm. But the ranges of pressures and volumes in the left heart decreased and them in the right heart increased, at 3700 rpm or 4300 rpm. Furthermore, both arterial pressure and cardiac output increased with the LVAD support. At 3700 rpm, the mean aortic pressure (MAP) was 100 mmHg and cardiac output in the left ventricle (LVCO) or right ventricle (RVCO) was 6 L/min, which was closer to that in normal condition. With the increase of the LVAD support level, the values of left or right ventricular end-systolic elastance, systemic vascular resistance, and heart rate decreased. At 3700 rpm LVAD operating speed, the heart rate was restored from 91 bpm in LV failure condition to the current 72 bpm.

3.4. Effect of pump speed on hemodynamic variables

Figure 7 shows the instantaneous pump flow rate during a linear increase in pump speed from 0 rpm to 5500 rpm. The pump flow rate had instantaneous or persistent regurgitation when the pump speed was between 0 and 3594 rpm. And the pump regurgitation disappeared when pump speed was higher than 3594 rpm. However, the amplitude of the pump flow rate signal suddenly increased when the speed reached 4712 rpm. It indicated that the left ventricle had an onset of suction. Therefore, the preliminary pump speed range should be between 3594 rpm and 4712 rpm.

Figure 7. Pump flow rate signal (Q_irp) as a function of linearly increasing pump speed.

DownLoad: Full-Size Img PowerPoint

For convenience of calculation, the constant speeds were varied from 3600 rpm to 4700 rpm with 100 rpm interval, and the corresponding changes in hemodynamic variables are shown in Figure 8.

Figure 8. The hemodynamic variables under LV failure with LVAD support as a function of constant speed varied from 3600 rpm to 4700 rpm with 100 rpm interval.

$P$ and

$V$ denote pressure and volume, LVCO and RVCO denote cardiac output in left ventricle and right ventricle, MAP and MPAP denote mean aortic pressure and mean pulmonary artery pressure, PP denotes aortic pulse pressure,

${Q_{irp, min}}$ denotes minimum pump flow rate over a heart period, subscripts

$lv$ ,

$la$ ,

$rv$ and

$ra$ denote left ventricle, left atrium, right ventricle and right atrium,

$ao$ and

$pa$ denote aorta and pulmonary artery,

$sys$ and

$dia$ denote systolic and diastolic,

$ed$ and

$es$ denote end-diastolic and end-systolic.

DownLoad: Full-Size Img PowerPoint

With the increase of pump speed, the pressures and volumes in the left heart decreased gradually. On the contrary, the pressures and volumes in the right heart gradually increased. Meanwhile, aortic pressure and pulmonary artery pressure were increasing. Moreover, the pulsatility of aortic pressure gradually weakened, and that of pulmonary artery pressure gradually enhanced. The increase of speed also enhanced the cardiac output and the minimum pump flow rate over a heart period. There was no negative value of the minimum pump flow rate within 3600−4700 rpm, which indicated that pump regurgitation disappeared when the speed was higher than 3600 rpm. However, when the pump operating speed is higher than 4400 rpm, the minimum pump flow rate decreased, the aortic pulse pressure increased, the left ventricular pressure approached 0 mmHg, and the left ventricular volume was less than 50 mL, which implied that the left ventricular suction occurred.

4. Discussion

In this study, a CLP model for describing the circulatory loop was developed to simulate hemodynamic responses. The CLP model was coupled with the baroreflex model and LVAD model, which provided a platform for the simulation of hemodynamics and the evaluation of LVAD operating condition.

LVADs have become the standard treatment for advanced heart failure. Nevertheless, pump regurgitation may occur under pump support at low speeds. Inadequate support may cause mean aortic pressure and cardiac output to be lower than the normal level, which is difficult to meet the patients' normal physiological perfusion. Pump support at high speed allows the left ventricle to unload, but it also increases right ventricular load (Figure 6). However, long-term overloading for the right ventricle may cause right heart failure ^[40].

The HeartMate III operating speed should be in the range of 3600−4400 rpm with the selected parameter settings to avoid pump regurgitation and ventricular suction. Generally, an optimal speed range is set for clinical continuous-flow LVADs, and it can be determined by monitoring cardiac index (cardiac output) and left ventricular size based on hemodynamic data ^[40]. In this study, cardiac output in the left ventricle and left ventricular end-diastolic volume are 5.5 L/min and 138 mL in normal conditions, respectively. In Figure 8, the cardiac output of 5.6 L/min at 3600 rpm is near to the normal cardiac output. And at 3800 rpm, the left ventricular end-systolic volume is 139 mL that is closer to the normal value of 138 mL. Hence, the low-end of the optimal speed range can be set to 3600 rpm and the high-end can be set to 3800 rpm.

In the physiological or pathological cardiovascular system, the baroreflex model objectively describes the changes in ventricular end-systolic elastances, systemic vascular resistance, and heart rate. Different levels of LVAD support also affect ventricular end-systolic elastances, systemic vascular resistance, and heart rate. However, heart rate is usually maintained as a constant in cardiovascular models without baroreflex regulation, which may not be in line with the actual situation.

A pulse pressure of greater than 15 mmHg is considered to be pulsatile flow while a pulse pressure of less than 15 mmHg is considered to be nonpulsatile flow ^[41]. In this study, the aortic pulse pressure is less than 15 mmHg for the operating speed between 3600 rpm and 4400 rpm in Figure 8. The continuous flow LVADs usually operate at a constant speed to provide the patient with adequate blood flow. However, long-term LVAD constant support diminishes pulse pressure and vascular pulsatility, which can lead to late complications such as aortic insufficiency, gastrointestinal bleeding, and hemorrhagic strokes ^[31,42]. Recently, the new generation of commercial products, such as HVAD (Medtronic, Inc., Minneapolis, MN, USA) and HeartMate III (Abbott Inc., Abbott Park, IL, USA), come with an artificial pulse. The artificial pulse is formed by intermittent speed change, and its main goal is to reduce the rate of pump thrombosis. Importantly, the cycle of variable speed is asynchronous with the native heart period. Moreover, the aortic pressure waveform becomes non-physiological ^[43]. Therefore, certain issues should be studied and solved, including modulation and optimization of pulsating speed with heart rate synchronization, hemodynamic analysis at pulsating speed support, and design of the physiological controller. These will be the subject of further research.

The limitations of this study are as follows.

From a hemodynamic viewpoint, the numerical model can reflect clinical observations under LV failure and LVAD support conditions. The numerical model cannot replace the importance and significance of in vivo research, and cannot reproduce all the expected clinical responses (such as neurohumoral response, myocardial remodeling). But it can prove the feasibility of some concepts to provide a valuable initial step for further in vivo research.

In this study, the valve model is assumed that valves can be opened and closed immediately, and it ignores the process of valve closure. Thus, the ideal valve model cannot simulate the valve flow rate during valve closure. The non-ideal diode modeling for the heart valves is a feasible way to simulate the dynamic process of valves ^[44,45,46].

The experimental data of HeartMate III are from earlier versions of the product design ^[4]. Since the pump data was first published, the HeartMate III design has gone through several versions. However, the method in this study is still valid.

In this study, the regulation of venous volume in the baroreflex model is ignored. This is because the systemic venous compliance is configured to be constant. But this does not affect the change in systemic venous pressure.

5. Conclusion

In this study, a numerical model had been developed to simulate the hemodynamics and to evaluate the operating condition of continuous-flow left ventricular assist device. The hemodynamic responses under various conditions (normal, left ventricular failure, and different levels of pump support) had been simulated successfully, and the hemodynamic variables under different conditions had been compared. The simulation results had shown that the developed numerical model provides a feasible simulation platform for reproducing hemodynamics. Also, the centrifugal pump HeartMate III had been taken as an example of continuous-flow LVAD, and its operating condition had been evaluated and discussed. The simulation results had shown that in the simulation system, HeartMate III is safe to operate between 3600 rpm and 4400 rpm, and that pump operating speeds are adjusted between 3600 rpm and 3800 rpm to provide optimal physiological perfusion. What makes this paper regret is that some system parameters are quoted or set as constants, it gives rise to some limitations on the final results. Actually, these parameters are variable because of the interpatient and intrapatient variations. The future work will simulate combining the patient data, and obtain the hemodynamic response for variable physiological or pathological parameters.

Acknowledgments

The study is supported by the Key Technology Research and Development Program of Shandong (CN) (No. 2017GSF221009).

Conflict of interest

The authors declare that there is no conflict of interest.

Appendix A

Table A1. Parameters for cardiovascular system in normal condition ^[7,14,31].

Parameters	Physiological Meaning	Value	Unit
${R_{ao}}$	Aortic Resistance	0.0398	${\rm{mmHg}} \cdot {\rm{s}} \cdot {\rm{m}}{{\rm{L}}^{{\rm{ - }}1}}$
${R_{svr}}(t)$	Systemic Vascular Resistance	Time-varying	${\rm{mmHg}} \cdot {\rm{s}} \cdot {\rm{m}}{{\rm{L}}^{{\rm{ - }}1}}$
${R_{sv}}$	Systemic Venous Resistance	0.023	${\rm{mmHg}} \cdot {\rm{s}} \cdot {\rm{m}}{{\rm{L}}^{{\rm{ - }}1}}$
${R_{pa}}$	Pulmonary Artery Resistance	0.005	${\rm{mmHg}} \cdot {\rm{s}} \cdot {\rm{m}}{{\rm{L}}^{{\rm{ - }}1}}$
${R_{pvr}}$	Pulmonary Vascular Resistance	0.12	${\rm{mmHg}} \cdot {\rm{s}} \cdot {\rm{m}}{{\rm{L}}^{{\rm{ - }}1}}$
${R_{pv}}$	Pulmonary Venous Resistance	0.011	${\rm{mmHg}} \cdot {\rm{s}} \cdot {\rm{m}}{{\rm{L}}^{{\rm{ - }}1}}$
${R_{{\rm{MV}}}}$	Mitral Valve Resistance	0.005	${\rm{mmHg}} \cdot {\rm{s}} \cdot {\rm{m}}{{\rm{L}}^{{\rm{ - }}1}}$
${R_{{\rm{AV}}}}$	Aortic Valve Resistance	0.008	${\rm{mmHg}} \cdot {\rm{s}} \cdot {\rm{m}}{{\rm{L}}^{{\rm{ - }}1}}$
${R_{{\rm{TV}}}}$	Tricuspid Valve Resistance	0.001	${\rm{mmHg}} \cdot {\rm{s}} \cdot {\rm{m}}{{\rm{L}}^{{\rm{ - }}1}}$
${R_{{\rm{PV}}}}$	Pulmonary Valve Resistance	0.003	${\rm{mmHg}} \cdot {\rm{s}} \cdot {\rm{m}}{{\rm{L}}^{{\rm{ - }}1}}$
${C_{lv}}(t)$	Left Ventricular Compliance	Time-varying	${\rm{mL}}/{\rm{mmHg}}$
${C_{la}}(t)$	Left Atrial Compliance	Time-varying	${\rm{mL}}/{\rm{mmHg}}$
${C_{rv}}(t)$	Right Ventricular Compliance	Time-varying	${\rm{mL}}/{\rm{mmHg}}$
${C_{ra}}(t)$	Right Atrial Compliance	Time-varying	${\rm{mL}}/{\rm{mmHg}}$
${C_{ao}}$	Aortic Compliance	0.08	${\rm{mL}}/{\rm{mmHg}}$
${C_{sar}}$	Systemic Arterial Compliance	1.33	${\rm{mL}}/{\rm{mmHg}}$
${C_{sv}}$	Systemic Venous Compliance	20.5	${\rm{mL}}/{\rm{mmHg}}$
${C_{pa}}$	Pulmonary Artery Compliance	0.18	${\rm{mL}}/{\rm{mmHg}}$
${C_{par}}$	Pulmonary Arteriolar Compliance	3	${\rm{mL}}/{\rm{mmHg}}$
${C_{pv}}$	Pulmonary Venous Compliance	20.5	${\rm{mL}}/{\rm{mmHg}}$
${L_{ao}}$	Aortic Inertance	0.0005	${\rm{mmHg}} \cdot {{\rm{s}}^2} \cdot {\rm{m}}{{\rm{L}}^{{\rm{ - }}1}}$
${L_{pa}}$	Pulmonary Artery Inertance	0.0005	${\rm{mmHg}} \cdot {{\rm{s}}^2} \cdot {\rm{m}}{{\rm{L}}^{{\rm{ - }}1}}$
${V_{lv, 0}}$	Left Ventricular Unstressed Volume	5	${\rm{mL}}$
${V_{la, 0}}$	Left Atrial Unstressed Volume	4	${\rm{mL}}$
${V_{rv, 0}}$	Right Ventricular Unstressed Volume	10	${\rm{mL}}$
${V_{ra, 0}}$	Right Atrial Unstressed Volume	4	${\rm{mL}}$
${E_{\max, lv}}(t)$	Left Ventricular End-systolic Elastance	Time-varying	${\rm{mmHg}}/{\rm{mL}}$
${E_{\min, lv}}$	Left Ventricular End-diastolic Elastance	0.04	${\rm{mmHg}}/{\rm{mL}}$
${E_{\max, rv}}(t)$	Right Ventricular End-systolic Elastance	Time-varying	${\rm{mmHg}}/{\rm{mL}}$
${E_{\min, rv}}$	Right Ventricular End-diastolic Elastance	0.04	${\rm{mmHg}}/{\rm{mL}}$
${E_{\max, la}}$	Left Atrial End-systolic Elastance	0.25	${\rm{mmHg}}/{\rm{mL}}$
${E_{min, la}}$	Left Atrial End-diastolic Elastance	0.15	${\rm{mmHg}}/{\rm{mL}}$
${E_{\max, ra}}$	Right Atrial End-systolic Elastance	0.25	${\rm{mmHg}}/{\rm{mL}}$
${E_{min, ra}}$	Right Atrial End-diastolic Elastance	0.15	${\rm{mmHg}}/{\rm{mL}}$

| Show Table

DownLoad: CSV

Table A2. Parameters for baroreflex system in normal condition ^[24].

Parameters
Afferent neural pathway
${\tau _p}{\rm{ = }}$ 2.076 ${\rm{s}}$		${f_{\min }} =$ 2.52 ${\rm{spikes/s}}$				${p_n} =$ 92 ${\rm{mmHg}}$
${\tau _z} =$ 6.37 ${\rm{s}}$		${f_{\max }} =$ 47.78 ${\rm{spikes/s}}$				${k_a} =$ 11.758 ${\rm{mmHg}}$
Efferent neural pathway
${f_{es, \infty }} =$ 2.1 ${\rm{spikes/s}}$	${f_{ev, \infty }} =$ 6.3 ${\rm{spikes/s}}$			${k_{es}} =$ 0.0675 ${\rm{s}}$				${f_{as, 0}} =$ 25 ${\rm{spikes/s}}$
${f_{es, 0}} =$ 16.11 ${\rm{spikes/s}}$	${f_{ev, 0}} =$ 3.2 ${\rm{spikes/s}}$			${k_{ev}} =$ 7.06 ${\rm{spikes/s}}$
Effector
${G_{{R_{svr}}}} =$ 0.36 ${\rm{mmHg}} \cdot {\rm{s}} \cdot {\rm{m}}{{\rm{L}}^{{\rm{ - }}1}} \cdot {\rm{spike}}{{\rm{s}}^{{\rm{ - }}1}}$			${\tau _{{R_{svr}}}} =$ 6 ${\rm{s}}$		${D_{{R_{svr}}}} =$ 2 ${\rm{s}}$		${f_{es, \min }} =$ 2.66 ${\rm{spikes/s}}$
${G_{{E_{\max, \, lv}}}} =$ 0.475 ${\rm{mmHg}} \cdot {\rm{s}} \cdot {\rm{m}}{{\rm{L}}^{{\rm{ - }}1}} \cdot {\rm{spike}}{{\rm{s}}^{{\rm{ - }}1}}$			${\tau _{{E_{\max, \, lv}}}} =$ 8 ${\rm{s}}$		${D_{{E_{\max, \, lv}}}} =$ 2 ${\rm{s}}$		${R_{svr0}} =$ 0.62 ${\rm{mmHg}} \cdot {\rm{s}} \cdot {\rm{m}}{{\rm{L}}^{{\rm{ - }}1}}$
${G_{{E_{\max, \, rv}}}} =$ 0.282 ${\rm{mmHg}} \cdot {\rm{s}} \cdot {\rm{m}}{{\rm{L}}^{{\rm{ - }}1}} \cdot {\rm{spike}}{{\rm{s}}^{{\rm{ - }}1}}$			${\tau _{{E_{\max, \, rv}}}} =$ 8 ${\rm{s}}$		${D_{{E_{\max, \, rv}}}} =$ 2 ${\rm{s}}$		${E_{\max, lv0}} =$ 1.8 ${\rm{mmHg}}/{\rm{mL}}$
${G_{{T_s}}} =$ -0.13 ${{\rm{s}}^2} \cdot {\rm{spike}}{{\rm{s}}^{{\rm{ - }}1}}$			${\tau _{{T_s}}} =$ 2 ${\rm{s}}$		${D_{{T_s}}} =$ 2 ${\rm{s}}$		${E_{\max, rv0}} =$ 0.6 ${\rm{mmHg}}/{\rm{mL}}$
${G_{{T_v}}} =$ 0.09 ${{\rm{s}}^2} \cdot {\rm{spike}}{{\rm{s}}^{{\rm{ - }}1}}$			${\tau _{{T_v}}} =$ 1.5 ${\rm{s}}$		${D_{{T_v}}} =$ 0.2 ${\rm{s}}$		${T_0} =$ 0.58 ${\rm{s}}$

| Show Table

DownLoad: CSV

Table A3. Partially normal parameters are adjusted to simulate LV failure.

Parameters	Normal	LV failure	Unit
${C_{sar}}$	1.33	1.02	${\rm{mL}}/{\rm{mmHg}}$
${R_{pvr}}$	0.12	0.14	${\rm{mmHg}} \cdot {\rm{s}} \cdot {\rm{m}}{{\rm{L}}^{{\rm{ - }}1}}$
${R_{pv}}$	0.011	0.027	${\rm{mmHg}} \cdot {\rm{s}} \cdot {\rm{m}}{{\rm{L}}^{{\rm{ - }}1}}$
${E_{\max, lv0}}$	1.8	0.05	${\rm{mmHg}}/{\rm{mL}}$
${R_{svr0}}$	0.62	0.71	${\rm{mmHg}} \cdot {\rm{s}} \cdot {\rm{m}}{{\rm{L}}^{{\rm{ - }}1}}$
${T_0}$	0.58	0.52	${\rm{s}}$
${k_a}$	11.758	14	${\rm{mmHg}}$
${G_{{R_{svr}}}}$	0.36	0.3	${\rm{mmHg}} \cdot {\rm{s}} \cdot {\rm{m}}{{\rm{L}}^{{\rm{ - }}1}} \cdot {\rm{spike}}{{\rm{s}}^{{\rm{ - }}1}}$
${G_{{E_{\max, \, lv}}}}$	0.475	0.15	${\rm{mmHg}} \cdot {\rm{s}} \cdot {\rm{m}}{{\rm{L}}^{{\rm{ - }}1}} \cdot {\rm{spike}}{{\rm{s}}^{{\rm{ - }}1}}$
${G_{{E_{\max, \, rv}}}}$	0.282	0.245	${\rm{mmHg}} \cdot {\rm{s}} \cdot {\rm{m}}{{\rm{L}}^{{\rm{ - }}1}} \cdot {\rm{spike}}{{\rm{s}}^{{\rm{ - }}1}}$
${G_{{T_s}}}$	-0.13	-0.07	${{\rm{s}}^2} \cdot {\rm{spike}}{{\rm{s}}^{{\rm{ - }}1}}$
${G_{{T_v}}}$	0.09	0.07	${{\rm{s}}^2} \cdot {\rm{spike}}{{\rm{s}}^{{\rm{ - }}1}}$

| Show Table

DownLoad: CSV

References

[1]	B. Han, W. M. Compton, C. M. Jones, E. B. Einstein, N. D. Volkow, Methamphetamine use, methamphetamine use disorder, and associated overdose deaths among US adults, JAMA psychiatry, 78 (2021), 1329–1342. https://doi.org/10.1001/jamapsychiatry.2021.2588 doi: 10.1001/jamapsychiatry.2021.2588
[2]	M. Shukla, B. Vincent, Methamphetamine abuse disturbs the dopaminergic system to impair hippocampal-based learning and memory: An overview of animal and human investigations, Neurosci. Biobehav. Rev., 131 (2021), 541–559. https://doi.org/10.1016/j.neubiorev.2021.09.016 doi: 10.1016/j.neubiorev.2021.09.016
[3]	V. Manja, A. Nrusimha, Y. Gao, A. Sheikh, M. McGovern, P. A. Heidenreich, et al., Methamphetamine-associated heart failure: A systematic review of observational studies, Heart, 109 (2023), 168–177. https://doi.org/10.1136/heartjnl-2022-321610 doi: 10.1136/heartjnl-2022-321610
[4]	C. W. Li, S. W. W. Ku, P. Y. Huang, L. Y. Chen, H. T. Wei, C. Strong, et al., Factors associated with methamphetamine dependency among men who have sex with men engaging in chemsex: Findings from the COMeT study in Taiwan, Int. J. Drug Policy, 93 (2021), 103119. https://doi.org/10.1016/j.drugpo.2021.103119 doi: 10.1016/j.drugpo.2021.103119
[5]	S. P. Xu, K. Zhang, T. Y. Luo, Development of the risk of relapse assessment scale for methamphetamine abusers in China, Drug Alcohol Depend., 227 (2021), 108992. https://doi.org/10.1016/j.drugalcdep.2021.108992 doi: 10.1016/j.drugalcdep.2021.108992
[6]	H. Yuan, X. Yu, X. Li, S. Qin, G. Liang, T. Bai, et al., Research on resting spontaneous brain activity and functional connectivity of acupuncture at uterine acupoints, Digital Chin. Med., 5 (2022), 59–67. https://doi.org/10.1016/j.dcmed.2022.03.006 doi: 10.1016/j.dcmed.2022.03.006
[7]	X. Li, B. Wei, T. Li, N. Zhang, MwoA auxiliary diagnosis via RSN-based 3D deep multiple instance learning with spatial attention mechanism, in 2020 11th International Conference on Awareness Science and Technology (iCAST), IEEE, (2020), 1–6. https://doi.org/10.1109/iCAST51195.2020.9319486
[8]	A. P. Daiwile, S. Jayanthi, J. L. Cadet, Sex differences in methamphetamine use disorder perused from pre-clinical and clinical studies: Potential therapeutic impacts, Neurosci. Biobehav. Rev., 137 (2022), 104674. https://doi.org/10.1016/j.neubiorev.2022.104674 doi: 10.1016/j.neubiorev.2022.104674
[9]	H. Mizoguchi, K. Yamada, Methamphetamine use causes cognitive impairment and altered decision-making, Neurochem. Int., 124 (2019), 106–113. https://doi.org/10.1016/j.neuint.2018.12.019 doi: 10.1016/j.neuint.2018.12.019
[10]	S. Sabrini, G. Y. Wang, J. C. Lin, J. K. Ian, L. E. Curley, Methamphetamine use and cognitive function: A systematic review of neuroimaging research, Drug Alcohol Depend., 194 (2019), 75–87. https://doi.org/10.1016/j.drugalcdep.2018.08.041 doi: 10.1016/j.drugalcdep.2018.08.041
[11]	S. J. Nieto, L. A. Ray, Applying the addictions neuroclinical assessment to derive neurofunctional domains in individuals who use methamphetamine, Behav. Brain Res., 427 (2022), 113876. https://doi.org/10.1016/j.bbr.2022.113876 doi: 10.1016/j.bbr.2022.113876
[12]	G. X. Liang, X. Li, H. Yuan, M. Sun, S. J. Qin, B. Z. Wei, Abnormal static and dynamic amplitude of low-frequency fluctuations in multiple brain regions of methamphetamine abstainers, Math. Biosci. Eng., 20 (2023), 13318–13333. https://doi.org/10.3934/mbe.2023593 doi: 10.3934/mbe.2023593
[13]	Z. X. Zhang, L. He, S. C. Huang, L. D. Fan, Y. N. Li, P. Li, et al., Alteration of brain structure with long-term abstinence of methamphetamine by voxel-based morphometry, Front. Psychiatry, 9 (2018), 722. https://doi.org/10.3389/fpsyt.2018.00722 doi: 10.3389/fpsyt.2018.00722
[14]	X. T. Li, H. Su, N. Zhong, T. Z. Chen, J. Du, K. Xiao, et al., Aberrant resting-state cerebellar-cerebral functional connectivity in methamphetamine-dependent individuals after six months abstinence, Front. Psychiatry, 11 (2020), 191. https://doi.org/10.3389/fpsyt.2020.00191 doi: 10.3389/fpsyt.2020.00191
[15]	L. Fan, Q. Zhang, S. Liang, H. Li, Z. He, J. Sun, et al., Imaging changes in brain microstructural in long-term abstinent from methamphetamine-dependence (in Chinese), J. Cent. South Univ. (Med. Sci.), 44 (2019), 491–500. https://doi.org/10.11817/j.issn.1672-7347.2019.05.004 doi: 10.11817/j.issn.1672-7347.2019.05.004
[16]	X. Y. Qi, Y. Y. Wang, Y. Z. Lu, Q. Zhao, Y. F. Chen, C. L. Zhou, et al., Enhanced brain network flexibility by physical exercise in female methamphetamine users, Cognit. Neurodyn., 2022 (2022), 1–17. https://doi.org/10.1007/s11571-022-09848-5 doi: 10.1007/s11571-022-09848-5
[17]	C. Y. Jia, Q. F. Long, T. Ernst, Y. Q. Shang, L. D. Chang, T. Adali, Independent component and graph theory analyses reveal normalized brain networks on resting-state functional MRI after working memory training in people with HIV, J. Magn. Reson. Imaging, 57 (2023), 1552–1564. https://doi.org/10.1002/jmri.28439 doi: 10.1002/jmri.28439
[18]	O. Sporns, Graph theory methods: applications in brain networks, Dialogues Clin. Neurosci., 20 (2018), 111–121. https://doi.org/10.31887/DCNS.2018.20.2/osporns doi: 10.31887/DCNS.2018.20.2/osporns
[19]	W. Li, L. Wang, Z. Lyu, J. J. Chen, Y. B. Li, Y. C. Sun, et al., Difference in topological organization of white matter structural connectome between methamphetamine and heroin use disorder, Behav. Brain Res., 422 (2022), 113752. https://doi.org/10.1016/j.bbr.2022.113752 doi: 10.1016/j.bbr.2022.113752
[20]	M. Siyah Mansoory, M. A. Oghabian, A. H. Jafari, A. Shahbabaie, Analysis of resting-state fMRI topological graph theory properties in methamphetamine drug users applying box-counting fractal dimension, Basic Clin. Neurosci., 8 (2017), 371–385. https://doi.org/10.18869/nirp.bcn.8.5.371 doi: 10.18869/nirp.bcn.8.5.371
[21]	Y. A. Zhou, Y. Hu, Q. J. Wang, Z. Yang, J. G. Li, Y. J. Ma, et al., Association between white matter microstructure and cognitive function in patients with methamphetamine use disorder, Hum. Brain Mapp., 44 (2023), 304–314. https://doi.org/10.1002/hbm.26020 doi: 10.1002/hbm.26020
[22]	M. S. Mansoory, A. Allahverdy, M. Behboudi, M. Khodamoradi, Local efficiency analysis of restingstate functional brain network in methamphetamine users, Behav. Brain Res., 434 (2022), 114022. https://doi.org/10.1016/j.bbr.2022.114022 doi: 10.1016/j.bbr.2022.114022
[23]	F. Miraglia, F. Vecchio, C. Pappalettera, L. Nucci, M. Cotelli, E. Judica, et al., Brain connectivity and graph theory analysis in Alzheimer's and Parkinson's Disease: The contribution of electrophysiological techniques, Brain Sci., 12 (2022), 402. https://doi.org/10.3390/brainsci12030402 doi: 10.3390/brainsci12030402
[24]	C. G. Yan, X. D. Wang, X. N. Zuo, Y. F. Zang, DPABI: Data processing & analysis for (resting-state) brain imaging, Neuroinformatics, 14 (2016), 339–351. https://doi.org/10.1007/s12021-016-9299-4 doi: 10.1007/s12021-016-9299-4
[25]	C. G. Yan, R. C. Craddock, Y. He, M. P. Milham, Addressing head motion dependencies for small-world topologies in functional connectomics, Front. Hum. Neurosci., 7 (2013), 910. https://doi.org/10.3389/fnhum.2013.00910 doi: 10.3389/fnhum.2013.00910
[26]	C. Yan, Y. Zang, DPARSF: A MATLAB toolbox for "pipeline" data analysis of resting-state fMRI, Front. Syst. Neurosci., 4 (2010), 13. https://doi.org/10.3389/fnsys.2010.00013 doi: 10.3389/fnsys.2010.00013
[27]	N. U. F. Dosenbach, B. Nardos, A. L. Cohen, D. A. Fair, J. D. Power, J. A. Church, et al., Prediction of individual brain maturity using fMRI, Science, 329 (2010), 1358–1361. https://doi.org/10.1126/science.1194144 doi: 10.1126/science.1194144
[28]	M. Koutrouli, E. Karatzas, D. Paez-Espino, G. A. Pavlopoulos, A Guide to conquer the biological network Era using graph theory, Front. Bioeng. Biotechnol., 8 (2020), 34. https://doi.org/10.3389/fbioe.2020.00034 doi: 10.3389/fbioe.2020.00034
[29]	B. T. T. Yeo, F. M. Krienen, J. Sepulcre, M. R. Sabuncu, D. Lashkari, M. Hollinshead, et al., The organization of the human cerebral cortex estimated by intrinsic functional connectivity, J. Neurophysiol., 106 (2011), 1125–1165. https://doi.org/10.1152/jn.00338.2011 doi: 10.1152/jn.00338.2011
[30]	G. Li, Y. D. Luo, Z. R. Zhang, Y. T. Xu, W. D. Jiao, Y. H. Jiang, et al., Effects of mental fatigue on < i > Small-World < /i > brain functional network organization, Neural Plast., 2019 (2019), 1716074. https://doi.org/10.1155/2019/1716074 doi: 10.1155/2019/1716074
[31]	L. J. Nestor, D. G. Ghahremani, E. D. London, Reduced neural functional connectivity during working memory performance in methamphetamine use disorder, Drug Alcohol Depend., 243 (2023), 109764. https://doi.org/10.1016/j.drugalcdep.2023.109764 doi: 10.1016/j.drugalcdep.2023.109764
[32]	M. Ahmadlou, K. Ahmadi, M. Rezazade, E. Azad-Marzabadi, Global organization of functional brain connectivity in methamphetamine abusers, Clin. Neurophysiol., 124 (2013), 1122–1131. https://doi.org/10.1016/j.clinph.2012.12.003 doi: 10.1016/j.clinph.2012.12.003
[33]	S. Arvin, A. N. Glud, K. Yonehara, Short- and long-range connections differentially modulate the dynamics and state of small-world networks, Front. Comput. Neurosci., 15 (2022), 783474. https://doi.org/10.3389/fncom.2021.783474 doi: 10.3389/fncom.2021.783474
[34]	Y. Liu, Q. Li, T. Y. Zhang, L. Wang, Y. R. Wang, J. J. Chen, et al., Differences in small-world networks between methamphetamine and heroin use disorder patients and their relationship with psychiatric symptoms, Brain Imaging Behav., 16 (2022), 1973–1982. https://doi.org/10.1007/s11682-022-00667-0 doi: 10.1007/s11682-022-00667-0
[35]	H. Khajehpour, B. Makkiabadi, H. Ekhtiari, S. Bakht, A. Noroozi, F. Mohagheghian, Disrupted resting-state brain functional network in methamphetamine abusers: A brain source space study by EEG, PLoS One, 14 (2019), e0226249. https://doi.org/10.1371/journal.pone.0226249 doi: 10.1371/journal.pone.0226249
[36]	F. X. Vollenweider, K. H. Preller, Psychedelic drugs: neurobiology and potential for treatment of psychiatric disorders, Nat. Rev. Neurosci., 21 (2020), 611–624. https://doi.org/10.1038/s41583-020-0367-2 doi: 10.1038/s41583-020-0367-2
[37]	B. Kim, J. Yun, B. Park, Methamphetamine-induced neuronal damage: Neurotoxicity and neuroinflammation, Biomol. Ther., 28 (2020), 381–388. https://doi.org/10.4062/biomolther.2020.044 doi: 10.4062/biomolther.2020.044
[38]	J. Zhang, H. Su, J. Y. Tao, Y. Xie, Y. M. Sun, L. R. Li, et al., Relationship of impulsivity and depression during early methamphetamine withdrawal in Han Chinese population, Addict. Behav., 43 (2015), 7–10. https://doi.org/10.1016/j.addbeh.2014.10.032 doi: 10.1016/j.addbeh.2014.10.032
[39]	A. Scalabrini, B. Vai, S. Poletti, S. Damiani, C. Mucci, C. Colombo, et al., All roads lead to the default-mode network-global source of DMN abnormalities in major depressive disorder, Neuropsychopharmacology, 45 (2020), 2058–2069. https://doi.org/10.1038/s41386-020-0785-x doi: 10.1038/s41386-020-0785-x
[40]	B. L. Foster, S. R. Koslov, L. Aponik-Gremillion, M. E. Monko, B. Y. Hayden, S. R. Heilbronner, A tripartite view of the posterior cingulate cortex, Nat. Rev. Neurosci., 24 (2023), 173–189. https://doi.org/10.1038/s41583-022-00661-x doi: 10.1038/s41583-022-00661-x
[41]	C. Caldinelli, R. Cusack, The fronto-parietal network is not a flexible hub during naturalistic cognition, Hum. Brain Mapp., 43 (2022), 750–759. https://doi.org/10.1002/hbm.25684 doi: 10.1002/hbm.25684
[42]	H. Zheng, Q. Zhou, J. J. Yang, Q. Lu, H. D. Qiu, C. He, et al., Altered functional connectivity of the default mode and frontal control networks in patients with insomnia, CNS Neurosci. Ther., 29 (2023), 2318–2326. https://doi.org/10.1111/cns.14183 doi: 10.1111/cns.14183
[43]	R. De Micco, N. Piramide, F. Di Nardo, M. Siciliano, M. Cirillo, A. Russo, et al., Resting-state network connectivity changes in drug-naive Parkinson's disease patients with probable REM sleep behavior disorder, J. Neural Transm., 130 (2023), 43–51. https://doi.org/10.1007/s00702-022-02565-7 doi: 10.1007/s00702-022-02565-7
[44]	S. Tikoo, F. Cardona, S. Tommasin, C. Giannì, G. Conte, N. Upadhyay, et al., Resting-state functional connectivity in drug-naive pediatric patients with Tourette syndrome and obsessive-compulsive disorder, J. Psychiatr. Res., 129 (2020), 129–140. https://doi.org/10.1016/j.jpsychires.2020.06.021 doi: 10.1016/j.jpsychires.2020.06.021
[45]	G. Dong, E. DeVito, J. Huang, X. Du, Diffusion tensor imaging reveals thalamus and posterior cingulate cortex abnormalities in internet gaming addicts. J. Psychiatr. Res., 46 (2012), 1212–1216. https://doi.org/10.1016/j.jpsychires.2012.05.015 doi: 10.1016/j.jpsychires.2012.05.015
[46]	Y. Katsumi, D. Putcha, R. Eckbo, B. Wong, M. Quimby, S. McGinnis, et al., Anterior dorsal attention network tau drives visual attention deficits in posterior cortical atrophy, Brain, 146 (2023), 295–306. https://doi.org/10.1093/brain/awac245 doi: 10.1093/brain/awac245
[47]	H. Y. Tan, T. Z. Chen, J. Du, R. J. Li, H. F. Jiang, C. L. Deng, et al., Drug-related virtual reality cue reactivity is associated with gamma activity in reward and executive control circuit in methamphetamine use disorders, Arch. Med. Res., 50 (2019), 509–517. https://doi.org/10.1016/j.arcmed.2019.09.003 doi: 10.1016/j.arcmed.2019.09.003
[48]	H. C. Zhao, M. J. Ge, O. Turel, A. Bechara, Q. H. He, Brain modular connectivity interactions can predict proactive inhibition in smokers when facing smoking cues, Addict. Biol., 28 (2023), e13284. https://doi.org/10.1111/adb.13284 doi: 10.1111/adb.13284
[49]	S. Yilmaz, Impaired biological rhythm in men with methamphetamine use disorder: The relationship with sleep quality and depression, J. Subst. Use, 28 (2023), 280–286. https://doi.org/10.1080/14659891.2022.2098847 doi: 10.1080/14659891.2022.2098847
[50]	L. Cerliani, M. Mennes, R. M. Thomas, A. Di Martino, M. Thioux, C. Keysers, Increased functional connectivity between subcortical and cortical resting-state networks in autism spectrum disorder, JAMA Psychiatry, 72 (2015), 767–777. https://doi.org/10.1001/jamapsychiatry.2015.0101 doi: 10.1001/jamapsychiatry.2015.0101
[51]	A. Moaveni, Y. F. Feyzi, S. T. Rahideh, R. Arezoomandan, The relationship between serum brain-derived neurotrophic level and neurocognitive functions in chronic methamphetamine users, Neurosci. Lett., 772 (2022), 136478. https://doi.org/10.1016/j.neulet.2022.136478 doi: 10.1016/j.neulet.2022.136478
[52]	R. Nusslock, G. H. Brody, C. C. Armstrong, A. L. Carroll, L. H. Sweet, T. Y. Yu, et al., Higher peripheral inflammatory signaling associated with lower resting-state functional brain connectivity in emotion regulation and central executive networks, Biol. Psychiatry, 86 (2019), 153–162. https://doi.org/10.1016/j.biopsych.2019.03.968 doi: 10.1016/j.biopsych.2019.03.968
[53]	M. Esposito, M. Tamietto, G. C. Geminiani, A. Celeghin, A subcortical network for implicit visuo-spatial attention: Implications for Parkinson's disease, Cortex, 141 (2021), 421–435. https://doi.org/10.1016/j.cortex.2021.05.003 doi: 10.1016/j.cortex.2021.05.003
[54]	Y. Z. Sun, Y. Ding, T. Y. Yu, C. Y. Chen, P. Li, X. G. Yang, et al., Effect of Tiaoshen acupuncture on anxiety after methamphetamine withdrawal (in Chinese), Chin. Acupunct. Moxibust., 42 (2022), 277–280. https://doi.org/10.13703/j.0255-2930.20210213-k0001 doi: 10.13703/j.0255-2930.20210213-k0001
[55]	R. K. Meleppat, C. R. Fortenbach, Y. F. Jian, E. S. Martinez, K. Wagner, B. S. Modjtahedi, et al., In vivo imaging of retinal and choroidal morphology and vascular plexuses of vertebrates using swept-source optical coherence tomography, Transl. Vision Sci. Technol., 11 (2022), 11. https://doi.org/10.1167/tvst.11.8.11 doi: 10.1167/tvst.11.8.11
[56]	R. K. Meleppat, K. E. Ronning, S. J. Karlen, M. E. Burns, E. N. Pugh, R. J. Zawadzki, In vivo multimodal retinal imaging of disease-related pigmentary changes in retinal pigment epithelium, Sci. Rep., 11 (2021), 16252. https://doi.org/10.1038/s41598-021-95320-z doi: 10.1038/s41598-021-95320-z
[57]	R. K. Meleppat, P. F. Zhang, M. J. Ju, S. K. Manna, Y. Jian, E. N. Pugh, et al., Directional optical coherence tomography reveals melanin concentration-dependent scattering properties of retinal pigment epithelium, J. Biomed. Opt., 24 (2019), 066011. https://doi.org/10.1117/1.Jbo.24.6.066011 doi: 10.1117/1.Jbo.24.6.066011
[58]	K. M. Ratheesh, L. K. Seah, V. M. Murukeshan, Spectral phase-based automatic calibration scheme for swept source-based optical coherence tomography systems, Phys. Med. Biol., 61 (2016), 7652–7663. https://doi.org/10.1088/0031-9155/61/21/7652 doi: 10.1088/0031-9155/61/21/7652

This article has been cited by:

1.	Jeongeun Son, Yuncheng Du, Model-Free Adaptive Control of the Failing Heart Managed by Mechanical Supporting Devices, 2022, 55, 24058963, 750, 10.1016/j.ifacol.2022.07.534
2.	Beatrice De Lazzari, Roberto Badagliacca, Domenico Filomena, Silvia Papa, Carmine Dario Vizza, Massimo Capoccia, Claudio De Lazzari, CARDIOSIM©: The First Italian Software Platform for Simulation of the Cardiovascular System and Mechanical Circulatory and Ventilatory Support, 2022, 9, 2306-5354, 383, 10.3390/bioengineering9080383
3.	José Emilio Traver, Cristina Nuevo-Gallardo, Inés Tejado, Javier Fernández-Portales, Juan Francisco Ortega-Morán, J. Blas Pagador, Blas M. Vinagre, Cardiovascular Circulatory System and Left Carotid Model: A Fractional Approach to Disease Modeling, 2022, 6, 2504-3110, 64, 10.3390/fractalfract6020064
4.	Jeongeun Son, Dongping Du, Yuncheng Du, 2021, Control of a Mechanical Blood Pump based on a Trade-off between Aortic Valve Dynamics and Cardiac Outputs, 978-1-7281-1179-7, 5555, 10.1109/EMBC46164.2021.9630864
5.	Xiangyu Wu, Yunpeng Zhang, Xiangyuan Zheng, Hongtao Liu, Houfu Wang, Numerical simulation for suction detection based on improved model of cardiovascular system, 2022, 77, 17468094, 103808, 10.1016/j.bspc.2022.103808
6.	Shahriar Sadatieh, Maryam Dehghani, Mohsen Mohammadi, Reza Boostani, Extremum-seeking control of left ventricular assist device to maximize the cardiac output and prevent suction, 2021, 148, 09600779, 111013, 10.1016/j.chaos.2021.111013
7.	Shulei Li, Donghai Jin, Xingmin Gui, Dynamic characteristic modeling of left ventricular assist devices based on hysteresis effects, 2023, 00104825, 106737, 10.1016/j.compbiomed.2023.106737
8.	Xiao-Cong Zhang, Qi Zhang, Gui-Fu Wu, Hai-Tao Hu, Ling Lin, Shuai Tian, Li-Ling Hao, Tong Wang, Evaluation of enhanced external counterpulsation for diabetic foot based on a patient-specific 0D-1D cardiovascular system model, 2024, 255, 01692607, 108333, 10.1016/j.cmpb.2024.108333
9.	Hongtao Liu, Shuqin Liu, A multi-objective physiological control system for continuous-flow left ventricular assist device with non-invasive physiological feedback, 2025, 100, 17468094, 107020, 10.1016/j.bspc.2024.107020
10.	Rajarshi Bhattacharjee, Shouvik Chaudhuri, Anindita Ganguly, 2023, Chapter 29, 978-981-99-2679-4, 331, 10.1007/978-981-99-2680-0_29

Reader Comments

Your name:*

Email:*
© 2023 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)

通讯作者: 陈斌, bchen63@163.com

1.
沈阳化工大学材料科学与工程学院沈阳 110142

Mathematical Biosciences and Engineering

3.9

Metrics

Article views(1977) PDF downloads(84) Cited by(1)

Preview PDF

Download XML

Export Citation

Article outline

Show full outline

Figures and Tables

Figures(8) / Tables(4)

Mathematical Biosciences and Engineering

Aberrant intrinsic functional brain topology in methamphetamine-dependent individuals after six-months of abstinence

Related Papers:

Abstract

1. Introduction

2. Materials and method

2.1. CLP Model

2.1.1. Ventricles

2.1.2. Atria

2.1.3. Valves

2.1.4. Systemic circulation and pulmonary circulation

2.2. Baroreflex model

2.2.1. Afferent neural pathway

2.2.2. Efferent neural pathways

2.2.3. Effectors

2.3. LVAD model

2.4. Simulation protocol

3. Results

3.1. Hemodynamic response under normal condition

3.2. Hemodynamic response under LV failure condition

3.3. Hemodynamic response under LVAD support

3.4. Effect of pump speed on hemodynamic variables

4. Discussion

5. Conclusion

Acknowledgments

Conflict of interest

Appendix A

References

This article has been cited by:

Reader Comments

通讯作者: 陈斌, bchen63@163.com

Metrics

Figures and Tables

Other Articles By Authors

Catalog

Mathematical Biosciences and Engineering

Aberrant intrinsic functional brain topology in methamphetamine-dependent individuals after six-months of abstinence

Related Papers:

Abstract

1. Introduction

2. Materials and method

2.1. CLP Model

2.1.1. Ventricles

2.1.2. Atria

2.1.3. Valves

2.1.4. Systemic circulation and pulmonary circulation

2.2. Baroreflex model

2.2.1. Afferent neural pathway

2.2.2. Efferent neural pathways

2.2.3. Effectors

2.3. LVAD model

2.4. Simulation protocol

3. Results

3.1. Hemodynamic response under normal condition

3.2. Hemodynamic response under LV failure condition

3.3. Hemodynamic response under LVAD support

3.4. Effect of pump speed on hemodynamic variables

4. Discussion

5. Conclusion

Acknowledgments

Conflict of interest

Appendix A

References

This article has been cited by:

Reader Comments

通讯作者: 陈斌, bchen63@163.com

Metrics

Figures and Tables

Other Articles By Authors

Related pages

Tools

Export File

Citation

Format

Content

Catalog