Cosmological singularities in f ( T, φ ) gravity

The pursuit of understanding the mysteries surrounding dark energy has sparked signiﬁcant interest within the ﬁeld of cosmology. While conventional approaches, such as the cosmological constant, have been extensively explored, alternative theories incorporating scalar ﬁeld-based models and modiﬁed gravity have emerged as intriguing avenues. Among these, teleparallel theories of gravity, speciﬁcally the f ( T, φ ) formulation, have gained prominence as a means to comprehend dark energy within the framework of teleparallelism. In this study, we investigate two well-studied models of teleparallel dark energy and examine the presence of cosmological singularities within these scenarios. Using the Goriely-Hyde procedure, we examine the dynamical systems governing the cosmological equations of these models. Our analysis reveals that both models exhibit Type IV singularities, but only for a limited range of initial conditions. These results could indicate


Introduction
Observations of the late-time acceleration of the Universe came as a surprise to the cosmological community [1].Since then, extensive efforts have been dedicated to explaining this expansion.Standard approaches, such as the Cosmological constant [1,2,3,4,5], as well as more exotic scenarios like Modified gravity theories [6,7,8], and recent proposals for the direct detection of dark energy [9], have been pursued.One fascinating avenue for understanding dark energy is through Quintessence, where a scalar field drives the late-time cosmic acceleration of the universe [10,11,12,13,14,15,16,17,18,19,20].Quintessence is particularly interesting as it represents the simplest scalar field dark energy scenario that avoids issues like ghosts or Laplacian instabilities.In quintessence models, the acceleration of the universe is driven by a slowly varying scalar field with a potential V (φ), similar to the mechanism of slow-roll inflation.However, in this case, contributions from nonrelativistic matter, such as baryons and dark matter, cannot be neglected.
It is worth noting that simple models of Quintessence have been shown to be in conflict with the current H0 tension [21,22,23], suggesting that simple Quintessence models may perform worse than Λ-CDM models in light of the current H0 data [24].This leads one to consider other more exotic possibilities for scalar field dark energy models and one such possibility is to consider models in Teleparallel gravity.Teleparallel gravity, a theory based on torsion, provides an alternative description of gravity [25,26,27,28,29,30,31,32,33], where gravitation is mediated by torsion.In this approach, the La-grangian density of Teleparallel Equivalent of General Relativity (TEGR) is proportional to the torsion scalar T .In TEGR, the tetrad field and spin connection pair replace the metric tensor and Levi-Civita connection, respectively, while the teleparallel connection replaces the usual connection [31,30].Consequently, at the level of the dynamical equations, curvature-based gravitational theories are equivalent to tensor-based theories [30,34].By introducing an arbitrary function f (T ) in place of the torsion scalar T , a generalization of TEGR known as f (T ) gravity is obtained [35,36,37,38,39,40,41], leading to new cosmological models.In this framework, the tetrad fields, which form the orthogonal basis for the tangent space, serve as the dynamical variables of teleparallel gravity.The torsion tensor is constructed using the first derivative of the tetrad product.The field equations are derived by varying the action with respect to the tetrad fields, while the spin connection preserves the local Lorentz invariance and contributes to the equations of motion.For further exploration of f (T ) gravity, refer to [42,43,44,45,46,47,48].The investigation of scalar-torsion theories with non-minimal coupling between the torsion scalar and a scalar field was carried out in the context of dark energy [49,50], including studies with arbitrary non-minimal coupling functions and tachyon terms for the scalar field [51,52].Another extension of f (T ) gravity is the generalized scalar-torsion f (T, φ) gravity, where φ represents the canonical scalar field, and the gravitational action incorporates a non-minimal coupling between the scalar field and torsion scalar [53].Additionally, within the covariant teleparallel framework, a new class of theories has been proposed, where the action depends on the scalar field and an arbitrary function of the torsion scalar [54].
Recently, a significant amount of research has been dedicated to exploring the various types of cosmological singularities that may occur in the present and distant future of the Universe [55,56,57,58,59,60,61,62,63,64,65,66,67,68].However, it is often challenging to classify and study cosmological singularities in highly unconventional cosmologies influenced by considerations of quantum gravity or phenomenology.Traditional methods may not be applicable in these cases.Therefore, alternative approaches are necessary to identify cosmological singularities in exotic cosmologies.In this context, the Goriely-Hyde procedure, a particular method in dynamical systems, can be extremely useful [69].Understanding the singularity structure of dynamical systems is an intriguing aspect, especially when these systems describe significant physical phenomena.Although various approaches have been proposed to inves-tigate the singularity structure of autonomous dynamical systems, the Goriely-Hyde procedure has proven particularly valuable for cosmological studies due to the abundance of interesting dynamical systems in cosmology [70].Previous applications of the Goriely-Hyde method have explored finite and non-finite time singularities in specific quintessence models [71,72,73].However, a comprehensive analysis of cosmological singularities teleparallel models of dark energy using this approach is still lacking, and our work aims to address this gap.The study of the cosmological dynamics and stability of f (t, φ) dark energy was conducted in [74], while an analysis of scalar perturbations was performed in [75].A recent full dynamical systems analysis of f (t, φ) dark energy for two particular models was done in [76] and we intend to use the dynamical systems approach developed there to pursue our singularity analysis.In Section II, we provide a concise overview of Teleparallel gravity method while in section III we provide a brief review of the Goriely-Hyde method.In Section IV we apply the Goriely-Hyde method to two particular models and demonstrate the diverse characteristics of singularities in f (T, φ) models, including both finite and infinite-time occurrences which can occur for both f (t, φ) models considered in [76] .Subsequently, in Section V, we consider two well-motivated ansatz for the Hubble parameter and classify the types of cosmological singularities (Types I-IV) that can arise within these regimes.Finally, we conclude our work in Section VI.

Teleparallel gravity
General relativity (GR) can account for most observed phenomena with appropriate modifications considered in the matter sector.In this context, the widely accepted concordance model combines Λ-CDM cosmology with inflation.However, the enigmatic nature of certain particle species remains a puzzle despite significant progress in physics beyond the standard model of particle physics.It is also plausible that the standard model of particle physics might not require substantial restructuring to address these observational challenges.Instead, it could be the gravitational sector that requires further examination.This could involve extensions of GR or modifications beyond GR as alternatives to its original formulation.The scientific literature has witnessed numerous proposals for new theories of gravity, motivated by various phenomena, theoretical approaches, or even quantum physics.One intriguing possibility that has garnered increasing attention in recent decades is teleparallel gravity, where torsion replaces curvature as the mechanism responsible for generating gravitational fields.This theory replaces the traditional Levi-Civita connection, which is curvature-based, with a teleparallel connection based on torsion.Numerous publications on this topic have emerged in the literature.Among the theories arising from torsion-based approaches to gravity is the teleparallel equivalent of general relativity (TEGR), which is dynamically equivalent to GR and thus indistinguishable from it through classical experiments.In teleparallel gravity (TG), one typically assumes an action of the form: Here, the gravitational part S g of the action depends on the tetrad e A µ and the spin connection ω A Bµ , while the matter part depends on the tetrad e A µ and arbitrary matter fields χ I , but not on the spin connection [77,78].This is because we assume that the hypermomentum vanishes, thereby preventing this coupling.Introducing a dependence on spin would effectively introduce a second matter tensor, resulting from the variation of the matter Lagrangian with respect to the spin connection.The variation of the matter part of the action, after integration by parts to eliminate derivatives acting on field variations, can be expressed as follows: Here, Ω I = 0 represents the matter field equations, and Θ A µ denotes the energy-momentum tensor.The corresponding variation of the gravitational action takes the form: The tensors W A µ and Y A Bµ arise from the variation and integration by parts, with their specific form depending on the particular theory under consideration.The explicit expression for W A µ can be found for several theories in [79].For brevity, we omit Y A Bµ here, as it turns out to be redundant in deriving the field equations, which can be entirely determined from W A µ alone.Furthermore, by varying with respect to the tetrad, one can derive the field equations: An alternative representation of the field equations, more commonly used, is obtained by transforming the first index into a spacetime index with the tetrad while lowering the second index: This yields the field equations in the form: However, deriving the field equations for the spin connection is more complex, as it must satisfy the conditions of being flat, R α βµν = 0, and metric-compatible, ∇ α g µν = 0, by definition.Various approaches exist to maintain these properties during the variation procedure [80,81].
For considering cosmological scenarios in teleparallel theories, it is helpful to consider an FLRW metric of the form [79] where N (t) and a(t) represent the lapse function and scale factor respectively.In the case of flat universes (k = 0) one can write and this results in the diagonal tetrad which turns out to be in the Weitzenbock gauge for the extensions to TEGR.An important remark here is that the above tetrad (with vanishing spin connection) is the only one that has the property that both the tetrad and the teleparallel connection obey cosmological symmetries for flat FLRW.One can also relax the condition that the teleparallel connection enjoys the symmetries of cosmology, but then, the corresponding cosmological equations would not respect the symmetries of cosmology.If we use the diagonal tetrad ( 9) in Cartesian coordinates one can, for example, find the modified FLRW equations for f (T ) gravity as where dots are derivatives with respect to time, so that ḟT = f T T Ṫ .One can further obtain modified FLRW equations for other teleparallel appraoches, like f (T, B) gravity being Or for the Teleparallel Gauss-Bonet models being (fixing the gauge such that N = 1) where While there are a multitude of approaches of dealing with such exotic cosmological systems like reconstruction methods or Noether-symmetries approaches, dynamical systems methods are also an efficient way to understand the dynamics of the models.Dynamical systems allow for the extraction of key features of the cosmology without solving the evolution equations directly (in an exact form).Thus, it then becomes possible to describe the overall nature of the gravitational theory and henceforth determine whether the model can generate a viable cosmological evolution.This therefore serves as a very useful tool especially in models where it is difficult to extract any cosmological solutions from directly solving the field equations such as in f(R) gravity.In the cases we have considered so far, one can for example write the equations (10) as an autonomous dynamical system using the variables These variables were considered in [82] and the authors considered the following scenarios: (i) absence of matter fluids and (ii) presence of dust and radiation components.Furthermore, the case when the parameter m = − Ḧ H 3 takes on constant values m = 0 (quasi-de Sitter evolution) and m = − 9 2 (matter dominated evolution) was explored.One can also write (11) in the dynamical systems method using the phase space variables The cosmological dynamics of f (T, T G ) gravity was investigated in [83].In particular, the model where α 1,2 = 0 are constants, was studied.In this case, the presence of a perfect dust fluid was assumed and the following dimensionless phasespace parameters were defined While we have briefly discussed several approaches to teleparallel gravity and their status quo in cosmology, what we are most interested in this work are the f (T, φ) models and we shall now discuss them in more detail.The action of TEGR (Teleparallel Equivalent of General Relativity) can be generalized to f (T ) gravity by introducing a scalar field φ.The action, including matter and radiation, can be expressed as [75,76]: Here, e = det[e A µ ] = √ −g represents the determinant of the tetrad field.The tetrad field, e A µ , A = 0, 1, 2, 3, is related to the metric tensor g µν and the Minkowski tangent space metric η AB as g µν = η AB e A µ e B ν , where η AB = (−1, 1, 1, 1).The tetrad satisfies the orthogonality condition e µ A e B µ = δ B A , and the spin connection is denoted by ω A Bµ .The function f (T, φ) represents an arbitrary function of the scalar field φ and the torsion scalar T , while X = −∂ µ φ∂ µ φ/2 represents the kinetic term of the field.This general action includes nonminimally coupled scalar-torsion gravity models with the coupling function f (T, φ), f (T ) gravity, and minimally coupled scalar field.
For a flat FLRW space-time background, the field equations derived from the action are [76]: In these equations, the Hubble parameter is denoted as H ≡ ȧ a , where an overdot represents a derivative with respect to cosmic time t.The energy density for matter and radiation are denoted as ρ m and ρ r respectively, and the pressure at the radiation era is p r .The torsion scalar T is given by T = 6H 2 .The non-minimal coupling function f (T, φ) is defined as [54]: where V (φ) is the scalar potential and G(T ) is an arbitrary function of the torsion scalar.
In the matter-dominated era, ω m = pm ρm = 0, and in the radiation era, ω r = pr ρr = 1/3.In this case, Eqs. ( 18)-( 20) reduce to: The modified Friedmann equations, taking into account dark energy, become: Comparing Eq. ( 22) with Eq. ( 25), and Eq. ( 23) with Eq. ( 26), we can extract the energy density (ρ de ) and pressure (p de ) for the dark energy sector: For simplicity we set P (φ) = 1 and consider the well studied exponential potential, V (φ) = V 0 e −λφ , where λ is a constant.In order to proceed further and really carry out the analysis we want to, we need a form for G(T ) and here we will be considering two forms which were studied in [76] and will be carrying out the Goriely-Hyde analysis on both of these models 1 .
1 It is worth discussing any effects of the separation of T from φ here in the action (21).If one, for example, considers actions like those in [50] where terms like T φ 2 come into play then one could potentially expect that results on singularities could be affected in some ways, as one would not be wrong to think that such a coupling between the torsion scalar and field terms could have some significant outcomes.Although a definitive answer on this aspect would need one to do a proper analysis similar to the one we have done for the models considered in our paper, we do feel very interesting results may await an endeavour like this.
3 The Goriely-Hyde Procedure The Goriely-Hyde technique [69] offers an elegant method for identifying finite-time singularities in dynamical systems.The procedure can be summarized as follows: -We begin by considering a dynamical system governed by n differential equations given by: where i = 1, 2, ..., n.Here, t represents time, but in quintessence models, it can be better represented as the number of e-foldings, denoted by N .We identify the parts of the equation f i that become significant as the system approaches the singularity.These significant parts are referred to as "dominant parts" [69].Each dominant part represents a mathematically consistent truncation of the system, denoted as f i.Consequently, the system can be expressed as: -Without loss of generality, the variables x i near the singularity can be represented as: where τ = t − t c , and t c is an integration constant.By substituting equation (4) into equation ( 3) and equating the exponents, we can determine the values of p i for different i, which collectively form the vector p = (p 1 , p 2 , ..., p n ).Similarly, we calculate the values of a i to form the vector a = (a 1 , a 2 , ..., a n ).
It is worth noting that if a comprises solely real entries, it corresponds to finite-time singularities.On the other hand, if a contains at least one complex entry, it may lead to non-finite-time singularities.Each (a i , p i ) set is referred to as a dominant balance of the system.-Next, we compute the Kovalevskaya matrix defined as: After obtaining the Kovalevskaya matrix, we evaluate it for different dominant balances and determine the eigenvalues.If the eigenvalues take the form (−1, r 2 , r 3 , ..., r n ), where r 2 , r 3 , ... > 0, then the singularity is regarded as general and will occur regardless of the initial conditions of the system.Conversely, if any of the eigenvalues r 2 , r 3 , ... are negative, the singularity is considered local and will only occur for certain sets of initial conditions.

Model I
In the first model, we consider a specific form for G(T ) as given in [84,76]: where β is a constant and T 0 represents the value of T at the initial epoch.This model, which has been investigated in [84], exhibits physically favorable critical points and offers an interesting approach for modeling the evolution of the Universe.By substituting this expression into Eqs.( 27)-( 28), the effective dark energy density and pressure terms are reduced to: In order to analyze the dynamics of the scalar-torsion f (T, φ) gravity model, [76] introduced a set of dimensionless phase space variables to represent the system in an autonomous form.These variables are defined as2 follows: These dimensionless variables allow for a simplified representation of the system's dynamics and facilitate the analysis of the scalar-torsion f (T, φ) gravity model.Using these variables, one can finally write the cosmological equations of this model as a dynamical system as follows : For our analysis , we would be considering λ to be a constant which is not equal to zero (which would again mean that we are considering an exponential potential form as we remarked earlier).Furthermore, we consider that 2u >> 3z −2 ( which can be justified considering the forms of z and u we have described earlier ).This would allow us to write the dynamical equations as Now we are in the right position to start off our singularity analysis.The first truncation that we consider is given by We can now write the Kovalevskaya matrix to be Using the dominant balances we introduced in (44), we can now plug them into the Kovalevskaya matrix (45) to get the eigenvalues to be We note that all the other eigenvalues besides the initial -1 are also negative, which means that according to the Goriely-Hyde method the singularities of this system with regards to this truncation can occur only for a limited set of initial conditions.Coupled with the fact that the dominant balances (44) have complex entries, this would mean that this truncation tells us that the singularities for this system may occur in non-finite time.
The second truncation that we consider is given by Using the ansatz of the Goriely-Hyde method, we get the exponents to be p = (−1, −1, −1, −3/2) from which we can get the dominant balances to be .
We again see that the balances have complex entries 3 while the the Kovalevskaya matrix can be written as Using the dominant balances we introduced in ( 48), we can now plug them into the Kovalevskaya matrix (49) to get the eigenvalues to be We note that as one of the other eigenvalues (-1.5) besides the initial -1 is also negative, according to the Goriely-Hyde method the singularities of this system with regards to this truncation can occur only for a limited set of initial conditions.Coupled with the fact that the dominant balances (48) have complex entries, this would mean that this truncation tells us that the singularities for this system may occur in non-finite time.
The third truncation that we consider is given by Using the ansatz of the Goriely-Hyde method, we get the exponents to be p = (1/2, 1/2, 1, 1/4) from which 3 At this point we would like to highlight that complex entries in â observed for the previous truncation and this one ( and which will be observed for Model II as well for a few truncations) are completely consistent with the fact that the system consists of expansion normalized variables which are real.As mentioned in section 2, complex entries for various a suggest that the singularities will be non-finite time in nature and hence these quantities taking up complex values is consistent with the analysis as shown in [69].Similar case has been for various cosmological systems (for example, see [71,72]) we can get the dominant balances to be .
We can now write the Kovalevskaya matrix to be Using the dominant balances we introduced in (52), we can now plug them into the Kovalevskaya matrix (53) to get the eigenvalues to be We note that as two of the other eigenvalues besides the initial -1 are also negative, according to the Goriely-Hyde method the singularities of this system with regards to this truncation can occur only for a limited set of initial conditions.But in this case we see something which we didn't in the previous two truncations ; the dominant balances (52) do not have complex entries.This means that this truncation tells us that it is definitely possible to have singularities occuring in finite-time for this particular model.While we can go on and evaluate more truncations, we find that in no other truncation would we see something which we have not observed already in these three truncations.Namely that there is no truncation for this system for which the eigenvalues besides -1 are all positive and so it does seem like for this model the singularities will not be general and can only happen for a limited set of initial conditions.

Model II
In this scenario, we consider the function G(T ) to be of the form G(T ) = T + αT 2 , where α is a constant [88].This represents a slight extension beyond the Teleparallel Equivalent of General Relativity (TEGR), as α = 0 corresponds to the TEGR model.For this particular G(T ), Eqs. ( 36)-( 37) can be expressed as follows: To establish an independent dynamical system, we introduce dimensionless variables defined as: Consequently, the corresponding dynamical system can be obtained as, We again consider λ to be a constant here, which would mean that we are interested in exponential potentials.Furthermore, we assume that 2u >> z − 1 which is again not hard to justify considering the definitions of these quantities in (57).By taking these considerations into account, the dynamical system takes the form We can now start with the Goriely-Hyde analysis of this system, with the first truncation that we consider being Using the ansatz of the Goriely-Hyde method, we get the exponents to be p = (−1, −1, −1, −1) from which we can get the dominant balances to be .
We can now write the Kovalevskaya matrix to be Using the dominant balances we introduced in (65), we can now plug them into the Kovalevskaya matrix (66) to get the eigenvalues to be As we have one of the eigenvalues besides -1 also being negative, this truncation tells us that the singularities that could appear for this model would also only be occuring for a limited set of initial conditions for the variables.Furthermore given that the dominant balances (65) all have real entries then this would mean that the singularities only appear in finite time.
The second truncation that we would be considering is given by f Using the ansatz of the Goriely-Hyde method, we get the exponents to be p = (−1, −1, −1, −1) from which we can get the dominant balances to be .
We can now write the Kovalevskaya matrix to be Using the dominant balances we introduced in (69), we can now plug them into the Kovalevskaya matrix (70) to get the eigenvalues to be r ∼ (−1, 1.6, 3.7, −0.2) We again see that there are eigenvalues besides -1 which are negative, which again suggests that the model may not have general singularities.Furthermore this truncation also suggests that singularities could take place in non-finite time as shown by the complex entries in the dominant balance (69).While we can again go on for more truncations, what we have found out is that the other truncations do not offer anything new other than what we have seen so far.Namely, no truncation suggests that the model can allow for general singularities and so we are not going to be evaluating for more truncations here.

Physical classification of the singularities
Until now, we have discussed the singularity structure within the dark energy scenario from a dynamical perspective.However, it is insufficient to merely acknowledge the existence of singularities in this system from a physical standpoint.Thus, it becomes necessary to appropriately classify the potential types of singularities that could occur in this model.Various types of physical singularities for cosmology at a specific time t = t s , where t s represents the occurrence of the singularities, can be classified as follows [57,89]: -Type I ("Big Rip"): In this case, the scale factor a, effective energy density ρ eff , and effective pressure density p eff diverge.-Type II ("Sudden/Quiescent singularity"): In this case, p eff diverges, as well as the derivatives of the scale factor beyond the second derivative.-Type III ("Big Freeze"): In this case, the derivative of the scale factor from the first derivative onwards diverges.-Type IV ("Generalized sudden singularities"): In this case, the derivative of the scale factor diverges from a derivative higher than the second.
Among these classifications, Type I singularities are considered strong singularities since they have the ability to distort finite objects, while singularities of Type II, Type III, and Type IV are regarded as weak singularities as they cannot be perceived as either the beginning or the end of the universe.Although there are other minor types of singularities, such as Type V singularities or "w" singularities, we will focus solely on Type I to Type IV singularities here.The most general form of the Hubble parameter for investigating singularities within the aforementioned classified types is expressed as [72]: Here, f 1 (t) and f 2 (t) are assumed to be nonzero regular functions at the time of the singularity, and similar conditions apply to their derivatives up to the second order.Additionally, ǫ is a real number.It is not mandatory for the Hubble parameter (34) to be a solution to the field equations; however, we will consider this case and explore the implications of this assumption on the singularity structure based on our dynamic analysis.First, we observe that none of the variables x, y, or z as defined in (10) can ever become singular for any cosmic time value.The singularities that can occur considering the Hubble parameter as defined in (34) are as follows: -For ǫ < −1, a big rip singularity occurs.
Another ansatz useful for classifying singularities was introduced in [61] whereby the scale factor was written as: where g(t) and f (t) and all their higher-order derivatives with respect to cosmic time are smooth functions of the cosmic time.For this ansatz, according to the values of the exponent ǫ, one can have the following singularities: -For ǫ < 0, a Type I singularity occurs.
Again, it is not mandatory for the scale factor in equation (73) to necessarily be a solution to the field equations, but we would like to consider this and equation (72) in order to gain a well-motivated understanding of the types of cosmological singularities we can encounter in the various models we have discussed so far.
To proceed further, we need to express the expansion normalized variables that we defined for both models in terms of the Hubble parameter alone.To do this, we realize that we need to express the potential and the derivative of the field parameter in each case in terms of the Hubble parameter as these are the quantities on which the expansion normalized variables really depend in both the scenario ( in this scenario we are talking about representing the x and y variables in both cases in terms of the Hubble parameter).For the model G(T ) = βT ln T T0 (33), we have While the potential for this case be written as For the model G(T ) = T + αT 2 , we have the same quantities to be Using these one can express the dynamical variables used in the Goriely-Hyde analysis of both the models (( 36),( 57)) completely in terms of the Hubble parameter ( we will not write the variables out explicitly here as they have quite long expressions) and now we can use both the ansatz ( 72)- (73) to see under what conditions will the variables blow up.Remember that we do not want the dynamical variables to blow up and the values of the exponents of the ansatz for which they do not blow up will tell us the singularities which will be possible for these models.The interesting conclusion that actually comes out when one puts both the ansatz into the dynamical variables is that only Type IV singularities are possible for both models.None of Type I, Type II or Type III singularities can occur for any of the models for any of the ansatz ( 72)-( 73) while Type IV singularities do take place for both the models, for any of the ansatz'.This is quite an interesting behaviour which to the best of our knowledge has only been shown in f (T, φ) theories, in that one is only observing Type IV singularities for both of the models considered.This leads one to speculate the possibility that f (T, φ) gravity may be better suited for cosmology than some of their other modified gravity counterparts as the theory is only admitting weak singularities.Furthermore, given the analysis from the Goriely-Hyde procedure, one is lead to conclude that such singularities can only occur for a limited set of initial conditions and may occur in finite or even non-finite time.

Concluding remarks
In this paper, we have considered a well studied formulation of teleparallel dark energy in the form of f (T, φ) gravity, where the scalar field drives the expansion of the universe.We considered two particular well studied models of this theory and probed cosmological singularities for both the scenarios.For this endeavor, we used a method pioneered by the works of Odintsov in recent years, in which we applied the Goriely-Hyde procedure to the various dynamical systems by which the cosmological equations of these three models could be described.This allowed us to make predictions about whether singularities in these scenarios would be strongly dependent on initial physical conditions and whether they could happen in finite or nonfinite times.After this, we employed two very well-motivated ansatz' for the Hubble parameter and the scale factor to reach the conclusion that one can only have Type IV singularities for both of the models considered in our work, that too only for a limited set of initial conditions.This work propels one to think in the direction that f (T, φ) theories may only allow for weak cosmological singularities, which may make them better placed than some of the other modified gravity based dark energy regimes which allow for more singularities and also those of the stronger types.