Self-assembly and mesophase formation in a non-ionic chromonic liquid crystal: insights from bottom-up and top-down coarse-grained simulation models

Accessing the electronic structure of liquid crystalline semiconductors with bottom-up electronic coarse-graining

Understanding the relationship between multiscale morphology and electronic structure is a grand challenge for semiconducting soft materials. Computational studies aimed at characterizing these relationships require the complex integration of quantum-chemical (QC) calculations, all-atom and coarse-grained (CG) molecular dynamics simulations, and back-mapping approaches. However, these methods pose substantial computational challenges that limit their application to the requisite length scales of soft material morphologies. Here, we demonstrate the bottom-up electronic coarse-graining (ECG) of morphology-dependent electronic structure in the liquid-crystal-forming semiconductor, 2-(4-methoxyphenyl)-7-octyl-benzothienobenzothiophene (BTBT). ECG is applied to construct density functional theory (DFT)-accurate valence band Hamiltonians of the isotropic and smectic liquid crystal (LC) phases using only the CG representation of BTBT. By bypassing the atomistic resolution and its prohibitive computational costs, ECG enables the first calculations of the morphology dependence of the electronic structure of charge carriers across LC phases at the ∼20 nm length scale, with robust statistical sampling. Kinetic Monte Carlo (kMC) simulations reveal a strong morphology dependence on zero-field charge mobility among different LC phases as well as the presence of two-molecule charge carriers that act as traps and hinder charge transport. We leverage these results to further evaluate the feasibility of developing mesoscopic, field-based ECG models in future works. The fully CG approach to electronic property predictions in LC semiconductors opens a new computational direction for designing electronic processes in soft materials at their characteristic length scales.

Maximum Entropy Theory of Multiscale Coarse-Graining via Matching Thermodynamic Forces: Application to a Molecular Crystal (TATB)

The MSCG/FM (multiscale coarse-graining via force-matching) approach is an efficient supervised machine learning method to develop microscopically informed coarse-grained (CG) models. We present a theory based on the principle of maximum entropy (PME) enveloping the existing MSCG/FM approaches. This theory views the MSCG/FM method as a special case of matching the thermodynamic forces from the extended ensemble described by the set of thermodynamic (relevant) system coordinates. This set may include CG coordinates, the stress tensor, applied external fields, and so forth, and may be characterized by nonequilibrium conditions. Following the presentation of the theory, we discuss the consistent matching of both bonded and nonbonded interactions. The proposed PME formulation is used as a starting point to extend the MSCG/FM method to the constant strain ensemble, which together with the explicit matching of the bonded forces is better suited for coarse-graining anisotropic media at a submolecular resolution. The theory is demonstrated by performing the fine coarse-graining of crystalline 1,3,5-triamino-2,4,6-trinitrobenzene (TATB), a well-known insensitive molecular energetic material, which exhibits highly anisotropic mechanical properties.

Investigating anionic surfactant phase diagrams using dissipative particle dynamics: development of a transferable model

Dissipative particle dynamics (DPD) provides a powerful coarse-grained simulation technique for the study of a wide range of soft matter systems. Here, we investigate the transferability of DPD models to the prediction of anionic surfactant phase diagrams, taking advantage of fast parameter sweeps to optimise the choice of DPD parameters for these systems. Parameters are developed which provide a good representation of the phase diagrams of SDS (sodium dodecyl sulfate) and three different isomeric forms of LAS (linear alkylbenzene sulfonates) across an extensive concentration range. A high degree of transferability is seen, with parameters readily transferable to other systems, such as AES (alkyl ether sulfates). Excellent agreement is obtained with experimentally measured quantities, such as the lamellar layer spacing. Isosurfaces are produced from the surfactant head group, from which the second moment M of the isosurface normal distribution is calculated for different phase structures. Lyotropic liquid crystalline phases are characterised by a combination of the eigenvalues of M, radial distribution functions, and visual inspections.

Molecular Simulation Approaches to the Study of Thermotropic and Lyotropic Liquid Crystals

Over the last decade, the availability of computer time, together with new algorithms capable of exploiting parallel computer architectures, has opened up many possibilities in molecularly modelling liquid crystalline systems. This perspective article points to recent progress in modelling both thermotropic and lyotropic systems. For thermotropic nematics, the advent of improved molecular force fields can provide predictions for nematic clearing temperatures within a 10 K range. Such studies also provide valuable insights into the structure of more complex phases, where molecular organisation may be challenging to probe experimentally. Developments in coarse-grained models for thermotropics are discussed in the context of understanding the complex interplay of molecular packing, microphase separation and local interactions, and in developing methods for the calculation of material properties for thermotropics. We discuss progress towards the calculation of elastic constants, rotational viscosity coefficients, flexoelectric coefficients and helical twisting powers. The article also covers developments in modelling micelles, conventional lyotropic phases, lyotropic phase diagrams, and chromonic liquid crystals. For the latter, atomistic simulations have been particularly productive in clarifying the nature of the self-assembled aggregates in dilute solution. The development of effective coarse-grained models for chromonics is discussed in detail, including models that have demonstrated the formation of the chromonic N and M phases.

All-atom simulations of bent liquid crystal dimers: the twist-bend nematic phase and insights into conformational chirality

The liquid crystal dimer 1,7-bis-4-(4'-cyanobiphenyl)heptane (CB7CB) is known to exhibit a nematic-nematic phase transition, with the lower temperature phase identified as the twist-bend nematic (N_TB) phase. Despite the achiral nature of the mesogen, the N_TB phase demonstrates emergent chirality through the spontaneous formation of a helical structure. We present extensive molecular dynamics simulations of CB7CB using an all-atom force field. The N_TB phase is observed in this model and, upon heating, shows phase transitions into the nematic (N) and isotropic phases. The simulated N_TB phase returns a pitch of 8.35 nm and a conical tilt angle of 29°. Analysis of the bend angle between the mesogenic units reveals an average angle of 127°, which is invariant to the simulated phase. We have calculated distributions of the chirality order parameter, χ, for the ensemble of conformers in the N_TB and N phases. These distributions elucidate that CB7CB is statistically achiral but can adopt chiral conformers with no preference for a specific handedness. Furthermore, there is no change in the extent of conformational chirality between the N_TB and N phases. Using single-molecule stochastic dynamics simulations in the gas phase, we study the dimer series CBnCB (where n = 6, 7, 8 or 9) and CBX(CH₂)₅YCB (where X/Y = CH₂, O or S) in terms of the bend angle and conformational chirality. We confirm that the bent molecular shape determines the ability of a dimer to exhibit the N_TB phase rather than its potential to assume chiral conformers; as |χ|_max increases with the spacer length, but the even-membered dimers have a linear shape in contrast to the bent nature of dimers with spacers of odd parity. For CBX(CH₂)₅YCB, it is found that |χ|_max increases as the bend angle of the dimer decreases, while the flexibility of the dimers remains unchanged through the series.

Self-assembling, structure and nonlinear optical properties of fluorescent organic nanoparticles in water

Owing to their intense emission, low toxicity and solubility in aqueous medium, fluorescent organic nanoparticles (FONs) have emerged as promising alternatives to inorganic ones for the realization of exogenous probes for bioimaging applications. However, the intimate structure of FONs in solution, as well as the role played by intermolecular interactions on their optical properties, remains challenging to study. Following a recent Second-Harmonic Scattering (SHS) investigation led by two of us [Daniel et al., ACS Photonics, 2015, 2, 1209], we report herein a computational study of the structural organization and second-order nonlinear optical (NLO) properties of FONs based on dipolar chromophores incorporating a hydrophobic triphenylamine electron-donating unit and a slightly hydrophilic aldehyde electron-withdrawing unit at their extremities. Molecular dynamics simulations of the FON formation in water are associated with quantum chemical calculations, to provide insight into the molecular aggregation process, the molecular orientation of the dipolar dyes within the nanoparticles, and the dynamical behavior of their NLO properties. Moreover, the impact of intermolecular interactions on the NLO responses of the FONs is investigated by employing the tight-binding version of the recently developed simplified time-dependent density functional theory (sTD-DFT) approach, allowing the all-atom quantum mechanics treatment of nanoparticles.

The Martini Model in Materials Science

The Martini model, a coarse-grained force field initially developed with biomolecular simulations in mind, has found an increasing number of applications in the field of soft materials science. The model's underlying building block principle does not pose restrictions on its application beyond biomolecular systems. Here, the main applications to date of the Martini model in materials science are highlighted, and a perspective for the future developments in this field is given, particularly in light of recent developments such as the new version of the model, Martini 3.

"Dividing and Conquering" and "Caching" in Molecular Modeling

Molecular modeling is widely utilized in subjects including but not limited to physics, chemistry, biology, materials science and engineering. Impressive progress has been made in development of theories, algorithms and software packages. To divide and conquer, and to cache intermediate results have been long standing principles in development of algorithms. Not surprisingly, most important methodological advancements in more than half century of molecular modeling are various implementations of these two fundamental principles. In the mainstream classical computational molecular science, tremendous efforts have been invested on two lines of algorithm development. The first is coarse graining, which is to represent multiple basic particles in higher resolution modeling as a single larger and softer particle in lower resolution counterpart, with resulting force fields of partial transferability at the expense of some information loss. The second is enhanced sampling, which realizes "dividing and conquering" and/or "caching" in configurational space with focus either on reaction coordinates and collective variables as in metadynamics and related algorithms, or on the transition matrix and state discretization as in Markov state models. For this line of algorithms, spatial resolution is maintained but results are not transferable. Deep learning has been utilized to realize more efficient and accurate ways of "dividing and conquering" and "caching" along these two lines of algorithmic research. We proposed and demonstrated the local free energy landscape approach, a new framework for classical computational molecular science. This framework is based on a third class of algorithm that facilitates molecular modeling through partially transferable in resolution "caching" of distributions for local clusters of molecular degrees of freedom. Differences, connections and potential interactions among these three algorithmic directions are discussed, with the hope to stimulate development of more elegant, efficient and reliable formulations and algorithms for "dividing and conquering" and "caching" in complex molecular systems.

Atomistic simulation studies of ionic cyanine dyes: self-assembly and aggregate formation in aqueous solution

Cyanine dyes are known to form large-scale aggregates of various morphologies via spontaneous self-assembly in aqueous solution, akin to chromonic liquid crystals. Atomistic molecular dynamics simulations have been performed on four cyanine dyes: pseudoisocyanine chloride (PIC), pinacyanol chloride (PCYN), 5,5',6,6'-tetrachloro-1,1',3,3'-tetraethylbenzimidazolylcarbocyanine chloride (TTBC) and 1,1'-disulfopropyl-3,3'-diethyl-5,5',6,6'-tetrachloro-benzimidazolylcarbocyanine sodium salt (BIC). Simulations employed an optimised general AMBER force field and demonstrate the organisation of the dyes into stacked structures at dilute concentrations. The thermodynamics of self-assembly was studied by calculating potentials of mean force for n-mers (n = 2, 3 or 4), from which the free energies of association are determined. We report binding free energies in the range of 8 to 15k_BT for dimerisation, concordant with typical values for ionic chromonics (7 to 14k_BT), and examine the enthalpic and entropic contributions to the aggregation process. The self-assembly of these dyes yields two distinct classes of structures. We observe the formation of H-aggregate stacks for PCYN, with further complexity in these assemblies for PIC; where the aggregates contain shift and Y junction defects. TTBC and BIC associate into a J-aggregate sheet structure of unimolecular thickness, and is composed of a brickwork arrangement between molecules. These sheet structures are characteristic of the smectic chromonic mesophase, and such assemblies provide a route to the emergence of nanoscale tubular architectures.