In this paper, we consider a new fractional differential system on an unbounded domain D α u(t) + ϕ(t, v(t), D γ 1 v(t)) = 0, t ∈ [0, +∞), α ∈ (2, 3], D β v(t) + ψ(t, u(t), D γ 2 u(t)) = 0, t ∈ [0, +∞), β ∈ (2, 3], subject to the conditions I 3-α u(t)| t=0 = 0, D α-2 u(t)| t=0 = h 0 g 1 (s)u(s) ds, D α-1 u(+∞) = Mu(ξ) + a, I 3-β v(t)| t=0 = 0, D β-2 v(t)| t=0 = h 0 g 2 (s)v(s) ds, D β-1 v(+∞) = Nv(η) + b. The nonlinear terms ϕ and ψ are dependent on the fractional derivative of lower order γ i ∈ (0, 1), i = 1, 2, which creates additional complexity to verify the existence of solutions. Moreover, a proper choice of Banach space allows the solutions to be defined on the half-line. From some standard fixed point theorems, sufficient conditions for the existence and uniqueness of solutions to boundary value problems are developed. Finally, the main result is applied to an illustrative example.