The ODE system in Eq. 14 does not have a straightforward analytical solution, therefore we provide a numerical solution for the system of equations using the molecular and pulse parameters from Ref. 34: ε 0 = 2.20 × 10 9 Vm -1 and s(t) = sin 2 πt tp , with t p = 9.67 fs the duration of the pulse. The circular molecule is chosen to be Mg-Porphyrin, with µ C = -1.84 ea 0 and ℏω = 2.21 eV, and the linear molecule to be Cy3 (µ L = 5.04 ea 0 and ℏω = 2.27 eV). The resulting populations, |c i | 2 , and coherences as given by the imaginary components of c i c j , i.e., I[c i c j ], between the four states are shown in Fig. and, for θ = 0 (solid) and θ = π/2 (dashed). As our model indicates, the phase will almost exclusively affect the formation of the linear exciton state, |L⟩. This effect is especially apparent in the exciton coherences (Fig. ). Note that in our model the linear exciton state populates more than the circular due to the linear chromophore having a larger transition dipole.