Simulating the Lyotropic Phase Behavior of a Partially Self-Complementary DNA Tetramer.

DNA oligomers in solution have been found to develop liquid crystal phases via a hierarchical process that involves Watson-Crick base pairing, supramolecular assembly into columns of duplexes, and long-range ordering. The multiscale nature of this phenomenon makes it difficult to quantitatively describe and assess the importance of the various contributions, particularly for very short strands. We performed molecular dynamics simulations based on the coarse-grained oxDNA model, aiming to depict all of the assembly processes involved and the phase behavior of solutions of the DNA GCCG tetramers. We find good quantitative matching to experimental data at both levels of molecular association (thermal melting) and collective ordering (phase diagram). We characterize the isotropic state and the low-density nematic and high-density columnar liquid crystal phases in terms of molecular order, size of aggregates, and structure, together with their effects on diffusivity processes. We observe a cooperative aggregation mechanism in which the formation of dimers is less thermodynamically favored than the formation of longer aggregates.

Molecularly induced order promotes charge separation through delocalized charge-transfer states at donor-acceptor heterojunctions.

The energetic landscape at the interface between electron donating and accepting molecular materials favors efficient conversion of intermolecular charge-transfer (CT) states into free charge carriers (FCC) in high-performance organic solar cells. Here, we elucidate how interfacial energetics, charge generation and radiative recombination are affected by molecular arrangement. We experimentally determine the CT dissociation properties of a series of model, small molecule donor-acceptor blends, where the used acceptors (B2PYMPM, B3PYMPM and B4PYMPM) differ only in the nitrogen position of their lateral pyridine rings. We find that the formation of an ordered, face-on molecular packing in B4PYMPM is beneficial to efficient, field-independent charge separation, leading to fill factors above 70% in photovoltaic devices. This is rationalized by a comprehensive computational protocol showing that, compared to the more amorphous and isotropically oriented B2PYMPM, the higher structural order of B4PYMPM molecules leads to more delocalized CT states. Furthermore, we find no correlation between the quantum efficiency of FCC radiative recombination and the bound or unbound nature of the CT states. This work highlights the importance of structural ordering at donor-acceptor interfaces for efficient FCC generation and shows that less bound CT states do not preclude efficient radiative recombination.

Multidimensional Potential Energy Surfaces Resolved at the RASPT2 Level for Accurate Photoinduced Isomerization Dynamics of Azobenzene.

We have used state-of-the-art ab initio restricted active RASPT2 computations using a 16 orbitals, 18 electrons active space to produce an extended three-dimensional map of the potential energy surfaces (PESs) of the ground and first nπ* excited states of azobenzene along CNNC torsion and the two CNN bending angles, which are the most relevant coordinates for the trans-cis photoisomerization process. Through comparison with fully unconstrained optimizations performed at the same level of theory, we show that the three selected coordinates suffice to correctly describe the photoisomerization mechanism and the S1-S0 crossing seam. We also provide a map of the nonadiabatic coupling between the two states in the region where they get closer in energy. Eventually, we show that treating the two CNN bending angles as independent coordinates is fundamental to break the symmetry and couple the two electronic states. The accuracy of the S0 and S1 PESs and couplings was validated with semiclassical dynamics simulations in the reduced space of the scanned coordinates, showing results in good agreement with published full-coordinate dynamics.