Let be a weak geometric p-rough path, , with increments indexed by words of length at most . Choose a control function for rough paths with . We use the standard sufficient vector-field hypothesis , , with globally bounded regularity norms; smooth bounded fields are a special case. Regard as the differential operator , and put .
A pathlevel solution of a rough differential equation is a continuous state path , with prescribed , satisfying a local expansion
for some on sufficiently small control intervals. One can take under the stated regularity. This definition specifies the state path, without requiring its full rough-path enhancement as additional data. For , the displayed expansion reads
Repeated indices are summed.
We use the following discrete form of the Davie lemma. For a finite ordered grid and a continuous superadditive control, suppose nonnegative errors vanish on adjacent grid points and, for every grid triple , satisfy
There are and , depending on but not on the grid, such that whenever . This is the discrete higher-order error estimate: splitting close to half the control size makes the two -power errors contract; a crossing single grid interval contributes zero adjacent error.
Define the local Euler map
For small it has Lipschitz constant at most . Taylor expansion, the Chen identity and the shuffle identity for iterated integrals give the composition estimate
Here is the cancellation underlying the estimate. Assign degree to a level- increment. Expanding each and grouping terms through degree gives the differential-operator coefficients multiplying . Products of coefficients from the first increment are combined by the shuffle identity, which is the iterated Leibniz rule. Chen's identity says that the sum of these products is . Thus every term of degree at most agrees with the corresponding term in . The remaining Taylor terms have control order at least ; bounded norms make the estimate uniform in . For smooth fields the remainder has order .
On a partition , construct the Euler sequence by . For grid endpoints let
It vanishes on adjacent points. Insert into the next local map. The Lipschitz and composition estimates give
The Davie lemma therefore yields a bound independent of on small control intervals. Together with the Euler expansion this gives . A finite subdivision into small-control intervals gives a uniform global bound as well.
Interpolate each Euler sequence linearly in time, and take partitions with mesh tending to zero. Continuity of the control makes these paths uniformly equicontinuous: the grid estimates apply to the interval enlarged to its neighboring grid endpoints, and the enlargement tends to zero uniformly with the mesh. The paths are uniformly bounded. The Arzelà-Ascoli theorem supplies a uniformly convergent subsequence with continuous limit .
For arbitrary in a small control interval, choose neighboring grid endpoints converging to . Continuity of , the coefficient functions and the control lets us pass to the limit in the grid error estimate. We obtain the defining expansion with , and . This constructs a pathlevel RDE solution. No continuity theorem for the RDE solution map was used to obtain existence.