Agent-based model generating stylized facts of fixed income markets

We develop an agent-based model (ABM) of a financial market with multiple assets belonging either to the fixed income or equity asset classes. The aim is to reproduce the main stylized facts of fixed income markets with regards to the emerging dynamics of the yield curves. Our ABM is rooted in the market model of Kaizoji et al. (J Econ Behav Organ 112:289–310, 2015) formulated with two types of traders: the rational and risk-averse fundamentalist investors and the noise traders who invest under the influence of social imitation and price momentum. The investors involved in the present market model diversify their investments between a preferred stock equivalent to a perpetual bond and multiple bonds of selected maturities. Among those, a zero-coupon bond provides a constant rate of return, while the prices of the coupon-paying bonds are determined at each time step by the equilibrium between the investors’ demands and supplies. As a result, the ABM creates an evolving yield curve determined by the aggregate impact of the traders’ investments. In agreement with real markets, it also produces transient turbulent periods in the prices’ time series as well as a humped term structure of volatility. We compare the dynamics arising from different processes governing the risk-free rate with those of the historical US Treasury market. Introducing Vasicek’s model of interest rates to both synthetic and empirical rates demonstrates the capacity of our ABM in reproducing the main characteristics of the surface of autocorrelation of the volatilities of the yields to maturity of the US Treasury bonds for the selected time-frame.


Introduction
This contribution aims at unveiling the potential residing in the creation of agentbased models (ABMs) of the fixed income asset class involving fundamentalist and chartist traders. To the best of the authors' knowledge, we provide the first ABM of a financial market containing multiple assets belonging to the fixed income asset class. The value of this type of computational economic model lies in its capacity to produce transient market phenomena originated by the interactions of the individual agents. As pointed out in Sornette (2014), these models are freed from any stationary condition and can reproduce such periods of market unrest instigated by prices' rallies and their subsequent corrections. Examples of ABMs developed to describe stock markets can be found in Kyle (1985), Black (1986) and Samanidou et al. (2007), among others.
The literature provides a significant number of equilibrium models implemented to analyse the dynamics of yield curves and to forecast their evolution through the behaviour of the corresponding forward rates. Such analyses are, e.g. found in Nelson and Siegel (1987), Cox et al. (1985), Vasicek (1977), Duffie and Kan (1996), Dai and Singleton (2002), Duffee (2002) and Heath et al. (1992). These models have, however, proven limited with regards to their capacity to account for out-of equilibrium conditions corresponding to the ever-recurring extreme transient phenomena observed in reality. Moreover, they do not consider the relationship between the equity and fixed income asset classes. The market model developed below aims at laying a foundation towards closing this gap. We thus propose a new approach to generate stylized facts of fixed income markets by a market model including the two previously mentioned asset classes. As exposed in Kaldor (1961), stylized facts represent "broad tendencies" understood as empirical truths observed in the markets of interest. The focus is therefore drawn on the emerging dynamics associated with the yields to maturity and evolving as a function of the traders' investments updated at each time step in conformity with the rules developed in this novel type of ABM.
The present formulation derives from the market model proposed by Kaizoji et al. (2015) introducing two types of agents investing in two assets. The risk-free asset provides a constant rate of return and the remaining asset's price is subjected to the market clearing process applied at each time step. This risky asset also pays dividends determined from a stochastic multiplicative process. The agents are either fundamentalists or noise traders and do not change their strategies in time. This is in contrast to earlier ABMs in which the traders may switch between predictor (Arthur et al. 1996) or between chartist and fundamentalist strategies (De Grauwe et al. 1995;Brock and Hommes 1997). The fundamentalist strategy is as described in Brock and Hommes (1998) and in Chiarella et al. (2009). The chartist strategy is derived from the developments of Lux and Marchesi (1999). In a nutshell, their investments are updated probabilistically at each time step under the influence of other agents and of the momenta in the prices' time series. At the aggregated level, Kaizoji et al. (2015)'s model proved able to reproduce fat-tail distributions of returns, slow decaying autocor-relations of absolute returns, fast decaying autocorrelations of signed returns, volatility clustering and transient faster-than-exponentially growing prices associated with bubbles. A market model characterized by the interactions of fundamentalist and chartist traders through the price-vector of dimension higher than two is provided in Xu et al. (2014). Other models introducing heterogeneous agents investing in a multi-asset market can be found in Borghesi and Bouchaud (2007), Chiarella et al. (2007), Fedyk et al. (2013) and Eckrot et al. (2016). More recent articles revealing a renewed interest in agent-based modelling are proposed by Gualdi et al. (2015), Baghestanian et al. (2015), Bouchaud (2013) and Hommes and LeBaron (2018). Concerning the application to fixed income, Braun-Munzinger et al. (2018) in particular present a model involving three distinct agents investing in a unique corporate bond.
Section 2 introduces the market framework involving the fixed income and equity asset classes as well as the two trader types investing in this framework. The market clearing process is defined subsequently before the unveiling of the parameter selection and the initialization of the variables achieved in Sect. 3. The latter section concludes with the analyses of the time series generated in two representative simulations, the first one corresponding to a quiet market regime and the other one exhibiting the characteristics of a turbulent regime. Section 4 presents the analyses comparing the dynamics emerging from the market model with those of the US Treasury market realized between November 1993 and June 2020. The section starts with the implementation of Vasicek's model of interest rates to calibrate the synthetic risk-free rates of the simulations. It continues with the presentation of the term structures of volatility and surfaces of autocorrelations of the yields to maturity associated with each model set-up and with the US Treasury market for comparison. Section 5 concludes.

Model set-up
Our market model is made of a stock and of multiple bonds of selected maturities. One of the bonds is a zero-coupon bond providing a constant reference rate throughout the simulations. The agents belong either to the fundamentalist or to the noise trader type. The fundamentalist formulates expectations of the future returns as well as of the risk associated with each asset. The fundamentalist's investments are ruled by the need to maximize a constant relative risk aversion (CRRA) expected utility function. The noise traders' investments are subjected to social imitation, adherence to momentum trading and to the influence of a time-varying herding propensity parameter. Their individual investments are updated probabilistically at each time step. Each trader type is represented by a unique agent expressing the aggregate excess demands for each asset at each time step. The market clearing condition imposing the equilibrium between the supplies and demands for the assets determines the endogenous asset prices.
In the following, we start with presenting the market framework considering the fixed income asset class. The assets' and wealth dynamics are then derived. The description of both aggregated traders is then given, starting with the fundamentalist and the associated generalized optimization problem. The threefold investment process governing the chartist's investments is subsequently defined before the presentation of the market clearing achieved by the Walrasian auctioneer (Walras 1874). The latter process provides the system of nonlinear equations ruling the updates of the asset prices at each time step. The selection of the parameter values initializing the market model follows, before the presentation of the time series generated by a configuration of the market model considering six coupon paying bonds and a unique preferred stock.

The fixed income market framework
Our aim is to constrain the model to involve traders following only "aggressive" investment strategies. The selection of such traders is justified by their catalysing effect on the market through the permanent rebalancing of their portfolio. This is in contrast to passive investors who "freeze" their investments to collect the regular coupon payments and benefit from the yields offered by their securities. This strategy is known as "riding the yield curve".
• Assumption 1: The times to maturity of the bonds evaluated at t are assumed to stay constant in [t; t + t] for t sufficiently small. • Assumption 2: The investors are supposed to replicate at t the exact same portfolio constructed at t − 1 before updating their investments. As soon as t is such that the earlier approximation does not hold anymore, each aggregated trader sells the corresponding bonds and buys new ones with the appropriate maturities. This replication has no impact on the bonds' prices as the sell orders of an agent replicating her portfolio are assumed to be compensated by the buy orders of another agent doing the same process. The sole impact on the prices hence results from the excess demands arising from the update of the agents' investments. • Assumption 3: The auctions of the coupon paying bonds take place at each time step and are ruled by the Walrasian market making process. They are not reopened, have no impact on the coupon rates and do not distinguish between competitive and non-competitive bidders. • Assumption 4: The coupon rates are defined by an exogenous stochastic process reflecting the state of the economy. The dependency between the excess demands for the bonds and the coupon rates is ensured by the traders' investment processes. • Assumption 5: The outstanding of coupon paying bonds remains constant in time. Table 1 presents the different types of assets involved in the market model and their corresponding maturities. The preferred stock is granted an infinite maturity and is hence considered as a perpetual bond. The model involves M = 8 assets, the first one being the risk-free bond and the last one the preferred stock.
The yields to maturity (YTM), implied spot rates, yield duration, convexity and dispersions are defined as follows.
1. Yields to maturity The market model is developed such that the bonds' yields are derived from the associated prices. The prices are updated at each time step through the market making process, with the updates of the corresponding yields being simultaneous. The YTMs are defined as internal rates of returns equalizing the present values of the future cash flows (coupons and principal payments) with As the time to maturity of A M tends to infinity, the associated yield YTM M is evaluated by: The consideration of the preferred stock as a perpetual bond is further justified by considering Eq. (89) provided in Appendix A.1, which emphasizes the fact that the price of the associated zero-coupon bond vanishes as the time to maturity tends to infinity. 2. Implied spot rates The implied spot rates (ISR) associated with each bond are introduced successively and correspond to the holding period returns realized over the entire life of each bond. They are particular forward rates covering the maturities of the bonds and thus provide approximations of the yields to maturity with the advantage to allow omitting the consideration of the rates at which the coupons are reinvested. These rates are evaluated in terms of the corresponding discount factors, as presented in Appendix A.1. 3. Yield durations and convexities The market model considers both modified (ModDur k t ) and Macaulay (MacDur k t ) durations associated with each asset. The Macaulay duration is defined as in Marrison (2002). Equations (99) to (101) of Appendix A.1 provide the approximations of these durations as well as of the corresponding convexities (Conv k t ). 4. Portfolio yields, durations, convexities and dispersions The aggregated portfolio statistics are defined successively in Appendix A.1, considering each portfolio as a unique entity equivalent to a bigger bond paying bundles of cash flows at each time step. Two approaches are applied to approximate the yields of the traders' portfolio: the market-value-weighted (MV-weighted) and basis-point-value-weighted (BPV-weighted) approaches. The latter replaces the traders' wealth fractions by the proportion of each asset's BPV to the sum of the BPVs of the assets in which the portfolios are invested. The BPVs are expressed as the change in market value of each asset originated by a variation of its yield by a basis point, i.e. by 0.01%. The modified market-value-weighted average portfolio and Macaulay durations, convexities and dispersions are further defined in the appendices. This concludes the presentation of the parameters involved in the fixed income framework.

Assets and wealth dynamics
This section defines the assets' returns and how the traders' wealth are updated. A particular attention is paid to the accrued interests in the returns of the coupon-paying bonds of finite maturities.
• Assets' returns -Risk-free asset As in Kaizoji et al. (2015), the risk-free asset provides a constant rate of return r f . A 1 is further defined as a zero-coupon bond having a maturity of one year. The selected "day-count" convention defines a year to correspond to 250 time steps. The issuance price PV 1 and face value FV 1 of A 1 are related by r f as follows: -Preferred stock The preferred stock A M pays periodical dividends to its holders. Its "one time step" return r M t consists of the sum of the price return where r d is the long-term growth rate, σ d the standard deviations of the multiplicative process and u M t ∼ N (0, 1) an i.i.d. random variable. -Coupon-paying bonds of finite maturities Assuming the absence of taxation and the simultaneous payments of the coupons and dividends, the returns realized on each bond A k for k ∈ {2, . . . , M − 1} are expressed as the sum of the price and coupon returns as follows: In particular, the coupon payment d k t inherently corresponds to the interest accrued between the purchase of the bond A k at t − 1 and its sale at t. This arises from the replication process exposed earlier and requiring to buy another bond at t having the same maturity as the bond bought at t − 1. The accrued interests are related to the annual coupon payments of C k t through: d k t = C k t × dt/250, where the time increment is set to dt = 1. The multiplicative process governing the determination of the accrued coupon payments is further expressed as: where σ k and ν k t , respectively, represent the standard deviation of the process and the associated i.i.d. random variable ∼ N (0, 1).

• Wealth dynamics
The wealth dynamics are expressed in terms of the wealth fractions of the corresponding traders and of the assets' returns described above. Let W i t represent the wealth level of trader i ∈ { f , c} at time step t. Also, let the number of shares z i,k t of asset k held by the corresponding trader at t be expressed in terms of the wealth fraction The wealth of the trader allocated to this asset is updated as: Summing over the 7 coupon-and dividend-paying assets, one obtains: The following trivial conditions are moreover applied on the agents' wealth fractions due to the absence of both borrowing and short selling: The fraction of wealth x i,r f t invested in the risk-free asset is thus: By summing the wealth increments given by Eq. (8) with the one associated with the risk-free asset, one can express the overall wealth dynamics between two consecutive time steps as: which is equivalent to the formulation proposed hereafter, considering the expressions of the total returns exposed earlier: where the coupon and dividend payments associated with each asset A k for k ∈ {2, . . . , 8} are: This concludes the presentation of the wealth updates and leads to the introduction of the fundamentalist's investment strategy.

Fundamentalist trader
The fundamentalist trader is defined in accordance with the descriptions provided in Chiarella et al. (2009), Kaizoji et al. (2015) and Xu et al. (2014). The trader is granted a constant relative risk-aversion (CRRA) utility function required to be maximized w.r.t. the wealth allocations. This objective is represented by a generalized myopic mean-variance optimization problem expressed in terms of the ex ante expectation of the utility function depending on the asset returns and wealth allocations introduced previously. Accounting for Eq. (12), the optimization problem proposed in Kaizoji et al. (2015) is generalized to the multi-asset case hereafter.
• Generalized optimization problem where U (W ) ∈ R is the CRRA utility function and x f t the vector containing the wealth fractions associated with each asset at t.
The derivation of the solution is provided in Appendix A.2 and consists in solving the Hamilton-Jacobi-Bellman partial differential equation (PDE) associated with the generalized optimization problem of concern. The equivalent formulation of the optimization problem as a stochastic optimal control problem is proposed below. This is achieved by considering Eq. (12) to define the wealth update in terms of a stochastic differential equation (SDE). To do so, the individual returns r k t+1 of each asset A k for k ∈ {2, . . . , 8} are expressed as in Xu et al. (2014) in terms of the expectations E t r k t+1 formulated at t of the returns to be realized at t +1, of the standard deviations σ k t and of the i.i.d. random variables ψ k t such that ψ k t ∼ N (0, 1): Repeating, we use discrete time equations and the slight abuse of notation of representing one time step by the symbol dt (which is equal to 1) and the corresponding one-step increment of the wealth of the fundamentalist trader is denoted dW f t . Considering Eqs. (12) and (15), the SDE accounting for the stochasticity of the wealth level of the fundamentalist is hence expressed as: which can be synthesized by the successive formulation proposed in Xu et al. (2014), where μ t represents the drift coefficient, σ t the standard deviation and t ∼ N (0, 1) an i.i.d. random variable. Considering the unit vector e ∈ R M−1 and the transpose operator (.) , one can further express μ t (W f ) and σ t (W ) in both scalar and vector notations: where we distinguish σ t (W f ) and σ t , the former accounting for the diffusion term of Eq. (17) and the latter representing the vector containing the standard deviations of the individual assets' returns. Combining Eqs. (17), (18) and (19), the resulting SDE is given by: where t represents a one-dimensional Brownian motion and where the covariance matrix is expressed as , one can finally write: The cost functional subjected to the previous wealth dynamics is subsequently introduced as follows: The problem hence requires to find the optimal control x ∈ R M−1 maximizing the cost functional subjected to the wealth defined by the stochastic system of Eq. (21).
• Generalized stochastic optimal control problem: The solution is given by the optimal state trajectory x * : [0, T ] → R M−1 . Introducing the risk aversion parameter γ = 1 (γ > 0 always holds), the cost-to-go function is given by: Consequently, the stochastic optimal control problem is expressed as follows, considering the differential form of the standard Brownian motion (t) ∈ R: subjected to: The solution of the generalized stochastic optimal control problem is obtained in Appendix A.2 and reads The remaining unknowns consist of the covariance matrix and of the vector containing the expectations of the future returns. The matrix = σ σ is defined to be diagonal and the individual variances remain constant. One obtains the following notations introducing the Kronecker delta δ kl : The expectations of the future returns are taken as the sum of a constant and a variable term for each of the yields to maturity and implied spot rates. This allows us to account for the trader's own opinion about the returns provided by each bond (constant term) while still considering the reality of the market (variable term). The constant terms consist of the initial values of the parameters. The variable terms are those evaluated at t − 1 as a consequence of the replication process explained in Sect. 2.1. Therefore, for each bond A k where k ∈ {2, . . . , M − 1}, one expresses: Fig. 1 Representation of the market as seen by the chartists: differentiation of the two pools of assets and corresponding steps of the investment process. The arrows represent the transit of the agents' wealth where ω 1 , ω 2 and ω 3 are the weights, respectively, associated with the initial values of the yields to maturity and implied spot rates as well as with their evaluations at t − 1. As explained in Appendix A.1, the implied spot rates cannot be computed for the preferred stock. Moreover, due to the particular identity of the asset, the varying term associated with the yield to maturity is neglected. For k = M, one hence obtains: The wealth fractions of the fundamentalist trader are ultimately expressed as: with the boundary condition given by The subsequent developments introduce the investment process of the chartist trader taking the other side of the trades.

Noise traders
As in Kaizoji et al. (2015), the individual chartist traders do not diversify their allocations but instead select an individual asset in which they invest their entire wealth. Their investments are updated probabilistically at each time step through an Ising model-like set-up under the influence of the momentum in the assets' prices' time series and the opinion of other noise traders. At the aggregated level, the representative noise trader allocates her wealth proportionally to the number of individual traders invested in the corresponding assets.
As shown in Fig. 1, the noise trader identifies two mutually exclusive pools of assets. The first one contains all the bonds including the risk-free one and is identified as P. The second one is solely constituted of the preferred stock A M . Built on this market representation, three distinct steps are defined below to decompose the investment process of the aggregated chartist. Figure 1 further illustrates which particular pool is concerned by each of these steps. The first one tackles the rebalancing of the aggregated chartist's portfolio between the two asset classes considered. The second one is concerned with the update of the investments of the individual agents inside the bonds of P. The third one finally reconciles the two previous updates simultaneously in order to derive the wealth fractions of the aggregated trader.
Step 1: Aggregated chartist's portfolio rebalancing between P and A M : The first step of the investment process of the aggregated chartist deals with the update of the allocations in either of the two groups of assets defined previously. It derives an estimate of the number of individual traders transiting between the two asset classes. Both social imitation and adherence to a momentum-following strategy influence the outcome of the present step.
The number of individual chartists invested in A M is N M and the number of remaining agents is N 1 M−1 . The total fraction of wealth invested in A M at t is hence: Following the definition of Kaizoji et al. (2015), the associated opinion index s c t is defined as: The parameter accounting for momentum is moreover defined as: introducing the momentum associated with the preferred stock: and the one related to P as: Parameter θ controls the time scale ∼ 1 1−θ over which the momentum is estimated. Equation (36) introduces the simple "rule of thumb" followed by the individual chartists to compare the price returns of the isolated asset with the "term-spreads" provided by P. As shown in Eq. (34), this spread is defined as the difference between the yields to maturity of the bonds situated at the extremities of the yield curve, fitted to the "one-time-step" returns of the preferred stock. This spread moreover represents a shorthand evaluation of the steepness of the yield curve giving insights about potential future returns obtained in excess of those given by the risk-free bond. Considering the fact that YTM 1 t = cst., the higher the yields of the longer-term bond, the more the aggregated chartist is influenced to allocate her wealth in the fixed income asset class.
As exposed in Kaizoji et al. (2015), the herding propensity κ t is introduced in order to represent the noise traders' susceptibility to herding and their propensity to adhere to a momentum following strategy. This parameter is moreover defined to be either constant or to follow a stochastic mean-reverting process. In the latter case, it incorporates the alternating regimes of pessimistic mood and exuberance observed in financial markets, as detailed in Shiller (2015) and in Sornette (2003). In the former case, it steers the model dynamics by accounting for stable market regimes. Section 2.6 provides a thorough definition of the two distinct processes ruling the dynamics of the parameter. Given the aforementioned rationale for its introduction in the model, the herding propensity is directly impacting the noise traders' dynamics through the transition probabilities presented in Eq. 37 and involving the opinion index and momentum parameter.
By definition, the transition probabilities p + t and p − t are associated with the following actions: an individual noise trader holding the preferred stock A M at t shifts her investment to P with a probability p + t and another noise trader invested in P at the same moment transits to A M with a probability p − t . These probabilities are consequently expressed as: where p + = cst. and p − = cst. control the average holding time associated with each group of assets. Setting p − > p + implies that an up-or downward increment of either of H step1 t or s step1 t is not engendering the same reaction of the agents. This is translated by a favour given to A M to the detriment of P.
Introducing the Bernoulli random variables ξ k ( p) taking the value of 1 with a probability of p and 0 otherwise, one can express the subsequent master equations: The net number of agents transiting between both pools of assets at t is finally estimated by the absolute value of the variable N +→− t defined as: Step 2: Aggregated chartist's investments within P: The second step of the allocation process tackles the allocations of the chartists inside P. Their updates are achieved independently for each pair of neighbouring assets (A i , A i+1 ) and (A 1 , A M−1 ). Again, both social imitation and adherence to a momentum-following strategy influence the outcome of these updates. The opinion index associated with the pair (A i , A i+1 ) is defined below, fitting the variable to the momenta introduced subsequently and introducing N i t as the number of chartists invested in A i at t: The momentum parameter associated with the same pair of assets is moreover expressed as: introducing the momentum parameters associated with each asset. These are further expressed in terms of the sub-parameters accounting for the momenta observed in the yields and in the prices: where α ∈ [0, 1]. In particular: where the yields are fitted to the "one time step" periods over which the price returns are determined. Equation (44) is moreover introducing the consideration of decreasing yields by the chartists. This enforces the assessment of positive realized price returns and contrasts with the previous step in which the agents valued higher yields implying higher potential price returns when looking at the overall market. Once invested inside the fixed income asset class, these agents naturally favour decreasing yields corresponding to positive realized price returns.
In accordance with the formulation of Eq. (37), the transition probabilities associated with each pair of consecutive assets are defined as: where p i,i+1 t + quantifies the probability to invest in A i when starting from A i+1 . Accordingly, p i,i+1 t − accounts for the probability to invest in A i+1 when starting from A i .
The number of agents transiting between two neighbouring assets in either direction is hence: Ultimately, the net number of agents leaving or entering each bond i is expressed as: Step 3: Execution of the update of the aggregated chartist's investments: The last step executes the update of the investments on the aggregated level. The following opinion index is first introduced for each of the M − 1 bonds of fixed maturities: The performance indicators accounting for the popularity of each bond and for the trend observed in their returns are defined as follows, considering the parameters H i t introduced in Eq. (42): Building on these parameters, the categorical probabilities are further defined as where Perf worst t and Perf best t are, respectively, min i∈{1,...,M−1} Perf i t and max i∈{1,...,M−1} Perf i t . The weights ω i t associated with each bond are finally defined as: As a matter of fact, ω i t corresponds to the share of the expected value of the Bernoulli trial associated with A i and endowed with a probability of success of p step3 t i over the sum of each of the M − 1 expectations. These weights are introduced to adjust the dispersion of the agents inside P as can be inferred from the expression of Eq. (52) presented below.
The number of agents transiting either to or from each bond of P is subsequently introduced as N step3 t i . By considering the result of the update achieved in the first step of the allocation process, one can hence write: where . represents the floor operator dropping the decimal part. With this variable, one can execute the update of the individual investments considering the outcome of the two previous steps: Accordingly, the updated number of agents invested in the preferred stock is given by: As a result, the updated share of wealth of the aggregated chartist trader invested in each asset A i for i ∈ {1, . . . , M} is ultimately expressed as: This concludes the presentation of the aggregated noise trader and leads to the formulation of the market making process.

Market clearing process
The present section tackles the market clearing process executed at each time step by the Walrasian auctioneer. First, the expressions of the excess demands formulated by both aggregated traders are developed. The system of nonlinear equations resulting from the satisfaction of the Walrasian equilibrium condition is provided subsequently.
On the one hand, the excess demand of the fundamentalist trader D where A i t and B i t are obtained from Eq. (31) presenting the expression of the wealth fractions associated with each asset: As a reminder and considering Eq. (29), i t for i ∈ {2, . . . , M −1} is further expressed as: And for i = M, M t is: On the other hand, the excess demand of the aggregated chartist trader associated with the same asset is expressed as: As has been exposed in the second step of the allocation process of the chartist trader presented in Sect. 2.4, the returns of the individualized bonds are evaluated from the auction prices. These being further set to P i auction = 1, the return of each bond realized at t is hence P i t /P i auction − 1 = P i t − 1. A M being further exempted from any auction, its returns are directly evaluated by P M t /P M t−1 . As a consequence, the following rule is applied in Eq. (61): For all i ∈ {2, . . . , M}, the Walrasian equilibrium condition is finally expressed as: The system of nonlinear equations involving the unknown prices remaining to be evaluated at each time step is finally obtained by inserting Eqs. (56) and (61) in Eq. (63). The resulting system is presented below and the detailed developments leading to this formulation are provided in Appendix A.3. The equation determining the prices reads where α ∈ R (M−1)×(M−1) , β ∈ R M−1 , ζ ∈ R M−1 , χ ∈ R (M−1)×(M−1) and λ ∈ R M−1 are defined as: (65) Ultimately, this system of M − 1 equations can be synthesized hereafter by the functions f i (t) ∈ R: An initial guess of the solution is provided by the previously realized prices (P 2 t−1 , . . . , P i t−1 , . . . P M−1 t−1 ) and is further refined by the following process introducing P, the vector notation of the prices of the coupon and dividend paying assets: where J is the Jacobian matrix defined by J i j = ∂ f i /∂ P j . This concludes the presentation of the market clearing process.

Dynamics of the herding propensity Ä t
Recall that the herding propensity κ t represents the noise traders' susceptibility to herding and their propensity to adhere to a momentum following strategy. It is defined such as to have a direct impact on the transition probabilities p + t and p − t given by expression (37) in Sect. 2.4. Following Kaizoji et al. (2015), we endow κ t with a dynamics capturing the changing nature of economic and social conditions that influence the behaviour of noise traders. As economic and social conditions are arguably not systematically constant, it is reasonable to consider market regimes corresponding to κ = cst., thereby introducing a source of stochastic variability inside the transition probabilities through the varying herding propensity. On the other hand, economic and social conditions happen to evolve relatively slowly during certain periods and are characterized by tendencies, trends and variability that can be considered approximately stationary. This scenario is accounted for by κ = cst., where the influence of the herding propensity on the transition probabilities remains steady. In fact, fixing the herding propensity to a constant value is analogous to neutralizing the exogenous influence of the social and economic conditions, thereby rendering the model a closed autonomous system. Hence, we will mainly consider the case where κ t follows a

Statistics
Instantaneous evaluation Long-run evaluation stochastic mean-reverting process defined as where the initial condition is given by κ 0 and where ν t is an i.i.d. discrete-time white noise process with a mean of 0 and a standard deviation of 1. The mean-reversion rate is further represented by η k , the mean value by μ k and the diffusion associated with the Wiener process ν t by σ k . The instantaneous and long-run evaluations of the expectation and variance of the herding propensity are successively given in Table 2 above. The process of concern is hence found to be stationary and normally distributed in the long-run, such that the following expression holds: The statistical estimation of the parameters of Eq. 68 is finally achieved as in Kaizoji et al. (2015) via the following expressions: where T is the time frame during which κ t is set to revert when being in the supercritical regime identified by the bounds located two standard deviations away from the mean.

Parameter selection
The initialization of the parameters introduced in the market model is subsequently exposed hereafter. The analyses of the time series generated in typical simulations ensue. The wealth fractions of the aggregated traders are first initialized below, starting with the fundamentalist trader. Given the initial wealth fraction x f ,M 0 , the constant variance σ 2 M and the initial expectation of the future return μ M 0 = E M r ,t +d 0 (1+r d )/P 0 associated with A M , the constant risk aversion parameter γ is defined a priori by the following expression: The expressions of the initial wealth fractions associated with the coupon paying bonds A i for i ∈ {2, . . . , M − 1} are obtained accordingly: by ensuring that the following condition is satisfied: The expectations of the future returns μ i 0 associated with these assets are defined below and the constant values of the variances σ 2 i are provided in Table 3. On the other hand, the initial wealth fractions of the aggregated chartist are obtained from the value of x c,M 0 through the following expression: Considering the floor operator . introduced earlier and the total amount of chartists N , the amount of individual agents initially invested in the coupon and dividend paying assets is hence given by: This further enables to initialize the opinion indices involved in the allocation process governing the aggregated chartist's investments: Standard deviation of the dividend processes Initial prices P 0 = (1, 1, 1, 1, 1, 1, 1) (4) (5) Initial coupons or dividend d 0 = (5, 5, 5, 5, 5, 5, 12) × 10 −5 Weight associated with the latter momenta Constants involved in the transition probabilities (37) Herding propensity Long-run average μ k = 0.98 Mean-reversion rate η k = 0.05 Diffusion of the associated Wiener process Table 3 lists the parameters and initial values of the variables involved in the market model. The values included in parentheses, respectively, correspond to each of the assets in terms of increasing order of maturities. The initialization of the herding propensity is achieved as in Kaizoji et al. (2015).
The annualized yields to maturity are successively initialized as follows: One can hence observe that the initial yield curve is flat as each yield initially equals 2.5%, with the exception of the yield of the perpetual bond initialized as YTM M 0 = 3%. This precise set-up is justified by the need to analyse the impacts of the agents' investments on a neutral configuration of the yield curve. This further explains the selection of the face values of the corresponding assets.
Building on these developments, the evaluations of the initial artificial changes of the market values, modified and Macaulay durations, convexities, discount factors and implied spot rates of each coupon-paying bond are straightforward. The according values are presented in Table 4 of Appendix A.4.
Finally, given the previous developments, the expectations of the future returns are initialized as: (81) This concludes the presentation of the initialization of the constant parameters and variables defined previously and leads to the subsequent analyses of the phenomena emerging from numerical simulations of the market model.

Time series analyses
The following section specifies the dynamics emerging in typical simulations. Figures Fig. 2 Time dependence of the variables associated with the market model. The simulation length is set to T sim = 5000 time steps and the herding propensity remains constant throughout the simulation and can be observed through bubbles in the corresponding prices' time series. These two different types of simulations are analysed hereafter.
The time series of the asset prices P i t show the divergence of the preferred stock's prices from those of the bonds. This divergence is governed by the growth rate of r d . The subcritical herding behaviour associated with the constant nature of the herding propensity results solely in subtle deviations from the average trajectory given by the previous rate. In contrast, the simulations for which the herding propensity varies are found to generate transient regimes when the preferred stock prices undergo turbulent patterns. Such periods are, e.g. observed for t ∈ [1300, 1700]. Analyses of the superexponential growth characterizing these bubbles are provided in Kaizoji et al. (2015), among others. Besides their prices, the bonds are characterized by their yields to maturity YTM i t . In both simulations, the yields fluctuate around their initial values with an amplitude related to the maturities of the bonds. These fluctuations reach higher values for bonds of low maturities. This constitutes an emerging phenomenon referred to as the humped "term structure of volatility". Moreover, according to the well-known price-yield relationship, the yields of the perpetual bond reach their lowest values during the peak of the price bubbles, e.g. for t ∈ [1400, 1700] in Fig. 3. The reverting of the prices of A M to levels directed by the growth rate of r d as shown t ∈ [3000, 3100] symmetrically drives the yields back to their initial value of 3%.
The noise traders switch their investments between the different assets with a probability expressed in terms of the momentum of the prices of the preferred stock and of the opinion of others. The dynamics of the momentum H M t of the preferred stock, of the probabilities p ± t (s  replicates the main trends of the preferred stock's prices. This is due to the high amplitudes of the exponentially weighted moving average of the returns of A M in comparison with the stable spreads proposed by the yield curve during the simulations. The latter amplitudes are nevertheless found to be of approximately an order of magnitude lower than those associated with the opinion index s step1 t . This explains the strong in influence given to the latter variable in the outcome of the Bernoulli trials. The transition probabilities fluctuate symmetrically around ( p + + p − )/2 and are bounded in [0, ≈ 0.15] for the simulation governed by a constant herding propensity. The smooth fluctuations observed in this simulation further explain the lack of abrupt changes in the preferred stock's price dynamics, as can be expected from the outcome of the Bernoulli trials involving these probabilities and generating no significant rebalancing of the aggregated chartist's portfolio between the two pools of assets. In contrast, the transition probabilities observed in the simulation involving the stochastic herding propensity increase sharply during the emergence or burst of the bubbles. The noise traders' opinion formation is illustrated by the following pairs of assets selected to avoid any redundancy in the analyses and to assess the integration of A 1 in P: (A 1 , A 2 ), (A 7 ,A 1 ) and (A 3 ,A 4 ). All three indices are found to fluctuate around a long-run average of zero with the same amplitudes, indicating an equilibrated dispersion of the agents in P. The risk-free asset is hence not introducing any bias in the allocation process. This statement is emphasized by the fact that the trends in the fluctuations of the indices s step2 t | 1,2 and s step2 t | 7,1 replicate each other and are approximately symmetrically mirrored by those of s step2 t | 3,4 . The emphasis is furthermore directed towards the ascertainment of the abrupt saturation of the opinion indices observed during periods of herding towards the preferred stock proposed by the simulation involving a varying herding propensity. Considering the fact that a bad performing bond is penalized twice: the first time by the loss of agents transiting to the neighbouring bonds and the second time by the outflow of agents towards A M , it appears that a very small amount of agents remain invested in the bonds during the culmination of the preferred stock's price bubbles. As a result, the opinion indices become very sensitive to the transition of agents as dictated by the outcome of the Bernoulli trials, explaining their abrupt variations at the corresponding time steps.
The transition probabilities define the aggregated noise trader's investment decisions as they determinate the update of the number of traders invested each asset. , H step1 t ) = 0. In this case, the reaching of a threshold value by the herding propensity satisfies the statement proposed in Sornette (1994) and granting the responsibility of the "sweeping of an instability" of the latter parameter to instigate a self-reinforcing loop leading to a gradual feedback between H M t and s step1 t provoking an outflow of agents from P to A M . At the apogee of this phenomenon, no agent is spurred to transit to the fixed income portfolio. This situation is, however, unsustainable and the flux of agents between the two asset classes is reversed as soon as p + t > 0, leading to the inevitable burst of the bubble. The market model is moreover found to be subjected to an initial phase of auto-regulation occurring at the beginning of each simulation. The quintessence of this observation lies in the increase of x c,M from 30% to a stable level of ≈ 50% on the long-run. As a consequence of the resulting increase of the preferred stock's price level, x f ,M is further reduced at the benefit of x f ,1 . This is due to the counter-cyclical reaction of the fundamentalist trader to the evolution of the market. More generally, this "price taker" characteristic is found to be revealed in the time series of the corresponding wealth fractions mirroring the fluctuations of the prices of A M .
Figures 4 and 5 provide a comparison of the fundamentalist's and aggregated chartist's portfolio statistics for simulations, respectively, involving a constant and a mean-reverting herding propensity. They show the evolution of the average portfolio yields AvgYieldMV i t and AvgYieldBPV i t , aggregate wealth levels W i t as well as approximated portfolio modified durations ApproxModDur i t , convexities ApproxConv i t and dispersions AvgDisp i t . First, the average market-value-weighted and basis-point-value-weighted portfolio yields to maturity (AvgYieldMV i t and AvgYieldBPV i t ) are analysed. The portfolio yields computed from the two different approaches provide different assessments of the attractivity of each investment strategy. The MV-weighted yields tend to replicate the patterns of the yields of the preferred stock less than is achieved by the BPVweighted yields. This is striking for the fundamentalist trader who constrains her endowments to A M during the simulations. The BPV-weighted yields contrast with the other yields in accordance to the expression of Eq. (113). In particular, as the perpetual bond is granted a much higher duration than the other bonds, its prices are The average portfolio modified durations AvgModDur i t tend to replicate the dynamics of the modified durations associated with the traders' predominant investment. The fundamentalist is hence found to have an almost stable portfolio modified duration over time with only few exceptions arising during the transient market unfolding. As a result, the latter variable is hence fluctuating between 20 and 30 years during periods of quietude. The portfolio modified duration of the fundamentalist is also more sensitive to increases of the perpetual bond's duration than to the according decrease of x f ,M t . This is illustrated by the humps in the time series of the corresponding variable generated in the simulation involving the varying herding propensity. The present market model moreover provides a striking observation during periods of extreme transient phenomena and unveiling the limitation of the well-known fixed income pricing formula expressed in Eq. (102). In fact, the relationship implying, among others, that simultaneous increases of the duration and decreases of the yields lead to an increase in the corresponding asset's market value is found not to hold in turbulent periods for the portfolio of the fundamentalist trader. As one can observe during any of the price bubbles associated with A M , the wealth of the fundamentalist decreases despite the fact that the portfolio modified duration increases and that the portfolio yield decreases. The convexity is besides found to increase, adding no interference to the previous statement. This sheds lights into the boundedness of the earlier mentioned pricing formula and explains why one might loose commonly used landmarks when managing a fundamentalist portfolio during periods of market unrest. Consequently, surfing partly on the price bubble might help the fundamentalist to contain the losses caused by the massive exodus of the chartists from the other assets to the perpetual bond subjected to the bubbles of concern. Following the tendencies observed for the portfolio modified durations, the average portfolio convexities AvgConv i t and dispersions AvgDisp i t propose patterns strikingly sensitive to the dynamics of the perpetual bond's prices. One can again observe the fact that the rebalancing executed by the fundamentalist as a reaction to the price deviations of A M is not sufficient to counteract the overall tendencies of the parameters of concern to follow the dynamics of the latter asset.
The comparison of the normalized wealth levels W f t of the fundamentalist and W c t of the aggregated chartist reveals the fact that the fundamentalist's strategy outperforms the chartist's one irrespectively of the process governing the evolution of the herding propensity. The notion of opportunity costs is however crucial to evaluate the reversal of the latter trends on the short-term during turbulent periods. As a matter of fact, such costs appear to be incurred by the fundamentalist during the build-up of the bubbles of the prices of A M . Nevertheless, the subsequent bursts instigated by the flee of individual chartists towards P bringing the variable x c,M t to anterior levels happen to cancel all the gains of the aggregated chartist, while the fundamentalist kept benefiting from the constant price return of the risk-free asset. As a result, the decreases of W f observed during the build-up of the price bubbles are adjusted by the still invariable price returns granted by A 1 left unimpaired by the variations of x f ,1 t and x c,1 t . This explains the slight increase of the difference of the wealth levels of the aggregated traders after each burst of a bubble.

Application of the Vasicek model of interest rates
The previous developments introduced a constant risk-free rate provided by the zerocoupon bond of one-year maturity. As a matter of fact, the rates provided by such risk-free assets are not constant in reality. The solution of the Vasicek model of interest rates developed by the eponym author in Vasicek (1977) is applied in the market model to provide a varying exogenous arbitrage-free risk-free rate. As a result, the comparison of the emerging dynamics associated with this new set-up with those of the one considering r f = cst. shall demonstrate the capacity of the new model to reproduce strikingly similar autocorrelations of the volatilities of the yields to maturity in comparison with those observed in the US Treasury market. The statistical estimation of the parameters involved in the expression of the varying risk-free rate is first achieved below, before the presentation of the emerging dynamics arising in typical simulations. The analyses of the autocorrelations are presented in Sect. 4.2.2.
The "one time step" rates are determined from the following closed-form expression: where ν t ∼ N (0, 1) is an i.i.d random variable and λ, θ and σ r are estimated in the same way as achieved for the herding propensity. The long-run average of the risk-free rate is equal to the constant value of the rate applied earlier, namely with θ = 8×10 −5 . The mean-reversion strength is subsequently estimated according to: where r max = 1 × 10 −4 and r critic = 9.5 × 10 −5 . The value attributed to θ is known and T = 60 time steps. One can hence evaluate the mean-reversion strength as λ ≈ 2 × 10 −3 . Finally, σ r is obtained as a function of the mean-reversion rate as where the variance of the risk-free rate is such that Var r t = 2 × 10 −5 , twice as low as the variance of the perpetual bond's returns. The initial value of the risk-free rate is finally set to r f ,0 = θ . Figure 6 presents the time series generated by a typical simulation involving a varying herding propensity as well as the stochastic process governing the evolution of r f ,t . The evolution of the risk-free rate has a strong reflection on x f ,1 , while the preferred stock's prices show a clear symmetric pattern of its dynamics. The evolution of the wealth levels further illustrates the indirect impact of the varying risk-free rate. The reason for the high sensitivity of the fundamentalist's wealth fractions to the variations of r f ,t is straightforward when considering Eqs. (30) and (31). In this regard, the fundamentalists' investments are hence found to have a major impact on the assets' prices, being mainly driven by the evolution of the risk-free rate and having a significant impact on the wealth levels of both aggregated traders. There is, however, no form of self-reinforcing feedback arising between the fundamentalist's investments and the corresponding impacts on the assets' prices, the latter agent having no incentive to follow series of positive returns and being predominantly reacting to the stochastic variations of the risk-free rate. This further explains why one cannot find clear trends in the evolution of the difference of both agent's normalized wealth levels, especially after the periods of extreme transient phenomena, where the earlier mention about the opportunity costs is not holding anymore as one is not guaranteed to obtain a stable  (respectively, mean-reverting) return in the money market instrument during the culmination of the preferred stock's price bubbles.

Dynamics of the fixed income parameters
The following developments focus on the emerging dynamics of the yields to maturity, implied spot rates and discount factors. For the sake of clarity, the l.h.s. (respectively, r.h.s.) of each illustration is associated with the simulation involving the constant herding propensity (respectively, varying).

Time dependencies of the fixed income parameters
Figures 7, 8 and 9 present the surfaces generated from the time series associated with the yields to maturity, implied spot rates and discount factors obtained from the two typical simulations presented in Sect. 3.2. The yields of the perpetual bond are represented separately in Fig. 7. Their distinction from the other yields is achieved for the sake of clarity as the corresponding asset is omitted in the computation of the yield curves. An ensemble of T sim = 5000 yield curves is generated in each simulation, i.e. one yield curve per time step. The impacts of the aggregated agents' investments are straightforward considering the initial flat and neutral configuration of the yield curve. The increasing volatile character of the yields to maturity in terms of decreasing maturities is striking in both simulations and further reveals the aptitude of the present market model to produce a humped term structure of volatility. This paradigm is subsequently tackled in Sect. 4.2.2.
The implied spot rates defined according to Eq. (97) provide accurate approximations of the yields to maturity. As a result, the surfaces presented in Fig. 8 are very similar to those of Fig. 7. The spikes appearing at t ∈ [2300, 2600] in the former very accurately replicate those emerging at the same period in the latter illustration. Figure 9 finally shows the surfaces generated by the linear interpolations of the discount factors associated with each of the bonds of fixed maturities. As one can observe, the This concludes the presentation of the evolution of the yields to maturity, implied spot rates and discount factors. Figure 10 presents the surfaces of volatilities of the yields to maturity. The volatilities are computed over moving windows of 125 time steps in both simulations. They are evaluated at each time step from the following definition,

Volatilities of the yields to maturity
where the sample average operator is represented by . and the selected moving window by τ = 125. Figure 11 shows the realized volatilities of the yields of US Treasury bonds evaluated between the 18th of November 1993 and the 6th of June 2020, considering moving windows of 250 business days. The selection of this moving window is justified by the need to amplify the trends of the fluctuations of the realized volatilities. Notwithstanding the fact that the volatilities of the yields of the US Treasury bonds are subjected to the exogenous policies of the Federal Reserve, the comparison between those realized endogenously in the present model with those of the US Treasury market enables to assess the capacity of the present market model to account for the main trends observed in reality.
The model produces a humped term structure of volatility in both simulations. The term structure of volatility is further decreasing when considering the bonds of maturities larger than one year. Such phenomena are found episodically in the US Treasury market, e.g. between 2004 and 2008, corresponding to t ∈ [3200, 3950] in Fig. 11. On the other hand, the exogenous imposition of low discount rates by the Federal Reserve are commonly known to be reflected on the yields of the US Treasury bonds of short maturities, explaining the boundedness of their volatilities during the corresponding time frames (t ∈ [4000, 5500]). The increases of these bounds hence  (2016) and Bertocchi et al. (2005), the volatilities of the US Treasury yield curve are found to be following such a hump-like shape on average, comforting the assertion of the capacity of the present market model to account for the main trends observed for the volatilities of the yields of US Treasury bonds.

Autocorrelations of the volatilities of the yields to maturity
Figures 12 and 13 present the surfaces generated from the coefficients of autocorrelation of the volatilities of the yields to maturity computed for each bond, excluding the risk-free one and including the perpetual one. The former and the latter, respectively, consider the simulations involving a constant and a varying risk-free rate. The surfaces are generated by computing the average of the coefficients of autocorrelations from a set of 300 simulations. The first 500 steps of each simulation are furthermore omitted in the computations due to the auto-regulation phenomenon commented in Sect. 3.2. The above-mentioned coefficients are determined for each simulation from Eq. (86), considering the volatilities σ YTM i ,t determined earlier and the lags l ∈ [0, 500]: The comparison between Figs. 12 and 13 reveals the fact that the surfaces associated with the market model involving a varying risk-free rate decay much faster than the others in terms of increasing lags. This emerging phenomenon is hence attributed to the stochastic character of the evolution of the risk-free rate. There is moreover no difference on the aggregated level for different selections of the processes governing the evolution of the herding propensity. The comparison of the surfaces presented in Figs. 12 and 13 with the surface of autocorrelation associated with the US Treasury market provided in Fig. 14 and computed from the earlier-mentioned volatilities further enables to distinguish the capacity of the simulations involving the varying risk-free rates to replicate the same emerging phenomena as they reproduce similar autocorrelation coefficients at the corresponding lags. More precisely, the latter simulations are capable to replicate the steeper gradients of the autocorrelation surfaces until the 250th lag as well as the milder ones for higher lags directing the surfaces towards null autocorrelations at the 500th lag.
Equation 86 defining the autocorrelations of the volatilities of the yields to maturity aims to uncover persistent dynamics revealing a self-dependency of the observed variable. For example, high autocorrelation values indicate that periods of high volatility induce high volatilities at a later point of time. As one can expect, letting the risk-free rate follow an independent stochastic process affects the dynamics of the yields to maturity through the transmission channel of the noise trader's investments presented in step 2 of Sect. 2.4. As a result, the residual self-dependencies of the aforementioned variable are found to diminish at a higher pace when increasing the lags, as one can expect the autocorrelations to be dampened. This is consistent with the empirical observations provided in Fig. 14. However, the latter illustration reveals the fact that the yield volatilities of the bonds of shorter maturities exhibit larger auto-correlations than those of longer maturities. This stylized fact is not reproduced by the present model, as one would need to consider the classical theories of the term structure of interest rates in the traders' investment processes. According to these theories, fixed income investors choose maturities based on their expectations of future short-term interest rates in order to maximize their overall horizon yield. The relative merits of expectations theory, segmented market theory or liquidity preference theory remains to be decided, as none of these theories have yet gained the upper hand, as exposed in Smith (2014).

Conclusion
The agent-based market model proposed in this research has extended the model proposed in Kaizoji et al. (2015) by introducing several bonds of fixed maturities. The inclusion of the various fixed income parameters allowed us to describe the investment processes of both aggregated traders to the peculiarities of the market in which they evolve. The market set-up initialized by a flat yield curve demonstrated its capacity to reproduce several stylized facts of the US Treasury market. The simulations involving a varying herding propensity showed that the well-known bond-pricing formula involving the durations, convexities and shifts in yields to maturity does not hold at the fundamentalist's portfolio level during periods of market unrest. Subsequently, the implementation of Vasicek's model of interest rates to simulate the evolution of the risk-free rate was found to allow for the linearly dependent fundamentalist's wealth fractions to replicate its dynamics and have a significant impact on the asset prices. The analyses of the volatilities of the yields to maturity further enabled to assess the capacity of the market model to produce a humped term structure of volatility as is observed on average in the US Treasury market. Finally, the market model introducing the mean-reverting process governing the dynamics of the risk-free rate also proved superior to the one involving r f = cst. with regards to its capacity to reproduce the main characteristics of the surface of autocorrelation of the volatilities of the yields to maturity of the US Treasury bonds for selected the time-frame.

Conflict of interest
The authors declare no conflicts of interest and certify that they have no affiliations with or involvement in any organization or entity with any financial interest or non-financial interest in the subject matter or materials discussed in this manuscript. No funding was received for conducting this study.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

• Security level
The following developments introduce the parameters mentioned in Sect. 2.1. First and foremost, the famous relationship between an asset's no-arbitrage price PV, face value FV, coupon payments C, time to maturity N and yield to maturity YTM is expressed as: Representing the discount factor associated with the payments occurring at maturity by δ(N ), the no-arbitrage price can be further expressed as: exacerbating the expression of the bonds' dirty prices in terms of those of a perpetuity and of a zero-coupon bond as follows: PV = PV perpetuity (1 − δ(N )) + δ(N )PV zero-coupon bond .  ApproxModDur 0 (0. 980583, 4.76207, 9.25537, 13.5134, 17.5685, 21.4533, 25.2023, 33.3704) ApproxMacDur 0 (1, 4.88112, 9.48676, 13.8513, 18.0078, 21.9896, 25.8324, 34.3715) ApproxConv 0 (1. 92309, 27.6518, 97.5912, 205.408, 347.599, 521.555, 725.624, 2224. By isolating the nonzero terms when inserting the Kronecker delta in the summations, one can finally express Equation (137) as: This concludes the presentation of the solvable set of M − 1 nonlinear equations expressed in terms of the M − 1 asset prices remaining to be evaluated at each time step of the simulations of the model. Table 4 below introduces the initial values of the parameters associated with the fixed income market framework. As one could anticipate, the initial values of the modified and Macaulay durations approximate the maturities of the bonds.