Next Article in Journal
Optimization of the Real-Time Response to Roadside Incidents through Heuristic and Linear Programming
Next Article in Special Issue
Dynamic Multicriteria Game with Pollution Externalities
Previous Article in Journal
Novel Methods in Multiple Criteria Decision-Making Process (MCRAT and RAPS)—Application in the Mining Industry
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Dynamic Multi-Objective Duopoly Game with Capital Accumulation and Pollution

1
Department of Economics, Université Panthéon-Assas, Paris II & CRED, EA 7321, 21 Rue Valette, 75005 Paris, France
2
Department of Econometrics and Operations Research & Center, Tilburg University, LE 5000 Tilburg, The Netherlands
3
Department of Economics, University of Antwerp, 2000 Antwerp, Belgium
*
Author to whom correspondence should be addressed.
Mathematics 2021, 9(16), 1983; https://doi.org/10.3390/math9161983
Submission received: 26 July 2021 / Revised: 14 August 2021 / Accepted: 15 August 2021 / Published: 19 August 2021
(This article belongs to the Special Issue Multicriteria Games)

Abstract

:
This paper studies a discrete-time dynamic duopoly game with homogenous goods. Both firms have to decide on investment where investment increases production capacity so that they are able to put a larger quantity on the market. The downside, however, is that a larger quantity raises pollution. The firms have multiple objectives in the sense that each one maximizes the discounted profit stream and appreciates a clean environment as well. We obtain some surprising results. First, where it is known from the continuous-time differential game literature that firms invest more under a feedback information structure compared to an open-loop one, we detect scenarios where the opposite holds. Second, in a feedback Nash equilibrium, capital stock is more sensitive to environmental appreciation than in the open-loop case.

1. Introduction

Corporate social responsibility (CSR) emphasizes that a firm does not just care about profit maximization, but has at the same time other objectives in mind that benefit society. According to Wirl et al. [1], examples of CSR projects are environmental reports, philanthropic support and sponsoring, energy and environmental management, mentoring and educational programs for workers, family-friendly workplaces and more. Wirl et al. [1] studied a dynamic model where a firm optimizes over CSR activities. Lambertini et al. [2] consider CSR activities in a strategic setting. In particular, a Cournot game is designed where one firm, as usual, maximizes profits, but where the other firm is a CSR firm in the sense that consumer surplus and pollution are also taken into account. Lambertini et al. [2] reach the surprising result that, provided the market is sufficiently large, the CSR firm earns higher profits. Yanase [3] also considers CSR activities but in that sense solely concentrates on the environment. The study investigates how an increase in the firms’ environmental consciousness affects the environment and economic welfare. Feichtinger et al. [4] take a similar approach by analyzing a dynamic oligopoly in which environmental externalities are taken into account in the firms’ objectives.
Our paper studies a dynamic duopoly in which the firms accumulate capital. The positive side is that with the capital stock firms can produce goods which can be sold on the market, resulting in firm profits. The negative side is that the production process of both firms causes pollution. As in the literature just mentioned, the firm takes both these elements explicitly into account. The new element of our research, however, is that we take a multi-objective approach, in that at the same time firms strive to maximize profits and minimize pollution (the different weights between the two objectives are not a priori fixed). In so doing, we obtain the new result that, compared to an open-loop information structure, in a feedback Nash equilibrium the firms are less sensitive to changes in environmental appreciation.
Our approach is close to albeit different from that of Rettieva [5]. She considers a dynamic, discrete time, two-player game model where the players use a common resource. Both players have two goals to optimize over an infinite horizon, namely maximizing profit from selling fish and minimizing catching costs (the value of those costs depends on both players’ catches). (She also studies the case where the players have different random horizons.) She introduces an original version of the Nash bargaining scheme (one for each goal), where status quo points are computed in different cases. (Those cases refer to different instances of zero-sum games.) Multi-criteria cooperative equilibria are considered by Rettieva [6].
Our contribution differs in several ways from that of Rettieva. First, we are interested in a pollution-control problem, and not in a resource-management problem. Second, we do not rely on the Nash bargaining approach; we follow an approach that can be loosely called Nash equilibrium with Pareto equilibrium. More precisely, in equilibrium each firm chooses a Pareto-optimal solution to its multi-objective problem. Third, among the objectives pursued by the players, one is especially environmental (whereas we may think that in Rettieva’s setting players actually have two economic objectives). (This paper is related to, but also differs from, the papers by Kuzyutin, Gromova and Pankratova [7,8] who develop a cooperative game theory approach of dynamic multi-criteria decision analysis and who pay special attention to the time-consistency issue).
Our paper is also connected to the work of Crettez and Hayek (2021) [9] who study a dynamic Cournot duopoly where production causes pollution, albeit without capital accumulation. They notably show that contrary to the case where the firms’ objective just concerns profit maximization, there exist Nash equilibria where production is lower than in the cooperative equilibrium. (In our setting, when firms have a unique objective, cooperating means behaving like a monopolist. A monopolist usually reduces production in order to raise the price of its product.) This stems from the fact that in a Nash equilibrium firms do not coordinate on the choice of the relative weight given to the environmental objective. This result highlights the fact that when firms pursue multiple and separate objectives, equilibrium behavior can display an over-reaction. In this paper, we also obtain this form of over-compliance. Here, however, open-loop equilibria differ from feedback equilibria, whereas in the work of Crettez and Hayek [9], open-loop equilibria are all feedback.
In addition to the literature that studies corporate social responsibility within a dynamic framework, our paper is also connected to straightforward capital accumulation models where firms maximize profits. A seminal paper in that area is that of Reynolds (1987) [10] who finds, as a major result, that firms invest more under a feedback information structure than under an open-loop one. If firms put the environmental weights to zero, we in fact have the same problem as in the work of Reynolds [10]. However, there is a difference, because where Reynolds [10] studies a continuous time framework, our model is in discrete time. Very surprisingly, it turns out that we find scenarios where in the long run the open-loop equilibrium admits steady state capital stocks that are bigger than in the feedback case.
Other approaches of multiple objective games and Stackelberg games can be found in recent publications such as [11,12,13,14,15].
The paper is organized as follows. Section 2 presents the model. Section 3 and Section 4 develop the open loop and the feedback Nash equilibria, respectively. The cooperative equilibrium is analyzed in Section 5. Section 6 compares the different equilibrium outcomes and Section 7 concludes the paper.

2. Model

Consider a dynamic duopoly where firms produce the same type of goods and where the inverse demand function at date t is given by
p ( q 1 t , q 2 t ) = max 0 , a b ( q 1 t + q 2 t ) , a > 0 , b > 0 .
In the expression above, q i t denotes the quantity supplied by firm i at date t, a is the maximum willingness to pay for the product and b represents the sensibility of the market price to the volume of production brought to the market. (We shall always assume that p ( q 1 t , q 2 t ) is strictly positive. We can also interpret a as being the maximum willingness to pay minus the average production cost (in that case, we assume that this average cost is constant)).
For simplicity, production is proportional to the stock of capital K i t held by each firm i at date t, and we set the proportionality coefficients equal to one: q i t = K i t . The law of motion of the capital stock is as follows
K i t + 1 = I i t + ( 1 δ ) K i t ,
where I i t is the investment rate of firm i at date t, and δ is a depreciation rate ( 0 < δ < 1 ). (This is a discrete-time variant of the model proposed by Reynolds [10] (see also Lambertini (2018) [16] p. 73).) We assume that it costs I i t 2 2 for firm i to decrease/increase its capital stock by an amount I i t (in the event that I i t < 0 , we have I i t ( 1 δ ) K i t ).
Production is a polluting activity. To simplify the analysis, we suppose that emissions are equal to total production, i.e., q 1 t + q 2 t , and that the dynamics of the pollution stock x t is written as follows
x t + 1 = η x t + q 1 t + q 2 t ,
where 0 < 1 η < 1 is the natural decay rate of pollution, and x 0 , the initial value of pollution, is given.
Firm i’s objective is given by the following multi-objective program.
max ( I i t ) t t = 0 β t a b K i t + K 3 i t K i t I i t 2 2 , 0 β t x t
subject to
K i t + 1 = I i t + ( 1 δ ) K i t ,
x t + 1 = η x t + K i t + K 3 i t ,
I i t ( 1 δ ) K i t ,
x 0 , K i 0 , and ( K 3 i t ) t being given ,
In the program above, ( K 3 i t ) t is a bounded sequence such that 0 K 3 i t a for all t. The meaning of the program above is that firm i tries to maximize the value of its profit 0 β t ( a b ( K i t + K 3 i t ) K i t I i t 2 2 ) where β is a discount factor such that 0 < β < 1 , and to minimize 0 β t x t , a measure of intertemporal pollution. For simplicity again, in this measure future values of pollution are discounted at the same rate as future receipts.

3. Open-Loop Nash Equilibria

We next state the definition of an open-loop equilibrium for our dynamic game.
Definition 1.
A Nash equilibrium is a sequence of decisions ( I ^ 1 t , I ^ 2 t ) t such that ( I ^ i t ) t is a solution to problems (3)–(7) when ( I ^ 3 i t ) t is given, i = 1 , 2 .
Notice that the sequence ( I ^ 3 i t ) t determines the sequence of ( K ^ 3 i t ) t that is supposed to be given in problem (3)–(7).
Before presenting optimality conditions for an open-loop Nash equilibrium, notice that in any equilibrium each firm chooses a Pareto-optimal solution to its multi-objective problem. (This is why some authors call Nash equilibria for multi-objective games Pareto-Nash equilibria. See, e.g., Qu and Ji (2016) [17] Section 1.2 for a quick literature review.) As a consequence, in any equilibrium each firm gives a relative weight to its environmental objective. In the sequel, we will denote by λ 2 the relative weight of the environmental objective of firm 1 and by γ 2 the corresponding relative weight of firm 2. Of course, not all pairs of weights ( λ 2 , γ 2 ) will be admissible in equilibrium.

3.1. Optimality Conditions

The next result determines the set of equilibrium weights given by the two firms to their environmental objectives as well as the equilibrium investment rates.
Proposition 1.
For each pair ( λ 2 , γ 2 ) of relative weights given to the environmental objectives of firms 1 and 2, respectively, satisfying the following conditions
β λ 2 1 β η 2 b + δ ( 1 β 1 + δ ) b β γ 2 1 β η < a b + δ ( 1 β 1 + δ ) ,
β γ 2 ( 1 β η ) 2 b + δ ( 1 β 1 + δ ) b β λ 2 ( 1 β η ) < a b + δ ( 1 β 1 + δ ) ,
there is an open-loop Nash equilibrium ( I 1 t , I 2 t ) t such that I i t > ( 1 δ ) K i t , i = 1 , 2 , with
I 1 t = K 10 2 r 1 t r 1 ( 1 δ ) + s 1 t s 1 ( 1 δ ) + D 1 1 s 1 t + 1 ( 1 δ ) ( 1 s 1 t ) + K 20 C 2 r 1 t r 1 ( 1 δ ) s 1 t s 1 ( 1 δ ) .
I 2 t = K 20 2 r 1 t r 1 ( 1 δ ) + s 1 t s 1 ( 1 δ ) + D 2 1 s 1 t + 1 ( 1 δ ) ( 1 s 1 t ) + K 10 C 2 r 1 t r 1 ( 1 δ ) s 1 t s 1 ( 1 δ ) ,
where
C = 2 a ( λ 2 + γ 2 ) β 1 β η δ ( 1 + δ + 1 / β ) + 3 b ,
D 1 = a b C λ 2 β 1 β η δ ( 1 + δ + 1 / β ) + b , D 2 = a b C γ 2 β 1 β η δ ( 1 + δ + 1 / β ) + b ,
r 1 = ( 1 δ ) 2 + 3 b + 1 / β Δ 2 ( 1 δ ) ,
s 1 = ( 1 δ ) 2 + b + 1 / β Δ 2 ( 1 δ ) ,
Δ = ( ( 1 δ ) 2 + 3 b + 1 / β ) 2 4 ( 1 δ ) 2 / β
Δ = ( ( 1 δ ) 2 + b + 1 / β ) 2 4 ( 1 δ ) 2 / β .
The equilibrium values of and K 1 t , K 2 t , x t are given by
K 1 t = K 10 ( r 1 t + s 1 t ) 2 + D 1 ( 1 s 1 t ) + ( K 20 C ) 2 ( r 1 t s 1 t ) ,
K 2 t = K 20 ( r 1 t + s 1 t ) 2 + D 2 ( 1 s 1 t ) + ( K 10 C ) 2 ( r 1 t s 1 t ) ,
x t = η t ( x 0 K 10 + K 20 C r 1 η C 1 η ) + r 1 t K 10 + K 20 C r 1 η + C 1 η .
Proof. 
From Theorem 4.2 in the work of Hayek [18], we know that if a bounded sequence ( I 1 t ) t is a solution to problem (3)–(7) such that I 1 t > ( 1 δ ) K 1 t , then there exist λ 1 0 , λ 2 0 , ( μ t ) t in 1 (the set of sequence ( x t ) t such that t = 0 | x t | is a real number), ( ψ t ) t in 1 , not all nil, such that any each date, ( I 1 t ) t maximizes the following Hamiltonian, given ( I 2 t , K 2 t ) t :
H ( I 1 t , K 1 t , x t , λ 1 , λ 2 , μ t , ψ t ) = λ 1 ( β t ( a b ( K 1 t + K 2 t ) K 1 t I 1 t 2 2 ) + λ 2 ( β t x t ) + μ t + 1 ( η x t + K 1 t + K 2 t ) + ψ t + 1 I 1 t + ( 1 δ ) K 1 t
The first-order conditions are given by
0 = λ 1 β t I 1 t + ψ t + 1 ,
ψ t = λ 1 β t a 2 b K 1 t b K 2 t + μ t + 1 + ( 1 δ ) ψ t + 1 ,
μ t = λ 2 β t + η μ t + 1 .
The last two equations refer to the dynamics of the shadow values ψ t and μ t of the capital stock and pollution, respectively.
Note that since ( μ t ) t in 1 , ( ψ t ) t in 1 , the familiar transversality condition is necessarily satisfied.
Since H is concave in ( I t , K t ) and since the familiar transversality condition is satisfied, the first-order conditions above are sufficient. We now study these conditions in detail.
Observe that λ 1 0 . Otherwise, we would have ψ t + 1 = 0 for all t, which in turn would yield μ t + 1 = 0 for all t, so that λ 2 = 0 . All those variables would be nil, which is impossible.
By a standard argument, we obtain
μ t = λ 2 β t 1 β η .
The following dynamic system obtained from Equations (1), (2) and (15)–(17) remains to be solved. (Its solution will be an open-loop Nash equilibrium since the Hamiltonians are concave and the transversality condition is satisfied).
ψ t = λ 1 β t a 2 b K 1 t b K 2 t λ 2 β t + 1 1 β η + ( 1 δ ) ψ t + 1
K 1 t + 1 = ψ t + 1 λ 1 β t + ( 1 δ ) K 1 t ,
x t + 1 = η x t + K 1 t + K 2 t .
Since λ 1 > 0 , we can set λ 1 = 1 . Defining ψ t β t = Ψ t , Equations (19) and (20) can be written as
Ψ t = a 2 b K 1 t b K 2 t β λ 2 1 β η + β ( 1 δ ) Ψ t + 1 ,
K 1 t + 1 = β Ψ t + 1 + ( 1 δ ) K 1 t ,
Combining the two equations above, we obtain
K 1 t ( 1 δ ) K 1 t 1 β = a 2 b K 1 t b K 2 t β λ 2 1 β η + β ( 1 δ ) K 1 t + 1 ( 1 δ ) K 1 t β .
Rearranging, we obtain
( 1 δ ) K 1 t + 1 K 1 t ( 1 δ ) 2 + 2 b + 1 / β + ( 1 δ ) K 1 t 1 / β = a + b K 2 t + β λ 2 1 β η
or
( 1 δ ) K 1 t + 2 K 1 t + 1 ( 1 δ ) 2 + 2 b + 1 / β + ( 1 δ ) K 1 t / β = a + b K 2 t + 1 + β λ 2 1 β η
This equation also holds for firm 3 i . By considering symmetry and taking the optimality conditions for both firms into account, we obtain
( 1 δ ) ( K 1 t + 2 + K 2 t + 2 ) ( K 1 t + 1 + K 2 t + 1 ) ( 1 δ ) 2 + 3 b + 1 / β + ( 1 δ ) ( K 1 t + K 2 t ) / β
= 2 a + β ( λ 2 + γ 2 ) 1 β η .
Set Z t = K 1 t + K 2 t . We obtain:
( 1 δ ) Z t + 2 Z t + 1 ( 1 δ ) 2 + 3 b + 1 / β + ( 1 δ ) Z t / β
= 2 a + β ( λ 2 + γ 2 ) 1 β η
Define:
Δ = ( 1 δ ) 2 + 3 b + 1 / β 2 4 ( 1 δ ) 2 / β = ( 1 δ ) 2 + 3 b + 1 / β 2 ( 1 δ ) / β ( 1 δ ) 2 + 3 b + 1 / β + 2 ( 1 δ ) / β > 0
since the expression ( 1 δ ) 2 + 3 b + 1 / β 2 ( 1 δ ) / β > 0 , being a decreasing function of β which is positive for β = 1 .
Solving the above difference equation, we obtain Z t = c 1 r 1 t + c 2 r 2 t + C , where
C : = 2 a + ( λ 2 + γ 2 ) β 1 β η ( 1 δ ) ( ( 1 δ ) 2 + 3 b + 1 / β ) + ( 1 δ ) / β = 2 a + ( λ 2 + γ 2 ) β 1 β η ( 1 δ ) ( ( 1 δ ) 2 + 3 b ) δ / β = 2 a + ( λ 2 + γ 2 ) β 1 β η δ ( 1 δ 1 / β ) 3 b = 2 a ( λ 2 + γ 2 ) β 1 β η δ ( 1 + δ + 1 / β ) + 3 b .
Moreover,
r 1 = ( 1 δ ) 2 + 3 b + 1 / β Δ 2 ( 1 δ )
= 2 ( 1 δ ) β [ ( 1 δ ) 2 + 3 b + 1 / β + Δ ]
r 2 = ( 1 δ ) 2 + 3 b + 1 / β + Δ 2 ( 1 δ )
Z t = c 1 r 1 t + c 2 r 2 t + C .
Observe that r 1 r 2 = 1 / β > 1 and r 2 > 1 . We also have r 1 < 1 .
To show that r 1 < 1 , we shall prove that ( 1 δ ) 2 + 3 b + 1 / β Δ < 2 ( 1 δ ) . We have
( 1 δ ) 2 + 3 b + 1 / β 2 ( 1 δ ) < Δ
( ( 1 δ ) 2 + 3 b + 1 / β 2 ( 1 δ ) ) 2 < Δ ( since ( 1 δ ) 2 + 3 b + 1 / β 2 ( 1 δ ) > 0
( ( 1 δ ) 2 + 3 b + 1 / β ) 2 4 ( 1 δ ) ( ( 1 δ ) 2 + 3 b + 1 / β ) + 4 ( 1 δ ) 2 < ( ( 1 δ ) 2 + 3 b + 1 / β ) 2 4 ( 1 δ ) 2 / β
δ ( 1 δ ) ( 3 b + 1 / β ) < ( 1 δ ) / β
δ ( 1 δ 1 / β ) 3 b < 0
In this case, since Z t is bounded, we have c 2 = 0 and Z t = c 1 r 1 t + C .
Notice that Z t > 0 if
C > 0 λ 2 + γ 2 < 2 a ( 1 β η ) β .
The right condition holds since it is implied by conditions (8) and (9).
Then, Z t = c 1 r 1 t + C = ( Z 0 C ) r 1 t + C = ( K 10 + K 20 C ) r 1 t + C = ( K 10 + K 20 ) r 1 t + C ( 1 r 1 t ) > 0 .
Now, setting K 2 t + 1 = Z t + 1 K 1 t + 1 in Equation (25), we obtain
( 1 δ ) K 1 t + 2 K 1 t + 1 ( 1 δ ) 2 + b + 1 / β + ( 1 δ ) K 1 t / β = a + b Z t + 1 + β λ 2 1 β η
or
( 1 δ ) K 1 t + 2 K 1 t + 1 ( 1 δ ) 2 + b + 1 / β + ( 1 δ ) K 1 t / β = a + b c 1 r 1 t + 1 + b C + β λ 2 1 β η .
Solving the above difference equation we obtain K 1 t = d 1 s 1 t + d 2 s 2 t + D t where
D t = a + b C + λ 2 β 1 β η ( 1 δ ) ( ( 1 δ ) 2 + b + 1 / β ) + ( 1 δ ) / β
+ b c 1 ( ( 1 δ ) r 1 2 ( ( 1 δ ) 2 + b + 1 / β ) r 1 + ( 1 δ ) / β ) r 1 t + 1 .
By definition of r 1 , we have
( 1 δ ) r 1 2 ( 1 δ ) 2 + 3 b + 1 β r 1 + ( 1 δ ) β = 0 .
Therefore, the denominator in the fraction above reduces to 2 b r 1 and thus
D t = a b C λ 2 β 1 β η δ ( 1 + δ + 1 / β ) + b + c 1 2 r 1 t = D 1 + c 1 2 r 1 t
where
D 1 = a b C λ 2 β 1 β η δ ( 1 + δ + 1 / β ) + b .
Notice that
Δ = ( 1 δ ) 2 + b + 1 / β 2 4 ( 1 δ ) 2 / β > 0 s 1 = ( 1 δ ) 2 + b + 1 / β Δ 2 ( 1 δ ) = 2 ( 1 δ ) β [ ( 1 δ ) 2 + b + 1 / β + Δ ] s 2 = ( 1 δ ) 2 + b + 1 / β + Δ 2 ( 1 δ ) .
We observe again that s 1 s 2 = 1 / β > 1 , s 2 > 1 and s 1 < 1 .
In this case, since K 1 t is bounded we have d 2 = 0 and K 1 t = d 1 s 1 t + D t and we have K 2 t = Z t K 1 t . More precisely, d 1 = K 10 D 0 and
D 0 = D 1 + c 1 2 ,
where c 1 = ( K 10 + K 20 C ) .
We therefore obtain:
K 1 t = d 1 s 1 t + D t = ( K 10 D 0 ) s 1 t + D 1 + c 1 2 r 1 t = ( K 10 ( D 1 + c 1 2 ) ) s 1 t + D 1 + c 1 2 r 1 t = ( K 10 D 1 ) s 1 t + D 1 + c 1 2 r 1 t c 1 2 s 1 t = K 10 s 1 t + D 1 ( 1 s 1 t ) + ( K 10 + K 20 C ) 2 ( r 1 t s 1 t ) = K 10 ( r 1 t + s 1 t ) 2 + D 1 ( 1 s 1 t ) + ( K 20 C ) 2 ( r 1 t s 1 t ) .
We now must make sure that K 1 > 0 . This implies that D 1 > 0 . Using the expressions of D 1 and C, we obtain that D 1 > 0 if and only if:
β λ 2 1 β η 2 b + δ ( 1 + δ + 1 β ) b β γ 2 1 β η < a b + δ ( 1 + δ + 1 β ) .
which is condition (8). We obtain K 2 t by subtracting K 1 t from Z t .
K 2 t = ( K 10 + K 20 ) r 1 t + C ( 1 r 1 t ) K 10 ( r 1 t + s 1 t ) 2 + D 1 ( 1 s 1 t ) + ( K 20 C ) 2 ( r 1 t s 1 t ) = K 20 ( r 1 t + s 1 t ) 2 + C D 1 ( 1 s 1 t ) + ( C ) 2 ( r 1 t + s 1 t ) + ( K 10 ) 2 ( r 1 t s 1 t ) = K 20 ( r 1 t + s 1 t ) 2 + ( C D 1 ) ( 1 s 1 t ) + ( K 10 C ) 2 ( r 1 t s 1 t ) = K 20 ( r 1 t + s 1 t ) 2 + D 2 ( 1 s 1 t ) + ( K 10 C ) 2 ( r 1 t s 1 t ) .
where D 2 = C D 1 = a b C γ 2 β 1 β η δ ( 1 + δ + 1 / β ) + b (where we have used Equations (26) and (32)).
We also have to make sure that K 2 > 0 . This implies that C > D 1 . Using the expressions of C and D 1 , this condition can be written as
β γ 2 ( 1 β η ) 2 b + δ ( 1 β 1 + δ ) b β λ 2 ( 1 β η ) < a b + δ ( 1 β 1 + δ ) .
which is condition (9).
Finally, we obtain x t + 1 = η x t + K 1 t + K 2 t = η x t + Z t = η x t + ( K 10 + K 20 C ) r 1 t + C so
x t = η t x 0 K 10 + K 20 C r 1 η C 1 η + r 1 t K 10 + K 20 C r 1 η + C 1 η .
The proposition above shows that there are multiple interior Nash equilibria, each characterized by a pair of relative weights ( λ 2 , γ 2 ) given their environmental objectives.
Condition (8) defines an upper bound for the environmental appreciation parameter λ 2 such that the firm wants to be an active producer in the long run, translating into K 1 t > 0 . Analogously, condition (9) is a similar condition for γ 2 with respect to the long run level of K 2 . Furthermore, conditions (8) and (9) imply
λ 2 + γ 2 < 2 a ( 1 β η ) β ,
This condition ensures that C > 0 , so that in the long run K 1 t + K 2 t is positive.
The above condition is clear: if both firms give too much weight to their environmental objective, then they do not need to accumulate capital, because doing so will lead to an increase in pollution. Condition (8) will hold for large a, because then the market is profitable. It will not hold for large λ 2 , because firm 1 does not like to pollute and therefore reduces investments. It will hold for large γ 2 . Then, firm 2 has less incentive to invest, which increases the output price for firm 1, making investment for firm 1 more profitable.

3.2. Sensitivity Analysis

Lemma 1.
Consider two open-loop Nash equilibria with different relative weights λ 2 given to the environmental objective of firm 1 and the same relative weight γ 2 given to the environmental objective of firm 2. Then, for every time instant, the capital stock of the first firm K 1 t is always lower in the equilibrium associated with the higher value of λ 2 .
Proof. 
We show below that K 1 t λ 2 < 0 .
Using Equation (34), we have
K 1 t λ 2 = D 1 λ 2 ( 1 s 1 t ) C λ 2 ( r 1 t s 1 t ) 2
where
C = 2 a ( ( λ 2 + γ 2 ) β 1 β η ) δ ( 1 + δ + 1 / β ) + 3 b , D 1 = a b C ( λ 2 β 1 β η ) δ ( 1 + δ + 1 / β ) + b .
Set G = δ ( 1 δ 1 / β ) b < 0 .
C λ 2 = β ( 1 β η ) [ δ ( 1 + δ + 1 / β ) + 3 b ] < 0 ,
D 1 λ 2 = [ b C λ 2 + β 1 β η ] / G .
When t is large, the sign of K 1 t λ 2 is the same as that of D 1 λ 2 . We obtain:
D 1 λ 2 = β ( 1 β η ) G 1 b δ ( 1 + δ + 1 β ) + 3 b < 0
since
1 > b δ ( 1 + δ + 1 β ) + 3 b .
Moreover, this negative effect exists both in the short and the long run (because C λ 2 ( r 1 t s 1 t ) 2 < 0 , since one can check that r 1 < s 1 ). Therefore, we always have K 1 t λ 2 < 0 . □
Notice that this is in contrast with the model without capital accumulation of Crettez and Hayek [9] where a change in λ 2 has an indeterminate effect on the production level.
Lemma 2.
Consider two open-loop Nash equilibria with different relative weights γ 2 given to the environmental objective of firm 2 and the same relative weight λ 2 given to the environmental objective of firm 1. Then, for every time instant, the capital stock of the first firm K 1 t is always higher in the equilibrium associated with the higher value of γ 2 .
Proof. 
Using Equation (34), we have
K 1 t γ 2 = D 1 γ 2 ( 1 s 1 t ) C γ 2 ( r 1 t s 1 t ) 2
We have
C γ 2 = β ( 1 β η ) [ δ ( 1 + δ + 1 / β ) + 3 b ] < 0 ,
D 1 γ 2 = 1 G b C γ 2 < 0 ,
where we recall that G = δ ( 1 δ 1 / β ) b < 0 . Since r 1 < s 1 , the Lemma follows. □
Notice that we obtain the same result as when there is no capital. Both results are somewhat connected. When γ 2 becomes larger, firm 2 invests less (which is predicted by Lemma 1 but now with application to firm 2) which increases the output price, making investment for firm 1 more profitable. We next study feedback Nash equilibria.

4. Feedback Nash Equilibria

Let us assume that the investment rate of firm 3 i is given by a function I 3 i t ( K i t , K 3 i t , x t ) that depends on the state variables (feedback rule). Then, firm 3 i ’s capital accumulation reads
K 3 i t + 1 = I 3 i t ( K i t , K 3 i t , x t ) + ( 1 δ ) K 3 i t ,
with K 3 i 0 also being given.
We set firm i’s problem as follows
max ( I i t ) t 0 β t a b ( K i t + K 3 i t ) K i t I i t 2 2 , 0 β t x t
subject to
K i t + 1 = I i t + ( 1 δ ) K i t ,
K 3 i t + 1 = I 2 ( K i t , K 3 i t , x t ) + ( 1 δ ) K 3 i t ,
x t + 1 = η x t + K i t + K 3 i t ,
K i t 0 , x 0 , K i 0 , K 3 i 0 .
This problem motivates the following definition.
Definition 2.
A feedback Nash equilibrium is a pair of decision rules ( I ^ 1 ( K 1 , K 2 , x ) , I ^ 2 ( K 1 , K 2 , x ) ) such that the sequence ( I i t ) t determined by I i t = I ^ i ( K 1 t , K 2 t , x t ) is a solution to problems (47)–(51) when I 3 i = I ^ 3 i , i = 1 , 2 .
Notice that any decision rule for the investment rate must select a Pareto-optimal solution to problems (47)–(51). As was observed in the study of the open-loop equilibrium, in any feedback equilibrium each firm gives a relative weight to its environmental objective. We shall denote again by λ 2 the relative weight of the environmental objective of firm 1 and by γ 2 the corresponding relative weight of firm 2.

4.1. Necessary Conditions for a Feedback Nash Equilibrium

If ( I 1 t ) t is a solution to problems (47)–(51), then according to Theorem 4.2 in the work of Hayek (2018) [18] there exist two non-negative real numbers λ 1 and λ 2 , sequences ( μ t ) t , ( ψ t ) t , ( θ t ) t in 1 , not all nil, such that for all t, I 1 t maximizes the following Hamiltonian
H ( I 1 t , K 1 t , x t , λ 1 , λ 2 , μ t , ψ t , θ t ) = λ 1 β t a b ( K 1 t + K 2 t K 1 t I 1 t 2 2 + λ 2 ( β t x t ) ( 52 ) + μ t + 1 η x t + K 1 t + K 2 t + ψ t + 1 I 1 t + ( 1 δ ) K 1 t ( 53 ) + θ t + 1 I 2 ( K 1 t , K 2 t , x t ) + ( 1 δ ) K 2 t
We now look for a solution to the optimality conditions by assuming that
I 2 t = τ 1 2 K 1 t + τ 2 2 K 2 t + τ 3 2 .
That is, we conjecture that the firms’ decisions do not depend on the pollution stock (because pollution enters in a linear way in the firms’ objectives). This assumption implies the following first-order conditions:
0 = λ 1 β t I 1 t + ψ t + 1 ,
ψ t = λ 1 β t a 2 b K 1 t b K 2 t + μ t + 1 + ( 1 δ ) ψ t + 1 + θ t + 1 τ 1 2 ,
θ t = λ 1 β t b K 1 t + μ t + 1 + θ t + 1 τ 2 2 + 1 δ
μ t = λ 2 β t + η μ t + 1
In the system above, ψ t and θ t stand for the shadow prices of firm 1 and firm 2’s capital stocks, respectively, and μ t is the shadow price of pollution. To solve the first-order conditions we shall look for a linear investment rule such as I 1 = τ 1 1 K 1 t + τ 2 1 K 2 t + τ 3 1 . To determine the coefficients of the investment rules used by the two firms, we proceed as follows. Firstly, we solve for the dynamics of the capital stocks which is completely determined by the linear investment rules used by both firms. Secondly, we look for a sequence of shadow prices that result from the dynamics of the capital stock and pollution. Thirdly, we look for conditions on the investment rules ensuring that these rules give the optimal decisions given the sequences of shadow prices. To put it differently, we look for linear rules that satisfy a fixed-point property (for both players). This property is displayed in the following lemma.
Lemma 3.
If a feedback Nash equilibrium exists with I 1 t = τ 1 1 K 1 t + τ 2 1 K 2 t + τ 3 1 , I 2 t = τ 1 2 K 1 t + τ 2 2 K 2 t + τ 3 2 with τ 2 1 0 , τ 1 2 0 then the coefficients of these rules satisfy the following conditions
τ 1 1 = β 2 b + ( 1 δ ) τ 1 1 ( τ 1 1 + 1 δ ) + b + ( 1 δ ) τ 2 1 τ 1 2 + b z 1 2 z 1 β τ 2 2 + 1 δ 1 z 2 ( τ 1 1 + 1 δ ) ) z 2 z 1 b z 2 2 z 2 β τ 2 2 + 1 δ 1 z 1 ( τ 1 1 + 1 δ ) z 2 z 1 τ 2 1 = β 2 b + ( 1 δ ) τ 1 1 τ 2 1 + b + ( 1 δ ) τ 2 1 ( τ 2 2 + 1 δ ) b z 1 2 z 1 β τ 2 2 + 1 δ 1 τ 2 1 z 2 z 1 + b z 2 2 z 2 β τ 2 2 + 1 δ 1 τ 2 1 z 2 z 1 τ 3 1 = β a + 2 b + ( 1 δ ) τ 1 1 τ 3 1 + b + ( 1 δ ) τ 2 1 τ 3 2 λ 2 β 1 β η + ( 1 δ ) τ 3 1 + β τ 1 2 λ 2 β ( 1 β η ) β τ 2 2 + 1 δ 1 + b C f β τ 2 2 + 1 δ 1 b z 1 2 z 1 β τ 2 2 + 1 δ 1 C f ( z 2 1 ) + τ 3 1 z 2 z 1 + b z 2 2 z 2 β τ 2 2 + 1 δ 1 C f ( z 1 1 ) + τ 3 1 z 2 z 1
and an analogous system for τ 1 2 , τ 2 2 and τ 3 2 , where
C f = τ 3 1 ( τ 2 2 δ ) + τ 2 1 τ 3 2 τ 2 1 τ 1 2 + ( τ 1 1 δ ) ( τ 2 2 δ ) z 1 = τ 1 1 + τ 2 2 + 2 ( 1 δ ) ( τ 1 1 τ 2 2 ) 2 + 4 τ 2 1 τ 1 2 2 , z 2 = τ 1 1 + τ 2 2 + 2 ( 1 δ ) + ( τ 1 1 τ 2 2 ) 2 + 4 τ 2 1 τ 1 2 2 .
Proof. 
Consider the first-order conditions. Relying on a reasoning used in the proof of Proposition 1, we can assume that λ 1 = 1 . The first-order condition can then be written as
0 = β t I 1 t + ψ t + 1 , ψ t β t = a 2 b K 1 t b K 2 t λ 2 β 1 β η + ( 1 δ ) ψ t + 1 β t + θ t + 1 τ 1 2 β t , K 1 t + 1 = ψ t + 1 β t + ( 1 δ ) K 1 t θ t β t = λ 1 b K 1 t λ 2 β 1 β η + θ t + 1 β t τ 2 2 + 1 δ , K 2 t + 1 = I 2 ( K 1 t , K 2 t , x t ) + ( 1 δ ) K 2 t , = τ 1 2 K 1 t + τ 2 2 K 2 t + τ 3 2 + ( 1 δ ) K 2 t , = K 2 t ( τ 2 2 + 1 δ ) + τ 1 2 K 1 t + τ 3 2 .
Set Ψ t = ψ t β t , Θ t = θ t β t , the above system can be rewritten into:
Ψ t + 1 = I 1 t β = τ 1 1 K 1 t + τ 2 1 K 2 t + τ 3 1 β
Ψ t = a 2 b K 1 t b K 2 t λ 2 β 1 β η + ( 1 δ ) β Ψ t + 1 + β τ 1 2 Θ t + 1 ,
K 1 t + 1 = Ψ t + 1 β + ( 1 δ ) K 1 t ,
Θ t = b K 1 t λ 2 β 1 β η + Θ t + 1 β τ 2 2 + 1 δ ,
K 2 t + 1 = K 2 t ( τ 2 2 + 1 δ ) + τ 1 2 K 1 t + τ 3 2 .
Hence, from Equations (58), (60) and (62), we obtain
K 1 t + 1 = ( τ 1 1 + 1 δ ) K 1 t + τ 2 1 K 2 t + τ 3 1
K 2 t + 1 = τ 1 2 K 1 t + K 2 t ( τ 2 2 + 1 δ ) + τ 3 2
So, K 2 t = K 1 t + 1 ( τ 1 1 + 1 δ ) K 1 t τ 3 1 τ 2 1 and substituting in Equation (64), we obtain after a few algebras
K 1 t + 2 ( τ 1 1 + 1 δ + τ 2 2 + 1 δ ) K 1 t + 1 + ( τ 2 1 τ 1 2 + ( τ 1 1 + 1 δ ) ( τ 2 2 + 1 δ ) ) K 1 t
= τ 3 1 τ 3 1 ( τ 2 2 + 1 δ ) + τ 2 1 τ 3 2 .
Solving the sequential equation above, we obtain
K 1 t = c 1 z 1 t + c 2 z 2 t + C f
where
z 1 = τ 1 1 + τ 2 2 + 2 ( 1 δ ) ( τ 1 1 τ 2 2 ) 2 + 4 τ 2 1 τ 1 2 2 ,
z 2 = τ 1 1 + τ 2 2 + 2 ( 1 δ ) + ( τ 1 1 τ 2 2 ) 2 + 4 τ 2 1 τ 1 2 2 ,
C f = τ 3 1 ( τ 2 2 δ ) + τ 2 1 τ 3 2 τ 2 1 τ 1 2 + ( τ 1 1 δ ) ( τ 2 2 δ ) > 0 .
Notice that we must have
0 z 1 1 , 0 z 2 1 .
Now we can obtain K 2 t from
K 2 t = K 1 t + 1 ( τ 1 1 + 1 δ ) K 1 t τ 3 1 τ 2 1 .
In what follows, we shall obtain Θ t from (61), then Ψ t from (59) and finally use (58) to identify τ j 1 ( j = 1 , 2 , 3 ) in terms of τ j 2 . Proceeding analogously with equations from player’s 2 FOC will give τ j 2 in terms of τ j 1 .
The general solution of the sequential equation, Equation (61), is given by
Θ t = d 1 β τ 2 2 + 1 δ t + λ 2 β ( 1 β η ) ( β τ 2 2 + 1 δ 1 ) + b C f β τ 2 2 + 1 δ 1
+ b c 1 z 1 β τ 2 2 + 1 δ 1 z 1 t + b c 2 z 2 β τ 2 2 + 1 δ 1 z 2 t
where
d = Θ 0 λ 2 β ( 1 β η ) ( β τ 2 2 + 1 δ 1 ) b C f β τ 2 2 + 1 δ 1
b c 1 z 1 β τ 2 2 + 1 δ 1 b c 2 z 2 β τ 2 2 + 1 δ 1
Notice that we must have d = 0 since lim t β t Θ t = lim t θ t = 0 (as ( θ t ) t is in 1 ) and it holds that 0 < τ 2 2 + 1 δ < 1 . Indeed if τ 2 2 + 1 δ > 1 , then K 2 t is unbounded.
Now using Equation (69) in Equation (58) we obtain
Ψ t = a 2 b K 1 t b K 2 t λ 2 β 1 β η + ( 1 δ ) β τ 1 1 K 1 t + τ 2 1 K 2 t + τ 3 1 β +
β τ 1 2 λ 2 β ( 1 β η ) ( β τ 2 2 + 1 δ 1 ) +
b C f β τ 2 2 + 1 δ 1 + b c 1 z 1 β τ 2 2 + 1 δ 1 z 1 t + 1 + b c 2 z 2 β τ 2 2 + 1 δ 1 z 2 t + 1
and
Ψ t + 1 = a 2 b K 1 t + 1 b K 2 t + 1 λ 2 β 1 β η + ( 1 δ ) β τ 1 1 K 1 t + 1 + τ 2 1 K 2 t + 1 + τ 3 1 β +
β τ 1 2 λ 2 β ( 1 β η ) ( β τ 2 2 + 1 δ 1 ) +
b C f β τ 2 2 + 1 δ 1 + b c 1 z 1 β τ 2 2 + 1 δ 1 z 1 t + 2 + b c 2 z 2 β τ 2 2 + 1 δ 1 z 2 t + 2
= a + 2 b + ( 1 δ ) τ 1 1 ( τ 1 1 + 1 δ ) K 1 t + τ 2 1 K 2 t + τ 3 1 + b + ( 1 δ ) τ 2 1 τ 1 2 K 1 t + K 2 t ( τ 2 2 + 1 δ ) + τ 3 2
λ 2 β 1 β η + ( 1 δ ) τ 3 1 + β τ 1 2 λ 2 β ( 1 β η ) ( β τ 2 2 + 1 δ 1 ) +
b C f β τ 2 2 + 1 δ 1 + b z 1 2 z 1 β τ 2 2 + 1 δ 1 c 1 z 1 t + b z 2 2 z 2 β τ 2 2 + 1 δ 1 c 2 z 2 t .
Therefore,
Ψ t + 1 = 2 b + ( 1 δ ) τ 1 1 ( τ 1 1 + 1 δ ) + b + ( 1 δ ) τ 2 1 τ 1 2 K 1 t
+ 2 b + ( 1 δ ) τ 1 1 τ 2 1 + b + ( 1 δ ) τ 2 1 ( τ 2 2 + 1 δ ) K 2 t
+ a + 2 b + ( 1 δ ) τ 1 1 τ 3 1 + b + ( 1 δ ) τ 2 1 τ 3 2 λ 2 β 1 β η + ( 1 δ ) τ 3 1
+ β τ 1 2 λ 2 β ( 1 β η ) ( β τ 2 2 + 1 δ 1 ) +
b C f β τ 2 2 + 1 δ 1 + b c 1 z 1 2 z 1 β τ 2 2 + 1 δ 1 z 1 t + b c 2 z 2 2 z 2 β τ 2 2 + 1 δ 1 z 2 t .
We shall now express z 1 t and z 2 t in terms of K 1 t , K 2 t . We have
K 1 t = c 1 z 1 t + c 2 z 2 t + C f τ 2 1 K 2 t = K 1 t + 1 ( τ 1 1 + 1 δ ) K 1 t τ 3 1 , = c 1 z 1 t + 1 + c 2 z 2 t + 1 + C f ( τ 1 1 + 1 δ ) ( c 1 z 1 t + c 2 z 2 t + C f ) τ 3 1 = ( c 1 z 1 c 1 ( τ 1 1 + 1 δ ) ) z 1 t + ( c 2 z 2 c 2 ( τ 1 1 + 1 δ ) ) z 2 t ( τ 1 1 δ ) C f τ 3 1 .
We shall then solve the following system
c 1 z 1 t + c 2 z 2 t = K 1 t C f ( c 1 z 1 c 1 ( τ 1 1 + 1 δ ) ) z 1 t + ( c 2 z 2 c 2 ( τ 1 1 + 1 δ ) ) z 2 t = τ 2 1 K 2 t + ( τ 1 1 δ ) C f + τ 3 1
We obtain
z 1 t = K 1 t z 2 ( τ 1 1 + 1 δ ) C f ( z 2 1 ) τ 2 1 K 2 t τ 3 1 c 1 ( z 2 z 1 ) ,
z 2 t = τ 2 1 K 2 t + τ 3 1 K 1 t z 1 ( τ 1 1 + 1 δ ) + C f ( z 1 1 ) c 2 ( z 2 z 1 ) .
Substituting the two above expression in Equation (70), we obtain
Ψ t + 1 = 2 b + ( 1 δ ) τ 1 1 ( τ 1 1 + 1 δ ) + b + ( 1 δ ) τ 2 1 τ 1 2 + b c 1 z 1 2 z 1 β τ 2 2 + 1 δ 1 z 2 ( τ 1 1 + 1 δ ) c 1 ( z 2 z 1 )
b c 2 z 2 2 z 2 β τ 2 2 + 1 δ 1 z 1 ( τ 1 1 + 1 δ ) c 2 ( z 2 z 1 ) K 1 t
+ 2 b + ( 1 δ ) τ 1 1 τ 2 1 + b + ( 1 δ ) τ 2 1 ( τ 2 2 + 1 δ ) b c 1 z 1 2 z 1 β τ 2 2 + 1 δ 1 τ 2 1 c 1 ( z 2 z 1 )
+ b c 2 z 2 2 z 2 β τ 2 2 + 1 δ 1 τ 2 1 c 2 ( z 2 z 1 ) K 2 t
+ a + 2 b + ( 1 δ ) τ 1 1 τ 3 1 + b + ( 1 δ ) τ 2 1 τ 3 2 λ 2 β 1 β η + ( 1 δ ) τ 3 1
+ β τ 1 2 λ 2 β ( 1 β η ) ( β τ 2 2 + 1 δ 1 ) +
b C f β τ 2 2 + 1 δ 1 b c 1 z 1 2 z 1 β τ 2 2 + 1 δ 1 C f ( z 2 1 ) + τ 3 1 c 1 ( z 2 z 1 ) + b c 2 z 2 2 z 2 β τ 2 2 + 1 δ 1 C f ( z 1 1 ) + τ 3 1 c 2 ( z 2 z 1 ) .
Since Ψ t + 1 = I 1 t β = τ 1 1 K 1 t + τ 2 1 K 2 t + τ 3 1 β , using the equation above, we obtain by identification that
τ 1 1 = β 2 b + ( 1 δ ) τ 1 1 ( τ 1 1 + 1 δ ) + b + ( 1 δ ) τ 2 1 τ 1 2 + b z 1 2 z 1 β τ 2 2 + 1 δ 1 z 2 ( τ 1 1 + 1 δ ) z 2 z 1 b z 2 2 z 2 β τ 2 2 + 1 δ 1 z 1 ( τ 1 1 + 1 δ ) z 2 z 1 τ 2 1 = β 2 b + ( 1 δ ) τ 1 1 τ 2 1 + b + ( 1 δ ) τ 2 1 ( τ 2 2 + 1 δ ) b z 1 2 z 1 β τ 2 2 + 1 δ 1 τ 2 1 z 2 z 1 + b z 2 2 z 2 β τ 2 2 + 1 δ 1 τ 2 1 z 2 z 1 τ 3 1 = β a + 2 b + ( 1 δ ) τ 1 1 τ 3 1 + b + ( 1 δ ) τ 2 1 τ 3 2 λ 2 β 1 β η + ( 1 δ ) τ 3 1 + β τ 1 2 λ 2 β ( 1 β η ) β τ 2 2 + 1 δ 1 + b C f β τ 2 2 + 1 δ 1 b z 1 2 z 1 β τ 2 2 + 1 δ 1 C f ( z 2 1 ) + τ 3 1 z 2 z 1 + b z 2 2 z 2 β τ 2 2 + 1 δ 1 C f ( z 1 1 ) + τ 3 1 z 2 z 1
By analogy we can deduce the system satisfied by τ 1 2 , τ 2 2 , τ 3 2 . □
An important question is related to how the environmental appreciation parameter λ 2 affects firm 1’s incentive to invest. First of all, we observe that it does not interfere with the way the capital stocks of both firms influence I 1 . The reason for this is that the marginal pollution effect of I 1 is not influenced by K 1 or K 2 , since, whatever the value of K 1 or K 2 is, an additional unit of I 1 always enlarges pollution by one unit. However, still the marginal effect of I 1 is more negative when λ 2 is larger. This is captured by the negative effect of λ 2 on τ 3 1 , the constant part of the investment rule.

4.2. Symmetric and Partially Symmetric Nash Equilibria

From now on, following Reynolds (1987), we shall focus on partially symmetric Nash equilibria where τ 1 1 = τ 2 2 , τ 2 1 = τ 1 2 , τ 3 1 τ 3 2 . That is, firms are symmetric except with regard to the choices of the relative weight associated with the environmental objective. Indeed, for firm 1, the conditions given in the Lemma above reduce to
τ 1 1 = β 2 b + ( 1 δ ) τ 1 1 ( τ 1 1 + 1 δ ) + b + ( 1 δ ) τ 2 1 τ 2 1 + 1 2 b z 1 2 z 1 β τ 1 1 + 1 δ 1 + 1 2 b z 2 2 z 2 β τ 1 1 + 1 δ 1 τ 2 1 = β 2 b + ( 1 δ ) τ 1 1 τ 2 1 + b + ( 1 δ ) τ 2 1 ( τ 1 1 + 1 δ ) b z 1 2 z 1 β τ 1 1 + 1 δ 1 τ 2 1 z 2 z 1 + b z 2 2 z 2 β τ 1 1 + 1 δ 1 τ 2 1 z 2 z 1 τ 3 1 = β a + 2 b + ( 1 δ ) τ 1 1 τ 3 1 + b + ( 1 δ ) τ 2 1 τ 3 2 λ 2 β 1 β η + ( 1 δ ) τ 3 1 + β τ 2 1 λ 2 β ( 1 β η ) β τ 1 1 + 1 δ 1 + b C f β τ 1 1 + 1 δ 1 b z 1 2 z 1 β τ 1 1 + 1 δ 1 C f ( z 2 1 ) + τ 3 1 z 2 z 1 + b z 2 2 z 2 β τ 1 1 + 1 δ 1 C f ( z 1 1 ) + τ 3 1 z 2 z 1 .
We can obtain τ 3 2 in firm 2’s investment rule similarly (see Equation (81)). We now have the following result.
Lemma 4.
In a partially symmetric feedback Nash equilibrium, we have ( 1 δ ) τ 1 1 δ and 1 / 2 τ 2 1 < 0 .
Proof. 
By definition, in a partially symmetric feedback Nash equilibrium τ 1 1 = τ 2 2 , τ 2 1 = τ 1 2 . However, from Equations (65) and (66) and the fact that 0 z 1 1 , 0 z 1 1 , we have
z 1 = τ 1 1 | τ 2 1 | + 1 δ , 0 τ 1 1 | τ 2 1 | + 1 δ 1 ,
z 2 = τ 1 1 + | τ 2 1 | + 1 δ , 0 τ 1 1 + | τ 2 1 | + 1 δ 1 .
Now, the condition 0 τ 1 1 | τ 2 1 | + 1 δ implies | τ 2 1 | τ 1 1 + 1 δ so ( 1 δ ) τ 1 1 . Moreover, the condition τ 1 1 + | τ 2 1 | + 1 δ 1 implies τ 1 1 δ . Hence, it holds that
( 1 δ ) τ 1 1 δ .
Thus,
0 τ 1 1 + ( 1 δ ) 1
Since | τ 2 1 | = 1 2 ( z 2 z 1 ) we have | τ 2 1 | 1 2 . Let us now show that τ 2 1 < 0 .
Consider the equation
τ 2 1 = β 2 b + ( 1 δ ) τ 1 1 τ 2 1 + b + ( 1 δ ) τ 2 1 ( τ 1 1 + 1 δ ) b z 1 2 z 1 β τ 1 1 + 1 δ 1 τ 2 1 z 2 z 1 + b z 2 2 z 2 β τ 1 1 + 1 δ 1 τ 2 1 z 2 z 1
Assume τ 2 1 > 0 . The above equation can be written as
τ 2 1 = β 2 b + ( 1 δ ) τ 1 1 τ 2 1 + b + ( 1 δ ) τ 2 1 ( τ 1 1 + 1 δ ) b z 1 2 z 1 β τ 1 1 + 1 δ 1 1 2 + b z 2 2 z 2 β τ 1 1 + 1 δ 1 1 2
τ 2 1 1 β ( 1 δ ) τ 1 1 + ( 1 δ ) ( τ 1 1 + 1 δ ) =
β ( 2 b τ 2 1 b ( τ 1 1 + 1 δ ) b z 1 2 z 1 β τ 1 1 + 1 δ 1 1 2 + b z 2 2 z 2 β τ 1 1 + 1 δ 1 1 2
τ 2 1 1 β ( 1 δ ) ( 2 τ 1 1 + 1 δ ) =
β ( 2 b τ 2 1 b ( τ 1 1 + 1 δ ) b z 1 2 z 1 β τ 1 1 + 1 δ 1 1 2 + b z 2 2 z 2 β τ 1 1 + 1 δ 1 1 2
Notice that ( 1 δ ) ( 2 τ 1 1 + 1 δ ) ( 1 δ ) ( 1 + δ ) = ( 1 δ ) 2 < 1 .
Notice also that b z 2 z β τ 1 1 + 1 δ 1 is decreasing with z (knowing that τ 1 1 < 1 β ( 1 δ ) ).
So b z 2 2 z 2 β τ 1 1 + 1 δ 1 τ 2 1 z 2 z 1 b z 1 2 z 1 β τ 1 1 + 1 δ 1 τ 2 1 z 2 z 1 since z 2 z 1 .
Thus the left-hand side of the equation is positive while the right-hand side is negative, which shows that τ 2 1 cannot be positive. □
Thus, an increase in own capital stock triggers a less than proportional increase in own investment. The same remark applies for an increase in the other firm’s capital stock. More precisely, each firm’s investment rate decreases with the capital stock of the other firm. That is because the higher the capital stock of the other firm, the lower the market price (this is in line with Reynolds’ finding [10]).
It is not easy to see how the investment rule changes with b, the parameter describing how the equilibrium price depends on the production sent to the market. Figure 1 illustrates that both τ 1 1 and τ 2 1 decrease with b (the simulations are robust to change in δ and β ). The larger b , is, the more negative the effect of an increase in K 1 or K 2 on the output price. This makes investment less profitable, which is reflected by τ 1 1 and τ 2 1 being decreasing in b.
Interestingly, we see that τ 1 1 is more negative than τ 2 2 . This is because the negative effect of K 2 on I 2 is just caused by the fact that if K 2 is larger the output price is lower, which negatively affects the profitability of investment. In the work of Dawid et al. [19], this is called the size effect, and, similarly to K 2 , this effect also holds for K 1 . In addition to this effect, a larger value of K 1 also raises the cannibalization effect of investment. The cannibalization effect says that investing raises capital stock and then a decrease in the output price makes the profitability lower, because the output price decrease is multiplied by a larger capital stock. So, we conclude that the negative value of τ 2 is caused by just the size effect, while the larger negative value (in absolute terms) of τ 1 is caused by the sum of the size and the cannibalization effect.
Up to now, we have only presented necessary conditions for the existence of a feedback Nash equilibrium. We now give a sufficient condition for the existence of feedback Nash equilibria.

4.2.1. Existence and Properties of Symmetric Feedback Nash Equilibria

We first address the existence of completely symmetric equilibria when firms give the same relative weights to their environmental objective—that is, τ 3 1 = τ 3 2 , and we have
Proposition 2.
Assume that b is small enough. For any pair ( λ 2 , γ 2 ) of relative weights given by the two firms to their environmental objectives, satisfying 0 λ 2 = γ 2 a ( 1 β η ) β , there exists a symmetric feedback Nash equilibrium.
Proof. 
Consider the conditions given for a symmetric feedback Nash equilibrium:
τ 1 1 = β [ 2 b + ( 1 δ ) τ 1 1 ( τ 1 1 + 1 δ ) + b + ( 1 δ ) τ 2 1 τ 2 1 + 1 2 b z 1 2 z 1 β τ 1 1 + 1 δ 1 + 1 2 b z 2 2 z 2 β τ 1 1 + 1 δ 1 ] τ 2 1 = β [ 2 b + ( 1 δ ) τ 1 1 τ 2 1 + b + ( 1 δ ) τ 2 1 ( τ 1 1 + 1 δ ) b z 1 2 z 1 β τ 1 1 + 1 δ 1 τ 2 1 2 | τ 2 1 | + b z 2 2 z 2 β τ 1 1 + 1 δ 1 τ 2 1 2 | τ 2 1 | ] τ 3 1 = β [ a + 2 b + ( 1 δ ) τ 1 1 τ 3 1 + b + ( 1 δ ) τ 2 1 τ 3 1 λ 2 β 1 β η + ( 1 δ ) τ 3 1 + β τ 2 1 λ 2 β ( 1 β η ) β τ 1 1 + 1 δ 1 + b C f β τ 1 1 + 1 δ 1 b z 1 2 z 1 β τ 1 1 + 1 δ 1 C f ( z 2 1 ) + τ 3 1 z 2 z 1 + b z 2 2 z 2 β τ 1 1 + 1 δ 1 C f ( z 1 1 ) + τ 3 1 z 2 z 1 ]
We have τ 2 1 < 0 . So, from Equations (65) and (66), we have z 1 = τ 1 1 + τ 2 1 + ( 1 δ ) and z 2 = τ 1 1 τ 2 1 + ( 1 δ ) . Summing the first two lines of the system above, we obtain
τ 1 1 + τ 2 1 = β 3 b + ( 1 δ ) ( τ 1 1 + τ 2 1 ) ( τ 1 1 + τ 2 1 + 1 δ ) + β b z 1 2 β z 1 ( τ 1 1 + 1 δ ) 1 ,
or, using the definition of z 1
z 1 ( 1 δ ) = β 3 b + ( 1 δ ) z 1 ( 1 δ ) z 1 + β b z 1 2 β z 1 ( τ 1 1 + 1 δ ) 1 .
Set
g ( z 1 ) = z 1 ( 1 δ ) β 3 b + ( 1 δ ) z 1 ( 1 δ ) z 1 β b z 1 2 β z 1 ( τ 1 1 + 1 δ ) 1 .
Observe that g ( 0 ) = ( 1 δ ) < 0 . Let us show that g ( 1 ) > 0 . We have
g ( 1 ) = δ β 3 b + ( 1 δ ) δ β b β ( τ 1 1 + 1 δ ) 1 = 3 b β + δ ( 1 β ( 1 δ ) ) β b β ( τ 1 1 + 1 δ ) 1 > 0
since 0 τ 1 1 + ( 1 δ ) 1 and 0 < β < 1 . So, there exists z 1 [ 0 , 1 ] such that g ( z 1 ) = 0
Now subtracting the second line from the first one of the system of equations, we obtain
τ 1 1 τ 2 1 = β 2 b + ( 1 δ ) τ 1 1 ( τ 1 1 τ 2 1 + 1 δ ) + β b + ( 1 δ ) τ 2 1 ( τ 1 1 + τ 2 1 ( 1 δ ) ) + β b z 2 2 β z 2 ( τ 1 1 + 1 δ ) 1 ,
or
τ 1 1 τ 2 1 = β ( b + ( 1 δ ) ( τ 1 1 τ 2 1 ) ( τ 1 1 τ 2 1 + 1 δ ) + β b z 2 2 β z 2 ( τ 1 1 + 1 δ ) 1 ,
Using the definition of z 2 in the equation above we obtain
z 2 ( 1 δ ) = β b + ( 1 δ ) ( z 2 ( 1 δ ) ) z 2 + β b z 2 2 β z 2 ( τ 1 1 + 1 δ ) 1 .
Set
h ( z 2 ) = z 2 ( 1 δ ) β b + ( 1 δ ) ( z 2 ( 1 δ ) ) z 2 β b z 2 2 β z 2 ( τ 1 1 + 1 δ ) 1 .
We have h ( 0 ) = ( 1 δ ) < 0 . Let us now show that h ( 1 ) > 0 . Moreover,
h ( 1 ) = δ β b + ( 1 δ ) δ β b β ( τ 1 1 + 1 δ ) 1 = β b + δ 2 β b β ( τ 1 1 + 1 δ ) 1 > 0
since 0 τ 1 1 + ( 1 δ ) 1 and 0 < β < 1 . So there exists z 2 [ 0 , 1 ] such that h ( z 2 ) = 0 .
Consider now the function F defined on [ ( 1 δ ) , δ ] by
F ( τ 1 1 ) = τ 1 1 β [ 2 b + ( 1 δ ) τ 1 1 ( τ 1 1 + 1 δ ) + b + ( 1 δ ) τ 2 1 ( τ 1 1 ) τ 2 1 ( τ 1 1 ) + 1 2 b ( z 1 ( τ 1 1 ) ) 2 z 1 ( τ 1 1 ) β τ 1 1 + 1 δ 1
+ 1 2 b ( z 2 ( τ 1 1 ) ) 2 z 2 ( τ 1 1 ) β τ 1 1 + 1 δ 1 ] .
Clearly, F is continuous on [ ( 1 δ ) , δ ] . Using 0 z 1 = τ 1 1 | τ 2 1 | + 1 δ 1 and 0 z 2 = τ 1 1 + | τ 2 1 | + 1 δ 1 , we see that when τ 1 1 = ( 1 δ ) , then τ 2 1 ( ( 1 δ ) ) = z 1 ( ( 1 δ ) ) = 0 , z 2 ( ( 1 δ ) ) = 0 . Thus, F ( ( 1 δ ) ) = ( 1 δ ) < 0 .
When τ 1 1 = δ , then reasoning as above we have τ 2 1 ( δ ) = 0 , z 1 ( δ ) = 1 , z 2 ( δ ) = 1 . Then, F ( δ ) = δ β ( 2 b + ( 1 δ ) δ + b β 1 1 2 + b β 1 1 2 = δ ( 1 β ( 1 δ ) ) + 2 b β β b β 1 > 0 . Hence, there exists τ 1 1 [ ( 1 δ ) , δ ] such that F ( τ 1 1 ) = 0 . It is now important to obtain τ 2 1 as | τ 2 1 | = τ 1 1 + 1 δ z 1 .
Therefore, we have shown the existence of τ 1 1 , τ 2 1 , z 1 , z 2 that satisfy the two equations that must hold in an equilibrium.
Finally, we can use the last equation to obtain an expression for τ 3 1 ( τ 3 2 = τ 3 1 in a completely symmetric equilibrium). However, first notice using the definitions of C f , namely Equation (67) and z 1 , and noting that C ˜ f is defined in a similar way as C f , that we have in a completely symmetric equilibrium
C ˜ f = C f = τ 3 1 ( τ 1 1 δ τ 2 1 ) ( τ 1 1 δ τ 2 1 ) ( τ 1 1 δ + τ 2 1 ) = τ 3 1 ( τ 1 1 δ + τ 2 1 ) = τ 3 1 z 1 1 = τ 3 1 1 z 1 > 0 .
Building on this observation, we can use the definition of τ 3 1 to obtain
τ 3 1 = β [ a + 2 b + ( 1 δ ) τ 1 1 τ 3 1 + b + ( 1 δ ) τ 2 1 τ 3 1 λ 2 β 1 β η + ( 1 δ ) τ 3 1 + β τ 2 1 λ 2 β ( 1 β η ) β τ 1 1 + 1 δ 1 + b C f β τ 1 1 + 1 δ 1 b z 1 2 z 1 β τ 1 1 + 1 δ 1 C f ( z 2 1 ) + τ 3 1 z 2 z 1 + b z 2 2 z 2 β τ 1 1 + 1 δ 1 C f ( z 1 1 ) + τ 3 1 z 2 z 1 ] .
Notice that C f ( z 1 1 ) + τ 3 1 = 0 and C f ( z 2 1 ) + τ 3 1 z 2 z 1 = C f . Using the equation above, we then obtain
τ 3 1 = β a λ 2 β 1 β η + β τ 2 1 λ 2 β ( 1 β η ) β τ 1 1 + 1 δ 1 S + T
where
S = 1 β 2 b + ( 1 δ ) τ 1 1 + ( 1 δ ) + b + ( 1 δ ) τ 2 1
= 3 β b + 1 β ( 1 δ ) ( τ 1 1 + τ 2 1 + 1 )
T = β τ 2 1 1 z 1 1 b β τ 1 1 + 1 δ 1 b z 1 2 z 1 β τ 1 1 + 1 δ 1 .
In particular, if (as in Reynolds, 1987) λ 2 = γ 2 = 0 we have
τ 3 1 = β a S + T .
We want S + T > 0 to have τ 3 1 > 0 , otherwise, firm i’s capital stock cannot always be non-negative.
Observe that S is positive since z 1 < 1 < 1 / ( β ( 1 δ ) ) δ = 1 β δ + β δ 2 β β δ Now, let us turn to T. Recall that z 1 1 < 0 . Moreover, b β τ 1 1 + 1 δ 1 b z 1 2 z 1 β τ 1 1 + 1 δ 1 < 0 since b z 1 2 z 1 β τ 1 1 + 1 δ 1 is decreasing with z 1 (knowing that τ 1 1 < 1 β ( 1 δ ) ). (In fact, when z 1 = 0 , b β τ 1 1 + 1 δ 1 b z 1 2 z 1 β τ 1 1 + 1 δ 1 < 0 , when z 1 = 1 , this expression equals 0). So
β 2 τ 2 1 1 ( τ 1 1 δ + τ 2 1 ) ( b ( β τ 1 1 + 1 δ 1 ) b z 1 2 ( z 1 β τ 1 1 + 1 δ 1 ) ) < 0 .
Now, we notice that S + T > 0 is always true when b = 0 . So, τ 3 1 > 0 if b is small enough and λ 2 = γ 2 = 0 .
If λ 2 = γ 2 but are different from zero then sufficient conditions to have τ 3 1 > 0 under the hypotheses of the proposition are that b is small enough and 0 λ 2 a ( 1 β η ) β . □
Proposition 2 extends the result of Reynolds [10] by taking into account firms’ environmental concerns. The next proposition shows precisely how firms’ investment rules change with the common environmental concern ( λ 2 ) in the completely symmetric case.
Proposition 3.
Suppose that the assumptions of Proposition 2 hold. Then, each firm’s investment rule is a decreasing function of their environmental concern.
Proof. 
In the completely symmetric case, the constant term of both firms’ investment rules is given by Equation (74)
τ 3 1 = β a λ 2 β 1 β η + β τ 2 1 λ 2 β ( 1 β η ) β τ 1 1 + 1 δ 1 S + T .
Under the assumptions of Proposition 2, the denominator S + T is positive. We have
τ 3 1 λ 2 = β ( S + T ) ( 1 β η ) 1 + β z 2 β τ 1 1 + 1 δ 1 .
The result follows since the second fraction in the expression above is positive. □
Recall that only the constant terms in the investments rule depend on the relative weights given to the environmental objectives. The proposition states that these constant terms decrease when firms become more environmentally concerned. This proposition is somewhat expected, and it is similar to what we found in the open-loop case. It is illustrated in Figure 2.

4.2.2. Existence and Properties of Partially Symmetric Equilibria

Here, we extend the previous section by allowing for the environmental appreciation parameters to be different, i.e., λ 2 γ 2 . We obtain the following results.
Proposition 4.
The difference τ 3 1 τ 3 2 is a decreasing function of λ 2 and an increasing function of γ 2 .
Proof. 
Using the expression of τ 3 1 given in Lemma 3, we know that:
τ 3 1 = β a + 2 b + ( 1 δ ) τ 1 1 τ 3 1 + b + ( 1 δ ) τ 2 1 τ 3 2 λ 2 β 1 β η + ( 1 δ ) τ 3 1 β τ 1 2 λ 2 β ( 1 β η ) β τ 2 2 + 1 δ 1 + b C f β τ 2 2 + 1 δ 1 b z 1 2 z 1 β τ 2 2 + 1 δ 1 C f ( z 2 1 ) + τ 3 1 z 2 z 1 + b z 2 2 z 2 β τ 2 2 + 1 δ 1 C f ( z 1 1 ) + τ 3 1 z 2 z 1
However, τ 2 2 = τ 1 1 . Thus,
τ 3 1 = β a + 2 b + ( 1 δ ) τ 1 1 τ 3 1 + b + ( 1 δ ) τ 2 1 τ 3 2 λ 2 β 1 β η + ( 1 δ ) τ 3 1 β τ 1 2 λ 2 β ( 1 β η ) β τ 1 1 + 1 δ 1 + b C f β τ 1 1 + 1 δ 1 b z 1 2 z 1 β τ 1 1 + 1 δ 1 C f ( z 2 1 ) + τ 3 1 z 2 z 1 + b z 2 2 z 2 β τ 1 1 + 1 δ 1 C f ( z 1 1 ) + τ 3 1 z 2 z 1
By considering symmetry, we have (with τ 2 1 = τ 1 2 )
τ 3 2 = β a + 2 b + ( 1 δ ) τ 1 1 τ 3 2 + b + ( 1 δ ) τ 2 1 τ 3 1 γ 2 β 1 β η + ( 1 δ ) τ 3 2 β τ 1 2 γ 2 β ( 1 β η ) β τ 1 1 + 1 δ 1 + b C ˜ f β τ 1 1 + 1 δ 1 b z 1 2 z 1 β τ 1 1 + 1 δ 1 C ˜ f ( z 2 1 ) + τ 3 2 z 2 z 1 + b z 2 2 z 2 β τ 1 1 + 1 δ 1 C ˜ f ( z 1 1 ) + τ 3 2 z 2 z 1
Thus,
τ 3 1 τ 3 2 = β 2 b + ( 1 δ ) τ 1 1 ( τ 3 1 τ 3 2 ) + b + ( 1 δ ) τ 2 1 ( τ 3 2 τ 3 1 ) β ( λ 2 γ 2 ) 1 β η + ( 1 δ ) ( τ 3 1 τ 3 2 ) β τ 1 2 β ( λ 2 γ 2 ) ( 1 β η ) β τ 1 1 + 1 δ 1 + b ( C f C ˜ f ) β τ 1 1 + 1 δ 1 b z 1 2 z 1 β τ 1 1 + 1 δ 1 ( C f C ˜ f ) ( z 2 1 ) + τ 3 1 τ 3 2 z 2 z 1 + b z 2 2 z 2 β τ 1 1 + 1 δ 1 ( C f C ˜ f ) ( z 1 1 ) + τ 3 1 τ 3 2 z 2 z 1
Now, we have
C f = τ 3 1 ( τ 1 1 δ ) + τ 2 1 τ 3 2 τ 2 1 τ 2 1 + ( τ 1 1 δ ) ( τ 1 1 δ ) > 0 .
C ˜ f = τ 3 2 ( τ 1 1 δ ) + τ 2 1 τ 3 1 τ 2 1 τ 2 1 + ( τ 1 1 δ ) ( τ 1 1 δ ) > 0 .
and thus
C f C ˜ f = ( τ 3 1 τ 3 2 ) ( τ 1 1 δ ) + τ 2 1 ( τ 3 2 τ 3 1 ) ( τ 2 1 ) 2 + ( τ 1 1 δ ) 2
( τ 3 1 τ 3 2 ) ( τ 1 1 + τ 2 1 δ ) ( τ 2 1 ) 2 + ( τ 1 1 δ ) 2
= ( τ 3 1 τ 3 2 ) τ 1 1 τ 1 2 δ
We also have
( C f C ˜ f ) ( z 2 1 ) + τ 3 1 τ 3 2 = ( τ 3 1 τ 3 2 ) τ 1 1 τ 1 2 δ ( τ 1 1 τ 2 1 δ ) + τ 3 1 τ 3 2
= 0 ,
and
( C f C ˜ f ) ( z 1 1 ) + τ 3 1 τ 3 2 = ( τ 3 1 τ 3 2 ) τ 1 1 τ 1 2 δ ( τ 1 1 + τ 1 2 δ ) + τ 3 1 τ 3 2
= ( τ 3 1 τ 3 2 ) ( τ 1 1 + τ 1 2 δ ) + τ 1 1 τ 1 2 δ τ 1 1 τ 1 2 δ
= ( τ 3 1 τ 3 2 ) ( τ 1 1 + τ 1 2 δ ) + τ 1 1 τ 1 2 δ τ 1 1 τ 1 2 δ
= ( τ 3 1 τ 3 2 ) 2 τ 2 1 τ 1 1 τ 1 2 δ .
So, using z 2 z 1 = 2 τ 2 1 , we obtain
τ 3 1 τ 3 2 = β 2 b + ( 1 δ ) τ 1 1 ( τ 3 1 τ 3 2 ) + b + ( 1 δ ) τ 2 1 ( τ 3 2 τ 3 1 ) β ( λ 2 γ 2 ) 1 β η + ( 1 δ ) ( τ 3 1 τ 3 2 ) β τ 1 2 β ( λ 2 γ 2 ) ( 1 β η ) β τ 1 1 + 1 δ 1 + b β τ 1 1 + 1 δ 1 ( τ 3 1 τ 3 2 ) τ 1 1 τ 1 2 δ + b z 2 2 z 2 β τ 1 1 + 1 δ 1 ( τ 3 1 τ 3 2 ) 1 τ 1 1 τ 1 2 δ .
Hence,
( τ 3 1 τ 3 2 ) 1 β 2 b + ( 1 δ ) τ 1 1 b + ( 1 δ ) τ 2 1 + ( 1 δ ) β τ 1 2 τ 1 1 τ 1 2 δ b β τ 1 1 + 1 δ 1 b z 2 2 z 2 β τ 1 1 + 1 δ 1 = β β ( λ 2 γ 2 ) 1 β η + β τ 1 2 β ( λ 2 γ 2 ) ( 1 β η ) β τ 1 1 + 1 δ 1
Regarding the right-hand side, observe that
β ( λ 2 γ 2 ) 1 β η + β τ 1 2 β ( λ 2 γ 2 ) ( 1 β η ) β τ 1 1 + 1 δ 1 = β ( λ 2 γ 2 ) 1 β η 1 + β τ 2 1 β τ 1 1 + 1 δ 1 = β ( λ 2 γ 2 ) 1 β η β τ 2 1 ( β ( τ 1 1 + 1 δ ) 1 ) β τ 1 1 + 1 δ 1 = β ( λ 2 γ 2 ) 1 β η β ( τ 2 1 τ 1 1 1 + δ ) + 1 β τ 1 1 + 1 δ 1 = β ( λ 2 γ 2 ) 1 β η β z 2 + 1 β τ 1 1 + 1 δ 1 = β ( λ 2 γ 2 ) 1 β η 1 β z 2 1 β τ 1 1 + 1 δ
Thus,
( τ 3 1 τ 3 2 ) 1 β 2 b + ( 1 δ ) τ 1 1 b + ( 1 δ ) τ 2 1 + ( 1 δ ) β b τ 1 2 τ 1 1 τ 1 2 δ 1 β τ 1 1 + 1 δ 1 z 2 2 z 2 β τ 1 1 + 1 δ 1 = β β ( λ 2 γ 2 ) 1 β η 1 β z 2 1 β τ 1 1 + 1 δ
Notice that
1 β 2 b + ( 1 δ ) τ 1 1 b + ( 1 δ ) τ 2 1 + ( 1 δ ) = β b + 1 β ( 1 δ ) z 2 > 0 .
Moreover,
1 β τ 1 1 + 1 δ 1 z 2 2 z 2 β τ 1 1 + 1 δ 1 = ( 1 z 2 ) β z 2 ( τ 1 1 + 1 δ ) ( 1 + z 2 ) β ( τ 1 1 + 1 δ ) 1 z 2 β ( τ 1 1 + 1 δ ) 1 < 0 .
Hence,
τ 3 1 τ 3 2 = β β ( λ 2 γ 2 ) 1 β η 1 β z 2 1 β τ 1 1 + 1 δ β b + 1 β ( 1 δ ) z 2 + β 2 b τ 1 2 z 2 1 ( 1 z 2 ) β z 2 ( τ 1 1 + 1 δ ) ( 1 + z 2 ) β ( τ 1 1 + 1 δ ) 1 z 2 β ( τ 1 1 + 1 δ ) 1
= β β ( λ 2 γ 2 ) 1 β η 1 β z 2 1 β τ 1 1 + 1 δ β b + 1 β ( 1 δ ) z 2 b β 2 τ 1 2 1 + z 2 β z 2 ( τ 1 1 + 1 δ ) ( 1 + z 2 ) β ( τ 1 1 + 1 δ ) 1 z 2 β ( τ 1 1 + 1 δ ) 1 .
As for the denominator, its sign is positive if (consider the factor of b):
1 β τ 1 2 β z 2 ( τ 1 1 + 1 δ ) ( 1 + z 2 ) β ( τ 1 1 + 1 δ ) 1 z 2 β ( τ 1 1 + 1 δ ) 1 .
Observe that
β ( τ 1 1 + 1 δ ) 1 z 2 β ( τ 1 1 + 1 δ ) 1 β τ 2 1 β z 2 ( τ 1 1 + 1 δ ) ( 1 + z 2 )
= z 2 β ( τ 1 1 + 1 δ ) 1 β ( τ 1 1 + 1 δ ) 1 β τ 1 2 + β τ 2 1 ( 1 z 2 )
= z 2 β ( τ 1 1 + 1 δ ) 1 β ( τ 1 1 τ 2 1 + 1 δ ) 1 + β τ 2 1 ( 1 z 2 )
= z 2 β ( τ 1 1 + 1 δ ) 1 β z 2 1 + β τ 2 1 ( 1 z 2 )
( 1 z 2 ) 1 z 2 β ( τ 1 1 + 1 δ ) + β τ 2 1 )
( 1 z 2 ) 1 β ( τ 1 1 + 1 δ + τ 2 1 )
= ( 1 z 2 ) ( 1 β z 2 ) > 0 .
Therefore, τ 3 1 τ 3 2 decreases with λ 2 and increases with γ 2 . □
Proposition 5.
Suppose that b is small enough and that
λ 2 < a ( 1 β η ) β ( 1 + 1 δ 2 ) ,
γ 2 < a ( 1 β η ) β ( 1 + 1 δ 2 ) .
Then, there exists a partially symmetric feedback Nash equilibrium.
Proof. 
The existence of τ 1 1 and τ 2 1 (which are equal to τ 2 2 and τ 1 2 , respectively) can be obtained as in the proof of Proposition 2. We can then use the proof of Proposition 4 to obtain an expression of both τ 3 1 and τ 3 2 . To achieve this we use the fact that
C f = τ 3 1 ( τ 2 1 τ 1 1 + δ ) τ 2 1 ( τ 3 1 τ 3 2 ) ( τ 1 1 δ ) 2 ( τ 2 1 ) 2
C f ( z 2 1 ) + τ 3 1 z 2 z 1 = τ 3 1 τ 1 1 δ + τ 2 1 + 1 2 τ 3 1 τ 3 2 τ 1 1 δ + τ 2 1
= 1 2 τ 3 1 τ 3 2 1 z 1 .
C f ( z 1 1 ) + τ 3 1 z 2 z 1 = 1 2 τ 3 1 τ 3 2 τ 1 1 δ τ 2 1
= 1 2 τ 3 1 1 z 2 .
Recall that for τ 3 1 we have
τ 3 1 = β a + 2 b + ( 1 δ ) τ 1 1 τ 3 1 + b + ( 1 δ ) τ 2 1 ( τ 3 2 τ 3 1 + τ 3 1 ) λ 2 β 1 β η + ( 1 δ ) τ 3 1 β τ 1 2 λ 2 β ( 1 β η ) β τ 1 1 + 1 δ 1 + b C f β τ 1 1 + 1 δ 1 b z 1 2 z 1 β τ 1 1 + 1 δ 1 C f ( z 2 1 ) + τ 3 1 z 2 z 1 + b z 2 2 z 2 β τ 1 1 + 1 δ 1 C f ( z 1 1 ) + τ 3 1 z 2 z 1
Therefore, we obtain
τ 3 1 1 β 2 b + ( 1 δ ) τ 1 1 + b + ( 1 δ ) τ 2 1 + ( 1 δ ) = β a λ 2 β 1 β η + b + ( 1 δ ) τ 2 1 ( τ 3 2 τ 3 1 ) β τ 1 2 λ 2 β ( 1 β η ) β τ 1 1 + 1 δ 1 + b β τ 1 1 + 1 δ 1 τ 3 1 ( τ 2 1 τ 1 1 + δ ) τ 2 1 ( τ 3 1 τ 3 2 ) ( τ 1 1 δ ) 2 ( τ 2 1 ) 2 b z 1 2 z 1 β τ 1 1 + 1 δ 1 τ 3 1 τ 1 1 δ + τ 2 1 + 1 2 τ 3 1 τ 3 2 τ 1 1 δ + τ 2 1 + b z 2 2 z 2 β τ 1 1 + 1 δ 1 1 2 τ 3 1 τ 3 2 τ 1 1 δ τ 2 1
or
τ 3 1 1 β 2 b + ( 1 δ ) τ 1 1 + b + ( 1 δ ) τ 2 1 + ( 1 δ ) + β τ 2 1 b β ( τ 1 1 + 1 δ ) 1 τ 2 1 τ 1 1 + δ ( τ 1 1 δ ) 2 ( τ 2 1 ) 2 + β τ 2 1 b z 1 2 z 1 β τ 1 1 + 1 δ 1 1 τ 1 1 δ + τ 2 1 = β a λ 2 β 1 β η + b + ( 1 δ ) τ 2 1 ( τ 3 2 τ 3 1 ) + β τ 1 2 λ 2 β ( 1 β η ) β τ 1 1 + 1 δ 1 + b β τ 1 1 + 1 δ 1 τ 2 1 ( τ 3 1 τ 3 2 ) ( τ 1 1 δ ) 2 ( τ 2 1 ) 2 b z 1 2 z 1 β τ 1 1 + 1 δ 1 1 2 τ 3 1 τ 3 2 τ 1 1 δ + τ 2 1 + b z 2 2 z 2 β τ 1 1 + 1 δ 1 1 2 τ 3 1 τ 3 2 τ 1 1 δ τ 2 1
Rearranging, we obtain
τ 3 1 1 β 2 b + ( 1 δ ) τ 1 1 + b + ( 1 δ ) τ 2 1 + ( 1 δ ) β τ 2 1 b 1 β ( τ 1 1 + 1 δ ) 1 1 z 1 + β τ 2 1 b z 1 2 1 z 1 β τ 1 1 + 1 δ 1 1 z 1 = β a λ 2 β 1 β η β τ 1 2 λ 2 β ( 1 β η ) 1 β τ 1 1 + 1 δ + β ( τ 3 1 τ 3 2 ) b + ( 1 δ ) τ 2 1 1 β τ 1 1 + 1 δ 1 + b β τ 1 2 ( 1 z 1 ) ( 1 z 2 ) 1 1 β τ 1 1 + 1 δ τ 2 1 1 2 z 1 2 ( 1 z 2 ) 1 z 1 β τ 1 1 + 1 δ + 1 2 z 2 2 ( 1 z 1 ) 1 z 2 β τ 1 1 + 1 δ
Now we have
lim b 0 τ 3 1 = β a β 1 β z 2 1 β ( τ 1 1 + 1 δ ) λ 2 ( 1 δ ) τ 2 1 ( λ 2 γ 2 ) 1 β ( 1 δ ) ( τ 1 1 + τ 2 1 + 1 )
Since λ 2 is upper bounded by a ( 1 β η ) β ( 1 + 1 δ 2 ) , the expression above is positive. Then, it remains so when b is small enough. The same argument applies to τ 3 2 . □
Proposition 6.
Assume that the assumptions of Proposition 5 holds. Then, τ 3 1 is a decreasing function of λ 2 and an increasing function of γ 2 and vice versa for τ 3 2 .
Proof. 
The assertion is true for τ 3 1 when b is small enough (see Equation (121)). □
This result generalizes what was seen in the symmetric case and in the open-loop case. Notice that a larger γ 2 reduces the investment rate of firm 2 and thus its capital K 2 . This increases the output price, which in turn raises the profitability of investment for firm 1. The proposition is illustrated in Figure 3.

5. Cooperative Equilibrium

There are several ways to define cooperation for multi-objective firms. For instance, one could consider that firms collectively address the problem including all the firms’ objectives. Here, we follow a slightly different road. We suppose that both firms maximize the aggregate sum of their profits, i.e., they behave as a monopoly and minimize the aggregate pollution level. This implies that what matters for them is the effect of the aggregate value Z t of the capital stock on the aggregate profits and pollution. More formally, we assume that firms solve the following problem
max ( I i t , I 3 i t ) t t = 0 β t ( a b Z t ) Z t I i t 2 2 I 3 i t 2 2 , 0 β t x t
subject to
Z t + 1 = I i t + I 3 i t + ( 1 δ ) Z t ,
x t + 1 = η x t + Z t , Z 0 = K i 0 + K 3 i 0 being given .
and the solution is the cooperative equilibrium.
Proposition 7.
For any weight ν given to the environmental objective of the firms that is such that ν β 1 β η < a , there exists a cooperative equilibrium Z t = ( K 10 + K 20 C ) v 1 t + C where
v 1 = 2 ( 1 δ ) β ( 1 δ ) 2 + 4 b + 1 / β + Δ , Δ = ( 1 δ ) 2 + 4 b + 1 / β 2 4 ( 1 δ ) 2 / β , C = 2 ( a ν β 1 β η ) δ ( 1 + δ + 1 / β ) + 4 b .
Proof. 
Optimality conditions. Proceeding as in the preceding optimization problem, if a sequence ( I 1 t , I 2 t ) is solution to problem (122)–(124), then from Theorem 4.2 in the work of Hayek (2018) there exist λ 0 , ν 0 , ( μ t ) t in 1 , ( ψ t ) t in 1 not all nil, such that any each date, ( I 1 t , I 2 t ) maximize the following Hamiltonian
H ( I 1 t , I 2 t , Z t , x t , λ , ν , μ t , ψ t ) = λ β t ( ( a b Z t ) Z t I 1 t 2 2 I 2 t 2 2 ) ν β t x t + μ t + 1 η x t + Z t + ψ t + 1 I 1 t + I 2 t + ( 1 δ ) Z t .
The first-order conditions are given by
0 = λ β t I 1 t + ψ t + 1 ,
0 = λ β t I 2 t + ψ t + 1 ,
ψ t = λ β t a 2 b Z t + μ t + 1 + ( 1 δ ) ψ t + 1 ,
μ t = ν β t + η μ t + 1 .
It holds that λ 0 . Otherwise, we would have ψ t + 1 = 0 for all t, which in turn would yield μ t + 1 = 0 for all t, and thus ν = 0 , so that all those multipliers would be nil, which is impossible. Then set λ = 1 . Notice that we obtain I 1 t = I 2 t = ψ t + 1 β t . This is because firms’ cost functions are similar. Now by a standard argument, we also obtain
μ t = ν β t 1 β η .
Let us now solve the following dynamic system
ψ t = β t a 2 b Z t + ( ν β t + 1 1 β η ) + ( 1 δ ) ψ t + 1 ,
Z t + 1 = 2 ψ t + 1 β t + ( 1 δ ) Z t ,
x t + 1 = η x t + Z t .
Setting ψ t β t = Ψ t in the system above, we obtain
Ψ t = a 2 b Z t ν β 1 β η + ( 1 δ ) β Ψ t + 1 ,
Z t + 1 = 2 β Ψ t + 1 + ( 1 δ ) Z t .
Using Equation (133) in (134), we obtain
( 1 δ ) Z t + 2 Z t + 1 ( 1 δ ) 2 + 4 b + 1 / β + ( 1 δ ) Z t / β = 2 ν β 1 β η 2 a
We have
Z t = w 1 v 1 t + w 2 v 2 t + C
where
v 1 = ( 1 δ ) 2 + 4 b + 1 / β Δ 2 ( 1 δ ) ,
v 2 = ( 1 δ ) 2 + 4 b + 1 / β + Δ 2 ( 1 δ ) ,
Δ = ( 1 δ ) 2 + 4 b + 1 / β 2 4 ( 1 δ ) 2 / β ,
C = 2 ( a ν β 1 β η ) δ ( 1 + δ + 1 / β ) + 4 b .
Notice that v 1 = 2 ( 1 δ ) β [ ( 1 δ ) 2 + 4 b + 1 / β + Δ ] < 1 , and as v 1 v 2 = 1 / β > 1 , v 2 > 1 . Since Z t is bounded, we have w 2 = 0 and thus Z t = ( K 10 + K 20 C ) v 1 t + C . □
This equilibrium is built in the same way as the open-loop equilibria and has similar properties. We study the differences between the different equilibria in the next section.

6. Comparison

Some comparisons of equilibria seem to only be amenable by way of a numerical analysis. We shall, however, only focus on analytical results. We begin by studying the differences between open-loop and feedback equilibria and then we consider the differences between cooperative and non-cooperative open-loop equilibria.

6.1. Comparing the Open-Loop and the Feedback Nash Analysis in the Completely Symmetric Case

Let us fist contrast the open-loop and the feedback cases when we have full symmetry and when firms are not environmentally concerned.
Case λ 2 = γ 2 = 0 (no environmental concern)
The long-run equilibrium values of the capital stock in the open-loop and feedback Nash equilibria are:
D 1 = a δ ( 1 δ 1 / β ) 3 b = a δ ( 1 / β ( 1 δ ) ) + 3 b
C f = a ( 1 z 1 ) ( 1 / β ( 1 δ ) ( z 1 + δ ) ) + 3 b + T / β )
where T is a function that goes to zero when so does b. We have the following result
Proposition 8.
Suppose that τ 1 1 < 0 . Then there is a threshold b ̲ such that whenever b is lower than b ̲ , the long-run value C f of the (total) capital stock in the feedback Nash equilibrium is lower than the long-run value D 1 of the capital stock in the open-loop equilibrium. When τ 1 1 > 0 , C f can be either greater of lower than D 1 1 .
Proof. 
First of all we rearrange the expression of the long-run value of the capital stock in a feedback Nash equilibrium. From Equation (73) and the assumption that firms are symmetric, we obtain
C f = τ 3 1 τ 1 1 + τ 2 1 δ = τ 3 1 1 z 1 .
Likewise, from Equation (78) we have τ 3 1 = β a S + T and from Equation (77) we obtain
S + T = 1 β ( 3 b + ( 1 δ ) ( τ 1 1 + τ 2 1 + 1 ) )
+ β 2 τ 2 1 b ( τ 1 1 δ + τ 2 1 ) ( 1 ( β τ 1 1 + 1 δ 1 ) z 1 2 ( z 1 β τ 1 1 + 1 δ 1 ) )
= 1 β ( 3 b + ( 1 δ ) ( z 1 + δ ) ) + T
where T goes to zero when b goes to zero.
Let us study the sign of the difference d ˜ between the denominators of C f and D 1 , where
d ˜ = ( 1 z 1 ) ( 1 / β ( 1 δ ) ( z 1 + δ ) + 3 b + T / β ) δ ( 1 / β ( 1 δ ) ) 3 b = ( 1 z 1 δ ) 1 / β + ( 1 δ ) ( δ ( 1 z 1 ) ( z 1 + δ ) ) + ( 1 z 1 ) T / β 3 b z 1 = ( 1 z 1 δ ) 1 / β + ( 1 δ ) ( δ z 1 δ + z 1 2 + z 1 δ ) + ( 1 z 1 ) T / β 3 b z 1 = ( 1 z 1 δ ) 1 / β ( 1 δ ) z 1 ( 1 z 1 δ ) + ( 1 z 1 ) T / β 3 b z 1 = ( 1 z 1 δ ) ( 1 / β ( 1 δ ) z 1 ) + ( 1 z 1 ) T / β 3 b z 1 .
To study the sign of
( 1 z 1 δ ) ( 1 / β ( 1 δ ) z 1 )
notice that z 1 = τ 1 1 + τ 2 1 + ( 1 δ ) , 1 z 1 δ = τ 1 1 τ 2 1 . So, when τ 1 1 < 0 we have ( 1 z 1 δ ) > 0 and hence
( 1 z 1 δ ) ( 1 / β ( 1 δ ) z 1 ) > 0
When b goes to zero, T goes to zero and the difference d ˜ > 0 ; thus, C f < D 1 .
However the sign is indeterminate when τ 1 1 > 0 . If τ 1 1 > | τ 2 1 | then d ˜ < 0 and thus C f > D 1 . If τ 1 1 < | τ 2 1 | and b goes to zero, T goes to zero and the difference d ˜ becomes positive. Therefore, C f < D 1 . □
The important conclusion from Proposition 8 is that scenarios exist where firms accumulate more capital under open loop than under feedback. This contradicts Reynolds [10] who finds that firms invest more under a feedback information structure. The difference between the two approaches is that Reynolds studies a continuous-time framework, whereas we work in discrete time. Apparently, choosing either of the two can generate substantially different results.
Case λ 2 = γ 2 0 (environmental concern)
Next, we determine the effect of a marginal increase in λ 2 on the steady-state value of the capital stock associated with the open-loop and feedback Nash equilibria.
Proposition 9.
In the long-run, the capital stock in a feedback Nash equilibrium is less sensitive to a change in λ 2 than in an open-loop equilibrium. Moreover a rise in environmental concern leads to a smaller decrease in the capital stock in the feedback Nash equilibrium than in the open-loop equilibrium.
Proof. 
As for the open-loop equilibrium, recall from (41) that
D 1 λ 2 = 2 b δ ( 1 δ 1 / β ) b δ ( 1 δ 1 / β ) β 1 β η 3 b δ ( 1 δ 1 / β )
As for the feedback Nash equilibrium, we obtain from (73) and (74)
C f λ 2 = 1 1 z 1 β [ β 1 β η + β τ 2 1 β ( 1 β η ) ( β τ 1 1 + 1 δ 1 ) ] S + T
with
S + T = 1 β ( 3 b + ( 1 δ ) ( τ 1 1 + τ 2 1 + 1 ) )
+ β 2 τ 2 1 b ( τ 1 1 δ + τ 2 1 ) ( 1 ( β τ 1 1 + 1 δ 1 ) z 1 2 ( z 1 β τ 1 1 + 1 δ 1 ) ) .
From what was seen in Proposition 8 when either τ 1 < 0 or τ 1 1 < τ 2 1 , and b is small enough, then ( 1 z 1 ) ( S + T ) 3 b δ ( 1 δ 1 / β ) .
Moreover, regarding the numerator of C f λ 2 , we have
1 + β τ 2 1 1 β ( τ 1 1 + 1 δ ) 1 = 1 + β ( τ 2 1 τ 1 1 ( 1 δ ) ) β ( τ 1 1 + 1 δ ) 1 .
But since 0 < τ 1 1 τ 2 1 + 1 δ < 1 , we have 0 > τ 1 1 + τ 2 1 ( 1 δ ) > 1 , and this implies that the numerator of C f λ 2 is negative. Then, it follows that
C f λ 2 < D 1 λ 2 .
From Proposition 9 we can draw the important conclusion that in a feedback Nash equilibrium capital stock is less sensitive to environmental appreciation than in the open-loop case.

6.2. Cooperative vs. Non-Cooperative Open-Loop Equilibria

Here, we ask whether pollution could be lower in an open-loop non-cooperative equilibrium, compared to the cooperative solution. To address this question we focus on the long run. Denote by Z n c the long-run total value of the capital stock when firms compete and Z c the total value of the capital stock when firms cooperate. Notice that to compare the long-run values of pollution under non-cooperation and cooperation, it is sufficient to compare the total long-run value of the capital stock, as pollution is a positive linear function of production, which is itself equal to the capital stock. After a little algebra, we obtain the following proposition.
Proposition 10.
The capital stock in a duopoly is lower than in the cooperative equilibrium whenever
2 a ( 1 β η ) b β + 2 ν δ ( 1 + δ + 1 β ) + 3 b δ ( 1 + δ + 1 β ) + 4 b < λ 2 + γ 2 .
Proof. 
The proposition follows from the comparison of the long-run value C of the capital stock in the open-loop Nash equilibria given by Equation (26) and the long-run value C of the capital stock in the cooperative equilibrium given by Equation (139). □
Recall, however, that, according to (38), λ 2 + γ 2 < 2 a ( 1 β η ) β . Thus, when λ 2 + γ 2 is large, but not too large, it is possible that firms pollute less when they compete than when they cooperate. To obtain this result it is necessary that ν < ( λ 2 + γ 2 ) / 2 .
What is the intuition for this result? Recall that firms do not coordinate their environmental concerns. So, they can be more concerned than in the cooperative case. Due to this, they can reduce their production more than in the cooperative case. This result extends the findings of Crettez and Hayek [9] to the case where firms accumulate capital.

7. Conclusions

This paper studies a dynamic duopoly in which firms accumulate capital and have environmental concerns. They both take profits and pollution explicitly into account. They thus have multiple objectives whose weights are not a priori fixed. We have studied open-loop and feedback Nash equilibria. In each kind of equilibrium, each firm chooses a Pareto-optimal solution to its multi-objective problem.
We have paid attention to the existence of completely or partially symmetric Nash equilibria and to what we have called the cooperative equilibrium. We have obtained the new result that, compared to an open-loop information structure, in a feedback Nash equilibrium the firms are less sensitive to changes in environmental appreciation. Moreover, whereas it is known from the continuous-time differential game literature that firms invest more in a feedback information structure compared to an open-loop one, we detect scenarios where the opposite holds (in the case where firms have no environmental concerns). The complexity of the equations made it difficult to analytically compare all the different cases. The fully asymmetric equilibria were not treated either for the same reason.
We have also shown that firms may over-reduce pollution in a non-cooperative equilibrium compared to the cooperative equilibrium. This result extends what was found in a setting overlooking capital accumulation to the case where capital accumulation is possible.
There are at least four issues that need further research. In this paper, we have assumed that the only way firms can reduce pollution is to slow or decrease capital accumulation. There is no denying that this is a strong assumption. It is, however, relevant is some cases where companies choose to close an airline for short hauls or a coal-fire power station. Yet, an avenue for further research is to study how firms would change their capital accumulation policies if they could mitigate pollution or rely on more efficient technologies.
Second, the demand-side of our model does not depend on the firms’ concerns for the environment. Yet, it is plausible that there exists a positive link between environmental consciousness of firms and the demand for their products (because consumers better value the products sold by environmentally friendly firms). Therefore, it might be worthwhile for firms to take this positive link into account when deciding on their investment programs. Third, one could study some sequential move games, like a Stackelberg duopoly, and compare the different equilibria. A fourth avenue for further research would be the consideration of the scenario in which the future values of pollution are discounted at different rates (see, for instance, (Cabo et al. [20])).

Author Contributions

Conceptualization; methodology; formal analysis; writing—original draft preparation; writing—review and editing: B.C.; N.H. and P.M.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wirl, F.; Feichtinger, G.; Kort, P.M. Individual Firm and Market Dynamics of CSR Activities. J. Econ. Behav. Organ. 2013, 86, 169–182. [Google Scholar] [CrossRef]
  2. Lambertini, L.; Palestini, A.; Tampieri, A. CSR in an Asymmetric Duopoly with Environmental Externality. South. Econ. J. 2016, 83, 236–252. [Google Scholar] [CrossRef]
  3. Yanase, A. Corporate Environmentalism in Dynamic Oligopoly. Strateg. Behav. Environ. 2013, 3, 223–250. [Google Scholar] [CrossRef]
  4. Feichtinger, G.; Lambertini, L.; Leitmann, G.; Wrzaczek, S. R & D for Green Technologies in a Dynamic Oligopoly: Schumpeter, Arrow and Inverted U’s. Eur. J. Oper. Res. 2016, 249, 1131–1138. [Google Scholar]
  5. Rettieva, A. Equilibria in Dynamic Multicriteria Games. Int. Game Theory Rev. 2017, 19, 1750002. [Google Scholar] [CrossRef]
  6. Rettieva, A. Dynamic Multicriteria Games with Finite Horizon. Mathematics 2018, 6, 156. [Google Scholar] [CrossRef] [Green Version]
  7. Kuzyutin, D.; Gromova, E.; Pankratova, Y. Sustainable Cooperation in Multicriteria Multistage Games. Oper. Res. Lett. 2018, 46, 557–562. [Google Scholar] [CrossRef]
  8. Kuzyutin, D.; Smirnova, N.; Gromova, E. Long-Term Implementation of the Cooperative Solution in a Multistage Multicriteria Game. Oper. Res. Perspect. 2019, 6, 100107. [Google Scholar] [CrossRef]
  9. Crettez, B.; Hayek, N. A Dynamic Multi-objective Duopoly Game with Environmentally Concerned Firms. Int. Game Theory Rev. 2021, 2150008. [Google Scholar] [CrossRef]
  10. Reynolds, S.S. Investment, Preemption and Commitment in an Infinite Horizon Model. Int. Econ. Rev. 1987, 28, 69–88. [Google Scholar] [CrossRef]
  11. Bhaumik, A.; Roy, S.K.; Li, D.-F. Analysis of Triangular Intuitionistic Fuzzy Matrix Games Using Robust Ranking. J. Intell. Fuzzy Syst. 2017, 33, 327–336. [Google Scholar] [CrossRef]
  12. Roy, S.K.; Bhaumik, A. Intelligent Water Management: A Triangular Type-2 Intuitionistic Fuzzy Matrix Games Approach. Water Resour Manag. 2018, 32, 949–968. [Google Scholar] [CrossRef]
  13. Bhaumik, A.; Roy, S.K.; Weber, G.W. Hesitant interval-valued intuitionistic fuzzy-linguistic term set approach in Prisoners’ dilemma game theory using TOPSIS: A case study on Human-trafficking. Cent. Eur. J. Oper. Res. 2020, 28, 797–816. [Google Scholar] [CrossRef]
  14. Bhaumik, A.; Roy, S.K.; Li, D.F. (α, β, γ)-cut set based ranking approach to solving bi-matrix games in neutrosophic environment. Soft Comput. 2021, 25, 2729–2739. [Google Scholar] [CrossRef]
  15. Roy, S.K.; Maiti, S.K. Reduction methods of type-2 fuzzy variables and their applications to Stackelberg game. Appl. Intell. 2020, 50, 1398–1415. [Google Scholar] [CrossRef]
  16. Lambertini, L. Differential Games in Industrial Economics; Cambridge University Press: Cambridge, UK, 2018. [Google Scholar]
  17. Qu, S.; Ji, Y. The Worst-Case Weighted Multi-Objective Game with an Application to Supply Chain Competitions. PLoS ONE 2016, 11, e0147341. [Google Scholar] [CrossRef] [PubMed]
  18. Hayek, N. Infinite-Horizon Multiobjective Optimal Control Problems for Bounded Processes. Discret. Contin. Dyn. Syst. Ser. S 2018, 11, 1121. [Google Scholar] [CrossRef] [Green Version]
  19. Dawid, H.; Kopel, M.; Kort, P.M. New Product Introduction and Capacity Investment by Incumbents: Effects of Size on Strategy. Eur. J. Oper. Res. 2013, 230, 133–142. [Google Scholar] [CrossRef]
  20. Cabo, F.; Martín-Herràn, G.; Martínez-García, M.P. Non-constant Discounting, Social Welfare and Endogenous Growth with Pollution Externalities. Environ. Resour. Econ. 2020, 76, 369–403. [Google Scholar] [CrossRef]
Figure 1. Effects of change of b on τ 1 1 and τ 2 1 , where tau(1) and tau(2) denote, respectively, τ 1 1 and τ 2 1 .
Figure 1. Effects of change of b on τ 1 1 and τ 2 1 , where tau(1) and tau(2) denote, respectively, τ 1 1 and τ 2 1 .
Mathematics 09 01983 g001
Figure 2. τ 3 1 = τ 3 2 depending on λ 2 = γ 2 , where tau3 denotes τ 3 1 = τ 3 2 .
Figure 2. τ 3 1 = τ 3 2 depending on λ 2 = γ 2 , where tau3 denotes τ 3 1 = τ 3 2 .
Mathematics 09 01983 g002
Figure 3. τ 3 1 and τ 3 2 depending on λ 2 and γ 2 , where tau31 and tau32 denote, respectively, τ 3 1 and τ 3 2 .
Figure 3. τ 3 1 and τ 3 2 depending on λ 2 and γ 2 , where tau31 and tau32 denote, respectively, τ 3 1 and τ 3 2 .
Mathematics 09 01983 g003
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Crettez, B.; Hayek, N.; Kort, P.M. A Dynamic Multi-Objective Duopoly Game with Capital Accumulation and Pollution. Mathematics 2021, 9, 1983. https://doi.org/10.3390/math9161983

AMA Style

Crettez B, Hayek N, Kort PM. A Dynamic Multi-Objective Duopoly Game with Capital Accumulation and Pollution. Mathematics. 2021; 9(16):1983. https://doi.org/10.3390/math9161983

Chicago/Turabian Style

Crettez, Bertrand, Naila Hayek, and Peter M. Kort. 2021. "A Dynamic Multi-Objective Duopoly Game with Capital Accumulation and Pollution" Mathematics 9, no. 16: 1983. https://doi.org/10.3390/math9161983

APA Style

Crettez, B., Hayek, N., & Kort, P. M. (2021). A Dynamic Multi-Objective Duopoly Game with Capital Accumulation and Pollution. Mathematics, 9(16), 1983. https://doi.org/10.3390/math9161983

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop