Conformational Energies of n-Alkanes and the Random Configuration of Higher Homologs Including Polymethylene

Enhancing Reverse Intersystem Crossing in Triptycene-TADF Emitters: Theoretical Insights into Reorganization Energy and Heavy Atom Effects

Thermally activated delayed fluorescence (TADF) emitters based on the triptycene skeleton demonstrate exceptional performance, superior stability, and low efficiency roll-off. Understanding the interplay between the luminescent properties of triptycene-TADF molecules and their assembly environments, along with their excited-state characteristics, necessitates a comprehensive theoretical exploration. Herein, we predict the photophysical properties of triptycene-TADF molecules in a thin film environment using the quantum mechanics/molecular mechanics method and quantify their substantial dependency on the heavy atom effects and reorganization energies using the Marcus-Levich theory. Our calculated photophysical properties for two recently reported molecules closely align with experimental values. We design three novel triptycene-TADF molecules by incorporating chalcogen elements (O, S, and Se) to modify the acceptor units. These newly designed molecules exhibit reduced reorganization energies and enhanced reverse intersystem crossing (RISC) rates. The heavy atom effect amplifies spin-orbit coupling, thereby facilitating the RISC process, particularly at a remarkably high rate of ∼109 s-1.

Thermally Activated Delayed Fluorescence of a Pyromellitic Diimide Derivative in the Film Environment Investigated by Combined QM/MM and MS-CASPT2 Methods

Arylene diimide compounds exhibit thermally activated delayed fluorescence (TADF), but its mechanism remains elusive. Herein we studied the TADF mechanism of a carbazole-substituted pyromellitic diimide derivative (CzPhPmDI) in poly(methyl methacrylate) (PMMA) film by using DFT, TD-DFT, and MS-CASPT2 methods within the QM/MM framework. We found that the TADF mechanism involves three electronic states (i.e., S₀, S₁, and T₁), but the T₂ state is not involved because its energy is higher than the S₁ state by 6.9 kcal/mol. By contrast, the T₁ state is only 3.2 kcal/mol lower than the S₁ state and such small energy difference benefits the reverse intersystem crossing (rISC) process from T₁ to S₁ thereto TADF. This point is seconded by relevant radiative and nonradiative rates calculated. At room temperature, the ISC rate from S₁ to T₁ is calculated to be 6.1 × 10⁶ s^-1, which is larger than the fluorescence emission rate, 2.2 × 10⁵ s^-1; thus, the dominant S₁ population converts to the T₁ state. However, in the T₁ state, the rISC process (1.8 × 10⁴ s^-1) becomes the most important channel because of the negligible phosphorescence emission rate (3.5 × 10^-2 s^-1). So, the T₁ population is still converted back to the S₁ state to fluoresce enabling TADF. Unfortunately, the rISC process is blocked in low temperature. Besides, we found that relevant Huang-Rhys factors have dominant contribution from low-frequency vibrational motion related to the torsional motion of functional groups. These gained insights could provide useful information for the design of organic TADF materials with excellent luminescence efficiency.

Conformationally Controlled Linear and Helical Hydrocarbons Bearing Extended Side Chains

Conformationally controlled flexible molecules are ideal for applications in medicine and materials, where shape matters but an ability to adapt to multiple and changing environments is often required. The conformation of flexible hydrocarbon chains bearing contiguous methyl substituents is controlled through the avoidance of syn-pentane interactions: alternating syn-anti isomers adopt a linear conformation while all-syn isomers adopt a helical conformation. From a simple diamond lattice analysis, larger substituents, which would be required for most potential applications, result in significant and unavoidable syn-pentane interactions, suggesting substantially reduced conformational control. Through a combination of computation, synthesis, and NMR analysis, we have identified a selection of substitution patterns that allow large groups to be incorporated on conformationally controlled linear and helical hydrocarbon chains. Surprisingly, when the methyl substituents of alternating syn-anti hydrocarbons are replaced with acetoxyethyl groups, the main chain of almost 95% of the population of molecules adopt a linear conformation. Here, the side chains adopt nonideal eclipsed conformations with the main chain, thus minimizing syn-pentane interactions. In the case of all-syn hydrocarbons, concurrent removal of some methyl groups on the main chain adjacent to the large substituents is required to maintain a high population of molecules adopting a helical conformation. This information can now be used to design flexible hydrocarbon chains displaying functional groups in a defined relative orientation for multivalent binding or cooperative reactivity, for example, in targeting the interfaces defined by disease-relevant protein-protein interactions.

Monte Carlo simulation of molecular and structural properties of random copolymer thin films

Molecular and structural properties of random copolymer thin films were studied by Monte Carlo simulation of coarse-grained copolymer model on the high coordination lattice. Random copolymer thin films with 50% comonomer fraction with varied interaction strength between comonomer units were studied. Intramolecular interaction was represented by Flory's rotational isomeric state (RIS) model of polyethylene (PE). The non-bonded interactions were treated by Lennard-Jones potential with different parameter sets for comonomer units. When the interaction strength of comonomer is increased, the density of thin films is higher in the inner region and significantly dropped near the free surface along the normal direction of the film. Higher bulk densities and narrower interfacial thicknesses are observed for random copolymer with stronger bead interaction. End beads were segregated at the free surface region and tend to orient perpendicularly to the surface. Copolymer size and shape were significantly changed as a function of bead interaction. For orientation of the whole chain, the largest molecular axis tends to orient along the film surface, but changed toward random orientation for weaker bead interaction. The intra-chain energies were decreased at the surface region while the non-bonded energies were increased. The stronger the comonomer bead attraction, the more change in energetics across the thin film.

Interfacial Modeling of Fibrinogen Adsorption onto LiNbO3 Single Crystal-Single Domain Surfaces

For the development of next-generation protein-based biosensor surfaces, it is important to understand how functional proteins, such as fibrinogen (FBG), interact with polar substrate surfaces in order to prepare highly sensitive points of medical care diagnostics. FBG, which is a fibrous protein with an extracellular matrix, has both positively and negatively charged regions on its 3-dimensional surface, which makes interpreting how it effectively binds to polarized surfaces challenging. In this study, single-crystal LiNbO₃ (LNO) substrates that have surface charges were used to investigate the adsorption of FBG protruding polar fragments on the positively and negatively charged LNO surfaces. We performed a combination of experiments and multi-scale molecular modeling to understand the binding of FBG in vacuum and water-solvated surfaces of LNO. XPS measurements showed that the FBG adsorption on LNO increased with increment in solution concentration on surfaces independent of charges. Multi-scale molecular modeling employing Quantum Mechanics, Monte Carlo, and Molecular Mechanics addressed the phenomenon of FBG fragment bonding on LNO surfaces. The binding simulation validated the experimental observation using zeta potential measurements which showed presence of solvated medium influenced the adsorption phenomenon due to the negative surface potential.

Role of Intrinsic Factors of Polyimides in Glass Transition Temperature: An Atomistic Investigation

Polyimides (PIs) are in high demand in the field of active matrix organic light-emitting diode displays because of their excellent heat resistance, chemical stability, and mechanical properties. However, the most critical key to their application is to further enhance their glass transition temperature (T_g), which directly affects the processing temperature of thin-film transistors on the PI films. Therefore, it is of great importance to study the factors that have an influence on the T_g of PIs. To accomplish this goal, PIs derived from pyromellitic acid dianhydride and three sets of isomeric imidazole-based diamines were investigated. The investigation, by computational methods, was to clarify the effect of intrinsic factors associated with the molecular structure of the PIs on their T_g and to construct a structure-T_g relationship for these PIs. For each model system, all-atom molecular dynamics simulations were used to identify and distinguish the effects of chain rigidity, fractional free volume (FFV), cohesive energy density, hydrogen-bonding interactions, and charge-transfer complex interactions on T_g. The results showed that the physical property, chain rigidity, has a direct impact on T_g regardless of the polymer backbone structure. A linear correlation between the increase of FFV and the decrease of T_g was not established due to the existence of hydrogen-bonding interactions, but the tendency was maintained. Furthermore, the formation of hydrogen bonds was found to have an indirect relationship with T_g. That is, the increase of intrachain hydrogen bonds would lead to a decrease in chain rigidity and consequently reduce the T_g value.

Structure-Property Relationships of Poly(ethylene carbonate) and Poly(propylene carbonate)

Conformational characteristics of poly(ethylene carbonate) (PEC) and poly(propylene carbonate) (PPC) have been revealed via molecular orbital (MO) calculations and nuclear magnetic resonance (NMR) experiments on model compounds with the same bond sequences as those of the polycarbonates. Bond conformations derived from the MO calculations on the models were in exact agreement with those from the NMR experiments. Both PEC and PPC were indicated to adopt distorted conformations including a number of gauche bonds and cover themselves with negative charges, thus failing to form a regular packing and remaining amorphous. The MO data were applied to the refined rotational isomeric state (RIS) calculations to yield configurational properties such as the characteristic ratio, its temperature coefficient, the configurational entropy, and average geometrical parameters of unperturbed PEC and PPC chains. In the RIS calculations on PPC, the regio- and stereosequences were generated according to the Bernoulli trial or Markov stochastic process. In consequence, it was shown that the configurational properties of PPC do not depend significantly on its regio- and stereoregularities. The internal energy contribution to rubberlike chain elasticity, calculated from the temperature coefficient of the characteristic ratio, has indicated the possibility that PEC and PPC will behave as elastomers. The practical applications and potential utilizations of the polycarbonates are discussed on the basis of the conformational characteristics and configurational properties.

Stretching and folding of 2-nanometer hydrocarbon rods

Linear alkanes CnH2n+2 in vacuum isolation are finite models for an infinite polyethylene chain. Using spontaneous Raman scattering in supersonic jet expansions for n = 13-21 in different spectral ranges, we determine the minimal chain length nh for the cohesion-driven folding of the preferred extended all-trans conformation into a hairpin structure. We treat fully stretched all-trans alkanes as molecular "nanorods" and derive Young's modulus E for the stretching of an isolated single-strand polyethylene fibre by extrapolating the longitudinal acoustic mode to infinite chain length. Two key quality parameters for accurate intra- and intermolecular force fields of hydrocarbons (nh = 18 ± 1, E = 305 ± 5 GPa) are thus derived with high accuracy from experimental spectroscopy.

Orientation and conformation of lipids in crystals of transmembrane proteins

Orientational order parameters and individual dihedral torsion angles are evaluated for phospholipid and glycolipid molecules that are resolved in X-ray structures of integral transmembrane proteins in crystals. The order parameters of the lipid chains and glycerol backbones in protein crystals are characterised by a much wider distribution of orientational order than is found in fluid lipid bilayers and reconstituted lipid-protein membranes. This indicates that the lipids that are resolved in crystals of membrane proteins are mostly not representative of the entire lipid-protein interface. Much of the chain configurational disorder of the membrane-bound lipids in crystals arises from C-C bonds in energetically disallowed skew conformations. This suggests configurational heterogeneity of the lipids at a single binding site: eclipsed conformations occur also in the glycerol backbone torsion angles and the C-C torsion angles of the lipid head groups. Conformations of the lipid glycerol backbone in protein crystals are not restricted to the gauche C1-C2 rotamers found invariably in phospholipid bilayer crystals. Lipid head-group conformations in the protein crystals also do not conform solely to the bent-down conformation, with gauche-gauche configuration of the phosphodiester, that is characteristic of phospholipid bilayer membranes. Stereochemical violations in the protein-bound lipids are evidenced by ester carboxyl groups in non-planar configurations, and even in the cis configuration. Some lipids have the incorrect enantiomeric configuration of the glycerol backbone, and many of the branched methyl groups in the phytanyl chains associated with bacteriorhodopsin have the incorrect S configuration.

Viscosity dependence of intramolecular excimer formation with 1,5-bis(1-pyrenylcarboxy)pentane in alkane solvents as a function of temperature

Intramolecular excimer formation with 1,5-bis(1-pyrenylcarboxy)pentane, (1PC(5)1PC) is studied as a function of temperature in a series of alkane solvents and in toluene, covering a wide range of solvent viscosities η, from 0.2 to 125 cP. The rate constant k(a) of the monomer → excimer reaction is determined from the effectively single exponential monomer fluorescence decays. For the viscosity dependence of k(a) in n-alkanes, the Stokes-Einstein relation k(a) ∼ η(-1.0) does not hold. Instead, k(a) is proportional to η(-α), with α increasing upon cooling, from 0.56 at 85 °C to 0.86 at -30 °C. The activation energy E(a) of excimer formation with 1PC(5)1PC, always larger than the activation energy E(T/η) of solvent viscous flow, decreases when the solvent viscosity becomes smaller, from 20.7 kJ/mol in n-hexadecane to 11.8 kJ/mol in n-butane, approaching a value of 11-12 kJ/mol for the low viscosity solvents. As the excimer formation process depends on the restricted diffusion of the 1PC end groups as well as on the C-O and C-C rotations in the -O(CH(2))(5)O- chain, the limiting barrier of 11-12 kJ/mol is attributed to the activation energy E(c) of the multiple bond rotations. This fractional viscosity dependence (α < 1.0) is caused by the multidimensional character of the barrier crossing in the excimer formation process. This multidimensional character should also be taken into account in investigations of polymers and biological media employing excimer formation.

Crucial roles of charged saccharide moieties in survival of gram negative bacteria against protamine revealed by combination of grazing incidence x-ray structural characterizations and Monte Carlo simulations

Grazing incidence x-ray scattering techniques and Monte Carlo (MC) simulations are combined to reveal the influence of molecular structure (genetic mutation) and divalent cations on the survival of gram negative bacteria against cationic peptides such as protamine. The former yields detailed structures of bacterial lipopolysaccharide (LPS) membranes with minimized radiation damages, while the minimal computer model based on the linearized Poisson-Boltzmann theory allows for the simulation of conformational changes of macromolecules (LPSs and peptides) that occur in the time scale of ms. The complementary combination of the structural characterizations and MC simulation demonstrates that the condensations of divalent ions (Ca2+ or Mg2+) in the negatively charged core saccharides are crucial for bacterial survival.

Photocurrent generation in nanostructured organic solar cells

Photocurrent generation in nanostructured organic solar cells is simulated using a dynamical Monte Carlo model that includes the generation and transport properties of both excitons and free charges. Incorporating both optical and electrical properties, we study the influence of the heterojunction nanostructure (e.g., planar vs bulk junctions) on donor-acceptor organic solar cell efficiencies based on the archetype materials copper phthalocyanine (CuPc) and C(60). Structures considered are planar and planar-mixed heterojunctions, homogeneous and phase-separated donor-acceptor (DA) mixtures, idealized structures composed of DA pillars, and nanocrystalline DA networks. The thickness dependence of absorption, exciton diffusion, and carrier collection efficiencies is studied for different morphologies, yielding results similar to those experimentally observed. The influences of charge mobility and exciton diffusion length are studied, and optimal device thicknesses are proposed for various structures. Simulations show that, with currently available materials, nanocrystalline network solar cells optimize both exciton diffusion and carrier collection, thus providing for highly efficient solar energy conversion. Estimations of achievable energy conversion efficiencies are made for the various nanostructures based on current simulations used in conjunction with experimentally obtained fill factors and open-circuit voltages for conventional small molecular weight materials combinations.

Min-map bias Monte Carlo for chain molecules: biased Monte Carlo sampling based on bijective minimum-to-minimum mapping

A novel Monte Carlo (MC) simulation scheme based on Theodorou's bijective mapping strategy [D. N. Theodorou, J. Chem. Phys. 124, 034109 (2006)] is introduced. This min-map bias Monte Carlo acts in combination with any other proper, bare MC. It carries over the bare MC move from the original configuration space Omega(0), where trial move acceptance may be low, to a different configuration space, Omega(1), where acceptance is higher. The bare MC move is then performed in Omega(1) and the resulting configuration is finally mapped back to Omega(0). Mappings between Omega(0) and Omega(1) entail weighted selection of trial configurations, the bias of which is subsequently removed in the overall acceptance criterion. The new method is applied, in conjunction with continuum configurational bias as bare MC scheme, to the simulation of explicit hydrogen linear alkanes in the canonical ensemble. Min-map bias MC is found to alleviate the pervasive problem of very low acceptance rates encountered when using an explicit molecular description.