The ability to prepare a specific superposition of electronic excited states leading to a transitory symmetry breaking of the electronic density in complex systems remains a challenging concern. We investigate how an initial coherence can be controlled by laser fields. The selected molecular system is a symmetric dimer of phenylene ethynylene presenting different interesting properties: Two bright nearly degenerate excited states coupled through a conical intersection are addressable by orthogonal transition dipole moments and are well separated from neighboring states. Creating a superposed state with equal weights corresponds to a right or left electronic localization as in the double-well system followed by a transitory oscillation between the two wells. To ensure a spectral bandwidth typically smaller than 0.25 eV, the pulse duration is in the tens of femtoseconds range, so nuclear motion cannot be neglected. Optimal control theory (OCT) is applied with guess fields that effectively create the target coherence in the absence of dephasing due to the vibrational baths. We analyze the field reshaping proposed by the control and we further fit a sequence of pulses on the optimal field. The overall result is efficient and robust disymmetry control over reasonable timescales of few tens of femtoseconds, exceeding the pulse duration. The monotonically convergent algorithm is combined with the hierarchical equations of motion (HEOM) able to treat strongly coupled non-Markovian dynamics. We also check the implementation of the combined OCT-HEOM approach in the tensor-train representation with propagation using the time-dependent variational method.