Next Article in Journal
Advanced Multimodal Sentiment Analysis with Enhanced Contextual Fusion and Robustness (AMSA-ECFR): Symmetry in Feature Integration and Data Alignment
Previous Article in Journal
An Algorithmic Evaluation of a Family of Logarithmic Integrals and Associated Euler Sums
Previous Article in Special Issue
Higher-Order Delay Differential Equation with Distributed Deviating Arguments: Improving Monotonic Properties of Kneser Solutions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Oscillatory Features of Fourth-Order Emden–Fowler Differential Equations with Sublinear Neutral Terms

1
Faculty of Science, Mansoura University, Mansoura 35516, Egypt
2
Department of Mathematical Sciences, College of Science, Princess Nourah bint Abdulrahman University, P.O. Box 84428, Riyadh 11671, Saudi Arabia
3
Department of Mathematics, College of Science, Qassim University, P.O. Box 6644, Buraydah 51452, Saudi Arabia
4
Department of Mathematics, Faculty of Science, Mansoura University, Mansoura 35516, Egypt
*
Authors to whom correspondence should be addressed.
Symmetry 2024, 16(7), 933; https://doi.org/10.3390/sym16070933
Submission received: 13 June 2024 / Revised: 12 July 2024 / Accepted: 19 July 2024 / Published: 22 July 2024

Abstract

:
This article examines the oscillatory characteristics of a fourth-order Emden–Fowler differential equation, specifically when it includes a sublinear neutral term. Our methodology centers on establishing multiple theorems that introduce innovative conditions to guarantee that there are no positive solutions to the examined equation. Due to the symmetry between non-oscillatory solutions, we obtain oscillation conditions by excluding only positive solutions. We employ the Riccati technique in various ways to achieve this objective. The criteria presented in this study complement and generalize many findings published in the literature. We support the efficiency of our findings by applying them to an example.

1. Introduction

In this study, we aim to investigate the fourth-order Emden–Fowler delay differential equations (DDE) with a sublinear neutral term, represented by
a ( ν ) x ν + g ν x α ξ ν γ + p ( ν ) x β ( λ ( ν ) ) = 0 ,
where ν ν 0 . We assume the following throughout this paper:
(H1)
0 < α 1 ,   β and γ are ratios of two positive odd integers;
(H2)
ξ ,   λ C 1 ( [ ν 0 , ) , R ) ,   ξ ( ν ) ν ,   λ ν ν ,   λ ( ν ) > 0 , and lim ν ξ ν = lim ν λ ν   = ;
(H3)
a C 1 ( [ ν 0 , ) , 0 , ) ,   a ν 0 ,   g ,   p C ( [ ν 0 , ) , 0 , ) ,   lim ν g ν = 0 , and
ν 0 ν 1 a 1 / γ s d s as ν .
A solution of (1) is defined as a function x C ( [ ν x , ) ) , ν x ν 0 , which satisfies x + g · x α ξ C 3 [ ν x , ) , a · x + g · x α ξ γ C 1 [ ν x , ) , and satisfies (1) on [ ν x , ) . Only solutions x of (1) that meet the condition
sup { | x ( ν ) | : ν ν * } > 0 , for all ν * ν x
are taken into consideration.
A solution of (1) is said to be non-oscillatory if it is eventually positive or eventually negative; otherwise, it is said to be oscillatory. If every solution to Equation (1) is oscillatory, then the equation is considered oscillatory.
Differential equations (DEs) have long been a milestone of mathematics and its application to the sciences. These equations serve as incredible tools for modeling a wide range of natural phenomena, ranging from population growth to the behavior of electrical circuits, disease spread, and celestial movement. Their importance is appraised by the fact that they offer a systematic approach to understanding and predicting change, making them indispensable in physics, engineering, biology, and numerous other disciplines; see [1,2,3,4,5].
The qualitative theory of differential inequalities focuses on the study of solution characteristics such as stability, periodicity, symmetry, and oscillation. Research into oscillatory phenomena in differential equations is a crucial and captivating area in both mathematics and applied sciences. Oscillations, characterized by periodic and repetitive patterns, are essential for understanding dynamic processes in various natural and engineered systems. These patterns appear in fields like physics, engineering, biology, and economics, offering valuable insights into system behaviors over time. This study delves into the intricate realm of oscillations in differential equations, aiming to uncover their fundamental principles, explore their applications, and contribute to the growing body of knowledge in this intriguing field (see [6,7,8,9,10]).
Oscillation criteria play a pivotal role in the study of nonlinear dynamics by providing sufficient criteria for the oscillatory features of solutions to DEs. In various applications, the ability to predict whether a system’s solutions will oscillate or remain steady is crucial. For instance, in mechanical and structural engineering, understanding oscillation helps in designing structures that can withstand periodic forces without leading to resonance, which could result in structural failure. In biological systems, oscillation criteria are essential for modeling phenomena such as cardiac rhythms or neural activity, where oscillatory patterns are indicative of healthy function or pathological conditions. In the realm of electrical engineering, oscillation criteria are fundamental for the design of oscillators and filters in circuits, ensuring that they perform their intended functions efficiently. Furthermore, in control systems engineering, these criteria aid in the design of controllers that manage the dynamic behavior of systems, preventing unwanted oscillations that can cause instability. Thus, oscillation criteria are indispensable tools in nonlinear dynamics, providing insights that help in the prediction, analysis, and control of complex systems across a wide range of scientific and engineering disciplines; see [11,12].
The Emden–Fowler equations, named after Robert Emden and Arthur Fowler, are central to the study of fourth-order differential models, particularly in mathematical physics. These equations are essential for modeling the structure and behavior of stellar objects such as white dwarfs, as they describe the complex relationships between pressure, density, and temperature. This modeling significantly contributes to the understanding of physical processes inside stars and their evolution. The applications of the Emden–Fowler equations extend to various fields. In thermodynamics, these equations describe temperature equilibria in complex thermal systems such as furnaces or geothermal models. In biological sciences and mathematical biology, they are used to model population dynamics and the spread of diseases, such as bacterial growth or virus transmission in different environments. In mechanical engineering and applied physics, they aid in modeling physical processes like fluid flow through porous media or material behavior under specific conditions. In geology, they assist in studying subsurface fluid movement and thermal distribution within the Earth’s crust, aiding in understanding geological formations and optimizing resource extraction; see [13,14,15,16].
The prior in-depth literature has enhanced our understanding of oscillatory properties of Emden–Fowler differential equations, with numerous influential authors shaping this area of study. Their significant contributions have left a lasting impact on the field, continuing to influence its direction and outcomes and attracting prominent attention. Below, we highlight some of these significant contributions:
Agarwal et al. [17] explored the oscillatory characteristics of DEs with a neutral term, described by
a ( ν ) x ν + g ν x α ξ ν + p ( ν ) x ( λ ( ν ) ) = 0 .
They established several sufficient conditions for the oscillations of these equations, considering the cases
ν 0 1 a s d s = ,
and
ν 0 1 a s d s < .
Baculíková et al. [18] examined the oscillatory properties of neutral differential equations (NDEs), represented by
a ( ν ) x ν g ν x α ξ ν + p ( ν ) x β ( λ ( ν ) ) = 0 .
They introduced new oscillatory criteria under the condition (3).
Zhang et al. [19] studied the oscillatory properties of a specific class of NDEs, given by
a ν κ ν α 1 κ ν + p ν x λ ν α 1 x λ ν = 0 ,
where κ ν = x ν + i = 1 m g i ν x ξ i ν . Their results simplify the analysis of these equations.
Tamilvanan et al. [20] investigated the oscillatory characteristics of the Emden–Fowler DEs, expressed as
a ( ν ) x ν + g ν x α ξ ν + p ( ν ) x β ( λ ( ν ) ) = 0 .
Wu et al. [21] investigated the oscillatory characteristics of Emden–Fowler DDEs, described by
a ( ν ) x ν + g ν x α ξ ν γ + p ( ν ) x β ( λ ( ν ) ) = 0 .
They introduced new criteria for determining oscillatory behavior under the conditions
ν 0 1 a 1 / γ s d s = ,
and
ν 0 1 a 1 / γ s d s < .
El-Nabulsi et al. [22] studied the oscillation of solutions to DEs characterized by
a ν x ν α + p ( ν ) f x ( λ ( ν ) ) = 0 ,
where f y / y α k > 0 for y 0 , and condition (2) holds. Related studies by Zhang et al. [10] and Moaaz et al. [23] further explored the oscillatory properties of (5) under condition (3).
Li and Rogovchenko [24] investigated the asymptotic behavior of solutions to a specific class of higher-order Emden–Fowler DDEs of the form
a ν x n ν α + p ( ν ) x β ( λ ( ν ) ) = 0 .
They introduced new theorems that refined and advanced the understanding of these equations.
Graef et al. [25] introduced a new technique for linearizing even-order nonlinear NDEs of the form
x ν + g ν x α ξ ν n + p ( ν ) x β ( λ ( ν ) ) = 0 , n 2 .
They developed new oscillation criteria by comparing these equations with first-order linear DDEs.
Baculikova et al. [26] and Xing [27] established oscillation criteria for a specific class of higher-order quasi-linear NDEs:
a ( ν ) x ν + g ν x ξ ν n 1 γ + p ( ν ) x γ ( λ ( ν ) ) = 0 ,
subject to condition (2). Alnafisah et al. [28] introduced refined properties that optimize the oscillatory characteristics of solutions to Equation (6) under condition (3).
Our research aims to extend the investigation of oscillatory behavior to a specific class of fourth-order NDEs, particularly those with sublinear neutral terms. The motivation of our work is to extend and generalize previous studies (for example, [21]) that explored the asymptotic properties and oscillatory nature of second-order NDEs.

2. Preliminaries

This section introduces some important preliminary results related to the monotonic properties of non-oscillatory solutions to the equation under study, which contribute significantly to obtaining our results. We start with the following notations:
κ ν : = x ν + g ν x α ξ ν ,
g ˜ ν : = 1 g ν c 1 1 α ,
and
p ˜ ν : = p ( ν ) g ˜ β λ ( ν ) .
Lemma 1
([29]). Suppose that ϱ C n [ ν 0 , ) , R + , where ϱ ( n ) ( ν ) has a constant sign and is non-zero on [ ν 0 , ) . Additionally, suppose there is ν 1 ν 0 such that ϱ ( n 1 ) ( ν ) ϱ ( n ) ( ν ) 0 for every ν ν 1 . If lim ν ϱ ( ν ) 0 , then for any δ ( 0 , 1 ) , there is ν ϵ [ ν 1 , ) such that
ϱ ( ν ) ϵ ( n 1 ) ! ν n 1 ϱ ( n 1 ) ( ν ) ,
for ϱ [ ν ϵ , ) .
Lemma 2
([30]). Let γ be a ratio of two odd positive integers and A and B be constants. Then,
B u A u γ + 1 / γ γ γ γ + 1 γ + 1 B γ + 1 A γ , a > 0 .
Lemma 3
([31]). Let ϱ C n [ ν 0 , ) , 0 , , ϱ i ν > 0 for i = 1 , 2 , , n , and ϱ n + 1 ν 0 , eventually. Then, eventually,  ϱ ν / ϱ ν ϵ ν / n for every ϵ 0 , 1 .
Lemma 4
([32]). Suppose that x represents an eventually positive solution of (1). In such a case, x will eventually fulfill the conditions of the following scenarios:
C 1 : κ > 0 , κ > 0 , κ > 0 , κ > 0 , a κ γ < 0 , C 2 : κ > 0 , κ > 0 , κ < 0 , κ > 0 , a κ γ < 0 ,
for ν ν 1 ν 0 .
Notation 1.
The category Ω i denotes the set of all solutions that eventually become positive, with the corresponding function satisfying C i for i = 1 , 2 . We note that all eventually positive solutions belong to Ω 1 Ω 2 .

3. Auxiliary Lemmas

In the following discussion, we will establish several novel conditions that are adequate for excluding positive solutions that satisfy Equation (1). We will then combine these conditions to derive oscillation criteria.
Lemma 5.
Let us consider x Ω 1 Ω 2 . Then, we have
x ν > g ˜ ν κ ν ,
and
a ( ν ) κ ( ν ) γ p ˜ ν κ β λ ( ν ) .
Proof. 
From (7), we have κ ν x ν . Using Lemma 4 for ν 1 ν 0 , we obtain
a ( ν ) κ ( ν ) γ < 0 .
Since λ ν ν , we can derive from (11) that
a ( ν ) κ ( ν ) γ a ( λ ν ) κ ( λ ν ) γ , ν ν 1 .
Notably, κ ( ν ) > 0 . Therefore, there exists c 1 > 0 such that κ ν c 1 for sufficiently large ν . By using (7), we conclude that
x v = κ v g v x α ξ v κ v g v κ α ξ v κ v g v κ α v = 1 g v κ α 1 v κ v .
Given that 0 < α 1 , and κ ν c 1 > 0 , we have
κ α 1 ν c 1 α 1 = 1 c 1 1 α .
Substituting this inequality into (13), we obtain
x ν 1 g ν c 1 1 α κ ν = g ˜ ν κ ν .
Finally, based on (1), we can establish
a ( ν ) κ ν γ = p ( ν ) x β ( λ ( ν ) ) p ( ν ) 1 g λ ( ν ) c 1 1 α β κ β λ ( ν ) = p ˜ ν κ β λ ( ν ) .
 □
Lemma 6.
Assume that (2) holds. If β γ and there is a nondecreasing function μ C 1 ( [ ν 0 , ) ,   0 , ) such that
lim sup ν ν 0 ν μ s p ˜ s 2 γ a λ s μ s γ + 1 γ + 1 γ + 1 c 2 β γ ϵ λ 2 s λ s μ s γ d s = ,
holds for every c 2 > 0 ,   ϵ 0 , 1 , then Ω 1 = .
Proof. 
We assume, for the sake of contradiction, that x ν Ω 1 . Now, define the function w ν by
w ν : = μ ν a ( ν ) κ ν γ κ β λ ν ,
which yields w ν > 0 ,   ν ν 1 ,
w ν = μ ν a ( ν ) κ ν γ κ β λ ν + μ ν a ( ν ) κ ν γ κ β λ ν β μ ν λ ν a ( ν ) κ ν γ κ λ ν κ β + 1 λ ν .
We see from (9), (10), (16), and (17) that
w ν μ ν μ ν w ν μ ν p ˜ ν β λ ν w ν κ λ ν κ λ ν .
Since κ ,   κ ,   κ , and κ are positive, and a κ γ < 0 , according to Lemma 1, we can deduce that
κ ν ϵ 2 ν 2 κ ν ,
and
κ λ ν ϵ 2 λ 2 ν κ λ ν ,
for all ϵ 0 , 1 and every sufficiently large ν . Substituting (19) into (18), we obtain
w ν μ ν μ ν w ν μ ν p ˜ ν ϵ 2 β λ 2 ν λ ν κ λ ν κ λ ν w ν = μ ν p ˜ ν + μ ν μ ν w ν ϵ 2 β λ 2 ν λ ν a 1 / γ λ ν a 1 / γ λ ν κ λ ν κ λ ν w ν .
Since a ( ν ) κ ( ν ) γ < 0 , we conclude that
a 1 / γ ν κ ν a 1 / γ λ ν κ λ ν .
Then,
w ν μ ν p ˜ ν + μ ν μ ν w ν ϵ β 2 λ 2 ν λ ν a 1 / γ λ ν w ν a 1 / γ ν κ ν κ λ ν = μ ν p ˜ ν + μ ν μ ν w ν ϵ β 2 λ 2 ν λ ν μ ν a λ ν 1 / γ κ λ ν β γ γ w γ + 1 γ ν .
Because κ > 0 and β γ , there exist constants c 2 > 0 and ν 2 ν 1 such that
κ λ ν c 2 ,
and
κ β / γ 1 λ ν c 2 β / γ 1 , ν ν 2 .
Thus, inequality (20) gives
w ν μ ν p ˜ ν + μ ν μ ν w ν ϵ γ c 2 β / γ 1 2 λ 2 ν λ ν μ ν a λ ν 1 / γ w γ + 1 γ ν .
Using Lemma 2 with
B = μ ν / μ ν , A = ϵ γ c 2 β / γ 1 λ 2 ν λ ν / 2 μ ν a λ ν 1 / γ
and u ν = w ν , we can derive the following inequality:
w ν μ ν p ˜ ν + 2 γ a λ ν μ ν γ + 1 γ + 1 γ + 1 c 2 β γ ϵ λ 2 ν λ ν μ ν γ .
Integrating (23) from ν 3 ν 2 to ν , one arrives at
ν 3 ν μ s p ˜ s 2 γ a λ s μ s γ + 1 γ + 1 γ + 1 c 2 β γ ϵ λ 2 s λ s μ s γ d s w ν 3 ,
which contradicts (15) as ν . This completes the proof. □
Lemma 7.
Suppose that (2) holds. If β γ and there is a μ 1 C 1 ( [ ν 0 , ) , 0 , ) such that μ 1 ( ν ) 0
lim sup ν ν 0 ν M 1 β γ μ 1 s p ˜ s λ s s 3 β / ϵ ˜ 2 γ γ + 1 γ + 1 a s μ 1 s γ + 1 ϵ s 2 μ 1 s γ d s = ,
holds for every ϵ , ϵ ˜ 0 , 1 , then Ω 1 = .
Proof. 
We assume, for the sake of contradiction, that x ν Ω 1 . Now, define a function w 1 ν by
w 1 ν : = μ 1 ν a ( ν ) κ ν γ κ γ ν ,
which yields w 1 ν > 0 . The derivative w 1 ν is given by
w 1 ν = μ 1 ν a ( ν ) κ ν γ κ γ ν + μ 1 ν a ( ν ) κ ν γ κ γ ν γ μ 1 ν a ( ν ) κ ν γ κ ν κ γ + 1 ν .
We see from (10), (25), and (26) that
w 1 ν μ 1 ν p ˜ ν κ β λ ( ν ) κ γ ν + μ 1 ν μ 1 ν w 1 ν γ κ ν κ ν w 1 ν .
From Lemma 3, we deduce that
κ ν ϵ ˜ 3 ν κ ν ,
and hence,
κ λ ν κ ν λ 3 ν ν 3 / ϵ ˜ .
From Lemma 1, we conclude that
κ ν ϵ 2 ν 2 κ ν ,
for all ϵ 0 , 1 . Thus, by (27)–(29), we have
w 1 ν μ 1 ν p ˜ ν κ β γ ν κ β λ ( ν ) κ β ν + μ 1 ν μ 1 ν w 1 ν ϵ γ 2 ν 2 κ ν κ ν w 1 ν = μ 1 ν p ˜ ν κ β γ ν λ 3 ν ν 3 β / ϵ ˜ + μ 1 ν μ 1 ν w 1 ν ϵ γ 2 ν 2 a ν μ 1 ν 1 / γ w 1 1 + γ / γ ν .
Since κ ν > 0 , there exist a ν 3 ν 2 and a constant M 1 > 0 such that
κ ν > M 1 .
Since β γ , then
κ β γ ν > M 1 β γ .
Thus, inequality (30) gives
w 1 ν M 1 β γ μ 1 ν p ˜ ν λ ν ν 3 β / ϵ ˜ + μ 1 ν μ 1 ν w 1 ν ϵ γ ν 2 2 a ν μ 1 ν 1 / γ w 1 1 + γ / γ ν .
Using Lemma 2, with
B = μ 1 ν / μ 1 ν , A = ϵ γ ν 2 / 2 a ν μ 1 ν 1 / γ ,
and u ν = w 1 ν , we can derive the following inequality:
w 1 ν M 1 β γ μ 1 ν p ˜ ν λ ν ν 3 β / ϵ ˜ + 2 γ γ + 1 γ + 1 a ν μ 1 ν γ + 1 ϵ ν 2 μ 1 ν γ .
Integrating (33) from ν 4 ν 3 to ν , one arrives at
ν 4 ν M 1 β γ μ 1 s p ˜ s λ s s 3 β / ϵ ˜ 2 γ γ + 1 γ + 1 a s μ 1 s γ + 1 ϵ s 2 μ 1 s γ d s w 1 ν 4 ,
which contradicts (15) as ν . This completes the proof. □
Lemma 8.
Assume that (2) holds. If 0 < β < γ and there is a nondecreasing function μ C 1 ( [ ν 0 , ) , 0 , ) such that
lim sup ν ν 0 ν μ s p ˜ s 2 β β + 1 β + 1 c 3 β γ a s μ s β + 1 ϵ μ s λ 2 s λ s β d s = ,
holds for every c 3 > 0 ,   ϵ 0 , 1 , then Ω 1 = .
Proof. 
We assume, for the sake of contradiction, that x ν Ω 1 . As in the proof of Lemma 5, we have (9) and (10). The function w ν is defined in (16), and then (17) holds. By (1), (16), and (17), we conclude that
w ν μ ν p ˜ ν + μ ν μ ν w ν β λ ν w ν κ λ ν κ λ ν .
By using (19), we observe that
w ν μ ν p ˜ ν + μ ν μ ν w ν ϵ β 2 λ 2 ν λ ν κ ν β γ / β κ ν γ / β κ λ ν w ν = μ ν p ˜ ν + μ ν μ ν w ν ϵ β 2 λ 2 ν λ ν μ ν a ν 1 / β κ ν β γ / β w β + 1 / β ν .
Note that 0 < β < γ and ( C 1 ) hold. Since a ν 0 , we can conclude that κ 4 ν 0 . This straightforwardly implies that κ ν is nonincreasing. Consequently, there exist constants c 3 > 0 and ν 3 ν 2 such that
κ ν c 3 ,
and
κ ν β γ / β c 3 β γ / β , ν ν 3 .
From (35) and (36), it follows that
w ν μ ν p ˜ ν + μ ν μ ν w ν ϵ β c 3 β γ / β 2 λ 2 ν λ ν μ ν a ν 1 / β w β + 1 / β ν .
Using Lemma 2 with
B = μ ν / μ ν , A = ϵ β c 3 β γ / β λ 2 ν λ ν / 2 μ ν a ν 1 / β ,
and u ν = w ν , It can be deduced from (37) that
w ν μ ν p ˜ ν + 2 β β + 1 β + 1 c 3 β γ a ν μ ν β + 1 ϵ μ ν λ 2 ν λ ν β .
By performing the integration of (38) over the interval ν 5 , ν , we can deduce that
ν 5 ν μ s p ˜ s 2 β β + 1 β + 1 c 3 β γ a s μ s β + 1 ϵ μ s λ 2 s λ s β d s w ν 5 ,
which contradicts (34) as ν . This completes the proof. □
Lemma 9.
If β γ and there is a μ ˜ C ν 0 , , 0 , such that
lim sup ν ν 0 ν M β / γ 1 μ ˜ s s 1 a u u p ˜ v λ v v β / ϵ 1 d v 1 / γ d u μ ˜ + s 2 4 μ ˜ s d s = ,
hold for some ϵ 1 0 , 1 , then Ω 2 = .
Proof. 
We suppose for contradiction that x ν Ω 2 . By integrating (10) from ν to ∞ and utilizing the property that r κ γ 0 , we derive
a ( ν ) κ ( ν ) γ ν p ( s ) g ˜ β λ ( s ) κ β λ ( s ) d s .
As κ > 0 ,   κ > 0 , and κ < 0 , Lemma 3 implies that κ ϵ 1 ν κ for all ϵ 1 0 , 1 . Upon integrating this inequality from λ ν to ν , we find
κ λ ν κ ν λ ν ν 1 / ϵ 1 .
Therefore, (40) becomes
a ν κ ν γ ν p ( s ) g ˜ β λ ( s ) λ s s β / ϵ 1 κ β s d s .
Since κ ν > 0 , then
a ν κ ν γ κ β ν ν p ˜ s λ s s β / ϵ 1 d s ,
or equivalently,
κ ν κ β / γ ν 1 a ν ν p ˜ s λ s s β / ϵ 1 d s 1 / γ .
Upon integrating this inequality from ν to ∞, the result is obtained as
κ ν κ β / γ ν 1 a u u p ˜ s λ s s β / ϵ 1 d s 1 / γ d u .
Now, define
G ν : = μ ˜ ν κ ν κ ν .
Then, G ν 0 , and
G = μ ˜ ν κ ν κ ν + μ ˜ ν κ ν κ ν μ ˜ ν κ ν 2 κ 2 ν = μ ˜ ν κ ν κ ν + μ ˜ ν μ ˜ ν G ν 1 μ ˜ ν G 2 ν .
Hence, by (41), we obtain
G ν μ ˜ ν κ β / γ 1 ν 1 a u u p ˜ s λ s s β / ϵ 1 d s 1 / γ d u + μ ˜ + ν μ ˜ ν G ν 1 μ ˜ ν G 2 ν .
Because κ ν > 0 and β γ , there exist constants M > 0 and ν 2 ν 1 such that
κ ν M , ν ν 2 ,
and
κ β / γ 1 ν M β / γ 1 .
Substituting (44) into (43), we have
G ν M β / γ 1 μ ˜ ν ν 1 a u u p ˜ s λ s s β / ϵ 1 d s 1 / γ d u + μ ˜ + ν μ ˜ ν G ν 1 μ ˜ ν G 2 ν .
Using Lemma 2 with B = μ ˜ + ν / μ ˜ ν ,   A = 1 / μ ˜ ν , and u ν = G ν , we obtain
μ ˜ + ν μ ˜ ν G ν 1 μ ˜ ν G 2 ν μ ˜ + ν 2 4 μ ˜ ν .
Consequently, (44) leads to
G ν M β / γ 1 μ ˜ ν ν 1 a u u p ˜ s λ s s β / ϵ 1 d s 1 / γ d u + μ ˜ + ν 2 4 μ ˜ ν .
Integrating (45) from ν 1 to ν , we have
ν 1 ν M β / γ 1 μ ˜ s s 1 a u u p ˜ v λ v v β / ϵ 1 d v 1 / γ d u μ ˜ + s 2 4 μ ˜ s d s G ν 1 ,
which contradicts (39) as ν . This completes the proof. □
Lemma 10.
If 0 < β < γ and there is a μ ˜ C ν 0 , , 0 , such that
lim sup ν ν 0 ν K β γ γ s ϵ ^ μ ˜ s s 1 a u u p ˜ v λ v v β / ϵ 1 d v 1 / γ d u μ ˜ + s 2 4 μ ˜ s d s = ,
hold for some ϵ 1 0 , 1 ,   ϵ ^ = β γ / ϵ 2 γ , then Ω 2 = .
Proof. 
We begin by assuming, for the sake of contradiction, that x ν Ω 2 . As demonstrated in the proof of Lemma 9, we have (9) and (10). The function G ν is defined in (42), and consequently, (43) holds, which can be expressed as
G ν μ ˜ ν κ β / γ 1 ν 1 a u u p ˜ s λ s s β / ϵ 1 d s 1 / γ d u + μ ˜ + ν μ ˜ ν G ν 1 μ ˜ ν G 2 ν .
Since κ > 0 ,   κ > 0 , and κ < 0 , Lemma 3 implies that
κ ν ϵ 2 ν κ ν ,
for all ϵ 2 0 , 1 , and
1 ϵ 2 ν κ ν κ ν .
Integrating (48) from ν 1 to ν , we obtain
ln ν ν 1 1 / ϵ 2 ln κ ν κ ν 1 ,
which implies
κ ν κ ν 1 ν ν 1 1 / ϵ 2 = κ ν 1 ν 1 1 / ϵ 2 ν 1 / ϵ 2 = K ν 1 / ϵ 2 ,
where K = κ ν 1 ν 1 1 / ϵ 2 . As 0 < β < γ , it follows that 0 < β / γ < 1 , which, combined with (49), yields
κ β / γ 1 ν K β / γ 1 ν β / γ 1 / ϵ 2 .
Hence, inequality (47) becomes
G ν K β / γ 1 μ ˜ ν ν β / γ 1 / ϵ 2 ν 1 a u u p ˜ s λ s s β / ϵ 1 d s 1 / γ d u + μ ˜ + ν μ ˜ ν G ν 1 μ ˜ ν G 2 ν .
Using Lemma 2 with B = μ ˜ + ν / μ ˜ ν ,   a = 1 / μ ˜ ν , and u ν = G ν , we conclude that
μ ˜ + ν μ ˜ ν G ν 1 μ ˜ ν G 2 ν μ ˜ + ν 2 4 μ ˜ ν .
Consequently, (51) leads to
G ν K β / γ 1 μ ˜ ν ν β / γ 1 / ϵ 2 ν 1 a u u p ˜ s λ s s β / ϵ 1 d s 1 / γ d u + μ ˜ + ν 2 4 μ ˜ ν .
Integrating this inequality from ν 2 to ν , we infer that
ν 2 ν K β / γ 1 μ ˜ s s β / γ 1 / ϵ 2 s 1 a u u p ˜ s λ v v β / ϵ 1 d v 1 / γ d u μ ˜ + s 2 4 μ ˜ s d s G ν 1 ,
which contradicts (39) as ν . Thus, we have completed the proof. □

4. Oscillation Criteria

This section builds on the conclusions from our earlier results to establish new conditions for analyzing the oscillatory properties of all solutions to Equation (1). By merging the criteria previously concluded to rule out positive solutions in both scenarios C 1 and C 2 , we can formulate the criteria presented in the following theorems. These criteria will help us ascertain the oscillatory features of the examined equation.
Theorem 1.
Suppose that
lim sup ν ν 0 ν μ s p ˜ s 2 γ a λ s μ s γ + 1 γ + 1 γ + 1 c 2 β γ ϵ λ 2 s λ s μ s γ d s = ,
and
lim sup ν ν 0 ν M β / γ 1 μ ˜ s s 1 a u u p ˜ v λ v v β / ϵ 1 d v 1 / γ d u μ ˜ + s 2 4 μ ˜ s d s = ,
hold. Then, (1) oscillates.
Proof. 
Assume that x Ω 1 Ω 2 . Lemma 4 presents two possible cases for κ and its derivatives. Using Lemmas 6 and 9, it is easy to see that cases ( C 1 ) and ( C 2 ) are excluded by conditions (15) and (39), respectively. This ends the proof. □
Theorem 2.
Suppose that
lim sup ν ν 0 ν M 1 β γ μ 1 s p ˜ s λ s s 3 β / ϵ ˜ 2 γ γ + 1 γ + 1 a s μ 1 s γ + 1 ϵ s 2 μ 1 s γ d s = ,
and
lim sup ν ν 0 ν M β / γ 1 μ ˜ s s 1 a u u p ˜ v λ v v β / ϵ 1 d v 1 / γ d u μ ˜ + s 2 4 μ ˜ s d s = ,
hold. Then, (1) oscillates.
Proof. 
This proof follows a method analogous to that of Theorem 1 and thus was omitted. □
Theorem 3.
Assume that
lim sup ν ν 0 ν μ s p ˜ s 2 β β + 1 β + 1 c 3 β γ a s μ s β + 1 ϵ μ s λ 2 s λ s β d s = ,
and
lim sup ν ν 0 ν K β γ γ s ϵ ^ μ ˜ s s 1 a u u p ˜ v λ v v β / ϵ 1 d v 1 / γ d u μ ˜ + s 2 4 μ ˜ s d s = ,
hold. Then, (1) is oscillatory.
Proof. 
This proof follows a method analogous to that of Theorem 1 and thus was omitted. □
Example 1.
Consider the NDE given by
1 v x ν + g 0 v x 1 / 3 ξ 0 ν + p 0 ν 5 x λ 0 ν = 0 , ν 1 ,
where β γ ,   ξ 0 ,   λ 0 0 , 1 , and p 0 > 0 . By comparing (1) and (52), we observe that
a ν = 1 / v , ξ ν = ξ 0 ν , and λ ν = λ 0 ν .
As a result, we have
g ˜ ν = 1 g 0 c 1 2 / 3 v ,
and
p ˜ ν : = 1 g 0 v c 1 2 / 3 p 0 ν 5 .
Let us consider the function μ ν = ν 4 . Condition (15) can then be expressed as
lim sup ν ν 0 ν μ s p ˜ s 2 γ a λ s μ s γ + 1 γ + 1 γ + 1 c 2 β γ ϵ λ 2 s λ s μ s γ d s = lim sup ν ν 0 ν s 4 1 g 0 s c 1 2 / 3 p 0 s 5 2 2 2 1 ϵ λ 0 s 4 2 s 6 λ 0 2 s 2 λ 0 s 4 d s = lim sup ν ν 0 ν p 0 s g 0 p 0 s 6 c 1 2 / 3 8 ϵ λ 0 4 1 s d s = p 0 8 ϵ λ 0 4 lim sup ν ln v v 0 + g 0 p 0 5 c 1 2 / 3 lim sup ν 1 s 5 g 0 p 0 5 c 1 2 / 3 v 0 5 = ,
which is satisfied when
p 0 > 8 ϵ λ 0 4 .
Similarly, if we consider the function μ 1 ν = ν 4 , condition (24) becomes
lim sup ν ν 0 ν M 1 β γ μ 1 s p ˜ s λ s s 3 β / ϵ ˜ 2 γ γ + 1 γ + 1 a s μ 1 s γ + 1 ϵ s 2 μ 1 s γ d s = lim sup ν ν 0 ν s 4 1 g 0 s c 1 2 / 3 p 0 s 5 λ 0 3 / ϵ ˜ 2 2 2 1 s 4 2 s 6 ϵ s 2 s 4 d s = lim sup ν ν 0 ν p 0 λ 0 3 / ϵ ˜ s g 0 p 0 λ 0 3 / ϵ ˜ s 6 c 1 2 / 3 8 ϵ 1 s d s = p 0 λ 0 3 / ϵ ˜ 8 ϵ lim sup ν ln ν ν 0 + g 0 p 0 λ 0 3 / ϵ ˜ 5 c 1 2 / 3 lim sup ν 1 v 5 g 0 p 0 λ 0 3 / ϵ ˜ 5 c 1 2 / 3 v 0 5 = ,
which is satisfied when
p 0 > 8 ϵ λ 0 3 / ϵ ˜ .
If we consider the function μ ˜ ν = ν , then condition (39) can be expressed as
lim sup ν ν 0 ν M β / γ 1 s s 1 a u u p ˜ v λ v v β / ϵ 1 d v 1 / γ d u μ ˜ + s 2 4 μ ˜ s d s = lim sup ν ν 0 ν s s u u 1 g 0 v c 1 2 / 3 p 0 ν 5 λ 0 1 / ϵ 1 d v 1 / γ d u 1 4 s d s = lim sup ν ν 0 ν s s u u p 0 ν 5 p 0 g 0 v 6 c 1 2 / 3 λ 0 1 / ϵ 1 d v d u 1 4 s d s = lim sup ν ν 0 ν s s p 0 4 u 3 p 0 g 0 5 u 4 c 1 2 / 3 λ 0 1 / ϵ 1 d u 1 4 s d s = lim sup ν ν 0 ν p 0 λ 0 1 / ϵ 1 8 s p 0 g 0 λ 0 1 / ϵ 1 15 s 2 c 1 2 / 3 1 4 s d s = p 0 λ 0 1 / ϵ 1 8 1 4 lim sup v ln v v 0 + p 0 g 0 λ 0 1 / ϵ 1 15 s 2 c 1 2 / 3 lim sup v 1 v 2 p 0 g 0 λ 0 1 / ϵ 1 15 v 0 2 c 1 2 / 3 = ,
which is satisfied when
p 0 > 2 λ 0 1 / ϵ 1 .
Now, according to Theorem 1, if conditions (53) and (55) are satisfied, then (52) exhibits oscillatory behavior. Similarly, Theorem 2 asserts that if conditions (54) and (55) are satisfied, then (54) exhibits oscillatory behavior.
Remark 1.
Consider the NDE given by
1 v x ν + 1 2 v x 1 / 3 0.3 ν + p 0 ν 5 x 0.9 ν = 0 , ν 1 ,
where p 0 > 0 . Now, conditions (53)–(55) are satisfied when p 0 > 24.387 ,   p 0 > 30.107 , and p 0 > 8 , respectively. Consequently, based on Theorem 1, we can deduce that Equation (56) exhibits oscillatory behavior if p 0 > 24.387 .

5. Conclusions

This work discussed the oscillatory features of solutions to a class of fourth-order Emden–Fowler DDEs with a sublinear neutral term. We were first able to obtain criteria that excluded all cases of the corresponding function of the positive solutions. Our approach, leveraging the Riccati technique in various ways, has allowed us to derive comprehensive oscillation criteria for the studied equation. As is clear from the remarks, the new findings are a development and extension of previous relevant findings.
Studying the oscillation of DDEs with a sublinear neutral term has many analytical problems. So, our study imposes some restrictions that limit the application of the results to a wider area of equations of this type. Therefore, we propose—as future work—obtaining an oscillation criterion that does not require constraints α 0 , 1 and lim t g t = 0 and is also considered the noncanonical case.

Author Contributions

Conceptualization, F.M. and W.A.; Methodology, F.M., W.A., O.M. and H.E.-M.; Formal analysis, F.M., W.A. and O.M.; Investigation, O.M. and H.E.-M.; Writing—original draft, F.M. and W.A.; Writing—review & editing, O.M. and H.E.-M. All authors have read and agreed to the published version of the manuscript.

Funding

Princess Nourah bint Abdulrahman University Researchers Supporting Project number (PNURSP2024R157), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia.

Data Availability Statement

Data are contained within the article.

Acknowledgments

The authors express their gratitude and sincere thanks to the editor and the anonymous reviewers for their valuable comments that contributed to improving the final version of the manuscript. The authors would like to thank Princess Nourah bint Abdulrahman University Researchers Supporting Project number (PNURSP2024R157), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Hale, J.K. Theory of Functional Differential Equations; Springer: New York, NY, USA, 1987; Volume 2014. [Google Scholar]
  2. Myshkis, A.D. On solutions of linear homogeneous differential equations of the first order of stable type with a retarded argument. Mat. Sb. 1951, 70, 641–658. [Google Scholar]
  3. Norkin, S.B. Introduction to the Theory and Application of Differential Equations with Deviating Arguments; Academic Press: Cambridge, MA, USA, 1973. [Google Scholar]
  4. Cooke, K.L. Differential Difference Equations; Academic Press: New York, NY, USA, 1963. [Google Scholar]
  5. Braun, M. Qualitative Theory of Differential Equations: Differential Equations and Their Applications; Springer: New York, NY, USA, 1993. [Google Scholar]
  6. Gyori, I.; Ladas, G. Oscillation Theory of Delay Differential Equations with Applications; Clarendon Press: Oxford, UK, 1991. [Google Scholar]
  7. Zafer, A. Oscillatory and Nonoscillatory Properties of Solutions of Functional Differential Equations and Difference Equations; Iowa State University: Ames, IA, USA, 1992. [Google Scholar]
  8. Grace, S.R.; Dzurina, J.; Jadlovska, I.; Li, T. On the oscillation of fourth-order delay differential equations. Adv. Differ. Equ. 2019, 2019, 118. [Google Scholar]
  9. Jiang, Y.; Jiang, C.; Li, T. Oscillatory behavior of third-order nonlinear neutral delay differential equations. Adv. Differ. Equ. 2016, 2016, 207. [Google Scholar] [CrossRef]
  10. Zhang, C.; Li, T.; Saker, S.H. Oscillation of fourth-order delay differential equations. J. Math. Sci. 2014, 201, 296–309. [Google Scholar] [CrossRef]
  11. Erbe, L.; Kong, Q.; Zhang, B.G. Oscillation Theory for Functional Differential Eqautions; Marcel Dekker: New York, NY, USA, 1995. [Google Scholar]
  12. Ladde, G.S.; Lakshmikantham, V.; Zhang, B. Oscillation Theory of Differential Equations with Deviating Arguments; Marcel Dekker: New York, NY, USA, 1987. [Google Scholar]
  13. Agarwal, R.P.; Bohner, M.; Li, T. Oscillation of second-order Emden-Fowler neutral delay differential equations. Ann. Math. 2014, 193, 1861–1875. [Google Scholar] [CrossRef]
  14. Chatzarakis, G.E.; Grace, S.R.; Jadlovská, I. Oscillation criteria for third-order Emden-Fowler differential equations with unbounded neutral coefficients. Complexity 2019, 2019, 5691758. [Google Scholar] [CrossRef]
  15. Tian, H.; Guo, R. Some Oscillatory Criteria for Second-Order Emden–Fowler Neutral Delay Differential Equations. Mathematics 2024, 12, 1559. [Google Scholar] [CrossRef]
  16. Feng, Q.; Zheng, B. Oscillation Criteria for Nonlinear Third-Order Delay Dynamic Equations on Time Scales Involving a Super-Linear Neutral Term. Fractal Fract. 2024, 8, 115. [Google Scholar] [CrossRef]
  17. Agarwal, R.P.; Bohner, M.; Li, T. Oscillation of second-order differential equations with a sublinear neutral term. Carpathian J. Math. 2014, 30, 1–6. [Google Scholar]
  18. Baculíková, B.; Sudha, B.; Thangavelu, K.; Thandapani, E. Oscillation of second order delay differential equations with nonlinear nonpositive neutral term. Math. Slovaca 2022, 72, 103–110. [Google Scholar] [CrossRef]
  19. Zhang, C.; Şenel, M.T.; Li, T. Oscillation of second-order half-linear differential equations with several neutral terms. J. Appl. Math. Comput. 2014, 44, 511–518. [Google Scholar] [CrossRef]
  20. Tamilvanan, S.; Thandapani, E.; Džurina, J. Oscillation of second order nonlinear differential equations with sub-linear neutral term. Differ. Equ. Appl. 2017. [Google Scholar] [CrossRef]
  21. Wu, Y.; Yu, Y.; Xiao, J.; Jiao, Z. Oscillatory behaviour of a class of second order Emden-Fowler differential equations with a sublinear neutral term. Appl. Math. Sci. Eng. 2023, 31, 2246098. [Google Scholar] [CrossRef]
  22. El-Nabulsi, R.A.; Moaaz, O.; Bazighifan, O. New Results for Oscillatory Behavior of Fourth-Order Differential Equations. Symmetry 2020, 12, 136. [Google Scholar] [CrossRef]
  23. Moaaz, O.; Muhib, A.; Zakarya, M.; Abdel-Aty, A. Delay differential equation of fourth-order: Asymptotic analysis and oscillatory behavior. Alex. Eng. J. 2022, 61, 2919–2924. [Google Scholar] [CrossRef]
  24. Li, T.; Rogovchenko, Y.V. On asymptotic behavior of solutions to higher-order sublinear Emden–Fowler delay differential equations. Appl. Math. Lett. 2017, 67, 53–59. [Google Scholar] [CrossRef]
  25. Graef, J.R.; Grace, S.R.; Tunç, E. Oscillatory behavior of even-order nonlinear differential equations with a sublinear neutral term. Opusc. Math. 2019, 39, 39–47. [Google Scholar] [CrossRef]
  26. Baculikova, B.; Dzurina, J.; Li, T. Oscillation Results for Even-Order Quasilinear Neutral Functional Differential Equations. Electron. Differ. Equ. 2011, 2011, 1–9. [Google Scholar]
  27. Xing, G.; Li, T.; Zhang, C. Oscillation of higher-order quasi-linear neutral differential equations. Adv. Differ. Equ. 2011, 2011, 45. [Google Scholar] [CrossRef]
  28. Alnafisah, Y.; Masood, F.; Muhib, A.; Moaaz, O. Improved Oscillation Theorems for Even-Order Quasi-Linear Neutral Differential Equations. Symmetry 2023, 15, 1128. [Google Scholar] [CrossRef]
  29. Agarwal, R.P.; Grace, S.R.; O’Regan, D. Oscillation Theory for Difference and Functional Differential Equations; Kluwer Academic: Dordrecht, The Netherlands, 2000. [Google Scholar]
  30. Zhang, C.; Agarwal, R.; Bohner, M.; Li, T. New results for oscillatory behavior of even-order half-linear delay differential equations. Appl. Math. Lett. 2013, 26, 179–183. [Google Scholar] [CrossRef]
  31. Kiguradze, I.T.; Chanturiya, T.A. Asymptotic Properties of Solutions of Nonautonomous Ordinary Differential Equations; Kluwer Academic Publishers: Dordrecht, The Netherlands, 1993. [Google Scholar]
  32. Masood, F.; Moaaz, O.; Santra, S.S.; Fernandez-Gamiz, U.; El-Metwally, H.A. Oscillation theorems for fourth-order quasi-linear delay differential equations. AIMS Math. 2023, 8, 16291–16307. [Google Scholar] [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Masood, F.; Albalawi, W.; Moaaz, O.; El-Metwally, H. Oscillatory Features of Fourth-Order Emden–Fowler Differential Equations with Sublinear Neutral Terms. Symmetry 2024, 16, 933. https://doi.org/10.3390/sym16070933

AMA Style

Masood F, Albalawi W, Moaaz O, El-Metwally H. Oscillatory Features of Fourth-Order Emden–Fowler Differential Equations with Sublinear Neutral Terms. Symmetry. 2024; 16(7):933. https://doi.org/10.3390/sym16070933

Chicago/Turabian Style

Masood, Fahd, Wedad Albalawi, Osama Moaaz, and Hamdy El-Metwally. 2024. "Oscillatory Features of Fourth-Order Emden–Fowler Differential Equations with Sublinear Neutral Terms" Symmetry 16, no. 7: 933. https://doi.org/10.3390/sym16070933

APA Style

Masood, F., Albalawi, W., Moaaz, O., & El-Metwally, H. (2024). Oscillatory Features of Fourth-Order Emden–Fowler Differential Equations with Sublinear Neutral Terms. Symmetry, 16(7), 933. https://doi.org/10.3390/sym16070933

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