Structure, Energetics, and Dynamics of Pedal-Like Motion in Stilbene from Molecular Simulation and ab Initio Calculations

α-Methylstilbene Isomers: Relationship of Structure to Photophysics and Photochemistry

Significant differences in the photochemical and photophysical behavior of trans-α-methylstilbene and trans-stilbene have been attributed to structural changes caused by the steric requirements of the methyl group. We present here the X-ray structures of cis- and trans-α-methylstilbene (c- and t-MeSt). This is the first X-ray structure of a cis-stilbene. Despite the pronounced departure from phenyl group coplanarity, the solid-state packing of t-MeSt resembles that of trans-stilbene in that both exhibit disorder with a bicycle pedal structural relationship, dynamic in t-St but static in t-MeSt. We compare the X-ray structures with calculated structures. We also compare our steady state and transient photochemical and spectroscopic results with predictions in a recent theoretical paper that anticipated some of our experiments. Deviations from planarity imposed by the methyl substitution account for the shorter lifetimes of the trans excited states. The rapid torsional relaxation of ¹t-MeSt* to the twisted intermediate ¹p*, k_tp = 2.9 × 10¹² s^-1, observed using fs transient absorption spectroscopy, explains the sharp decrease in the fluorescence quantum yield of t-MeSt. We correct misconceptions that have appeared in the literature concerning the shape of the stilbene potential energy surface in S₁. The nonplanarity due to methyl substitution leads to chirality issues that are relevant in biological molecules such as the protonated Schiff bases of retinal in the opsins.

Investigations into conformational transitions and solvation structure of a 7-piperidino-5,9-methanobenzo[8] annulene in water

Solvation shell structure of a 7-piperidino-5,9-methanobenzo[8] annulene (PMA) in water has been investigated in ambient conditions using both molecular dynamics (MD) and Car-Parrinello molecular dynamics (CPMD) calculations. From the MD calculations, we find that this molecule exists in three major conformational states out of which two are in twist-boat forms and one in chair form. Due to the limited time scale accessible in CPMD simulations, we have studied all the three conformational states separately using CPMD. The molecular geometry, electronic charge distribution and solvation structure for all three forms are investigated. The stability order of the chair and twist-boat conformations in water solvent has been reversed when compared to the gaseous phase results and in the case of polar aprotic solvents (J. Org. Chem., 1999, 61, 5979). From the radial distribution function, we find that the solvent density around the chair form is significantly lower, which has to be directly related to the smaller solvent accessible area for this conformation and this is in complete agreement with earlier reports. Among the findings are that the solvation shell structure around the nitrogen atom in the chair form of PMA is considerably different from the open conformational forms or the twist-boat forms. The dipole moment for the closed form is found to be significantly larger when compared to the twist-boat forms.

Orientational Melting and Reorientational Motion in a Cubane Molecular Crystal: A Molecular Simulation Study

Detailed molecular simulations are carried out to investigate the effect of temperature on orientational order in cubane molecular crystal. We report a transition from an orientationally ordered to an orientationally disordered plastic crystalline phase in the temperature range 425-450 K. This is similar to the experimentally reported transition at 395 K. The nature of this transition is first order and is associated with a 4.8% increase in unit cell volume that is comparable to the experimentally reported unit cell volume change of 5.4% (Phys. Rev. Lett. 1997, 78, 4938). An orientational order parameter, eta(T), has been defined in terms of average angle of libration of a molecular 3-fold axis and the orientational melting has been characterized by using eta(T). The orientational melting is associated with an anomaly in specific heat at constant pressure (C(P)) and compressibility (kappa). The enthalpy of transition and entropy of transition associated with this orientational melting are 20.8 J mol(-1) and 0.046 J mol(-1) K(-1), respectively. The structure of crystalline as well as plastic crystalline phases is characterized by using various radial distribution functions and orientational distribution functions. The coefficient of thermal expansion of the plastic crystalline phase is more than twice that of the crystalline phase.

High-pressure study of adamantane: variable shape simulations up to 26 GPa

We report simulations of adamantane by carefully combining ab initio and empirical approaches to enable simulations with internal degrees of freedom on crystalline adamantane up to a pressure of 26 GPa. Two sets of simulations, assuming the adamantane molecule as a rigid (RB) and flexible body (FB), have been carried out as a function of pressure up to 26 GPa to understand changes in the crystal structure as well as molecular structure. The flexible body simulations have been performed by including 6 lowest frequency internal modes (obtained from DFT calculations performed with Gaussian98) out of the total of 72. The calculated variation in c/a and V/V(0) from the RB and FB calculations as a function of pressure have been compared with the experimental curve. Other relevant thermodynamic and structural properties reported are the radial distribution functions, deviation in the position of a given type of atom with respect to its position at standard pressure, delta(s), cell parameters, volume, and energy. With an increase in pressure, three additional peaks are seen to develop gradually at three different pressures in the center of mass (com)-com radial distribution function (rdf). We attribute these changes to structural transformations (probably second-order phase transitions) which is consistent with the three phase transitions reported by Vijayakumar et al. for adamantane in the pressure range of 1 atm-15 GPa. Our simulations also show that these additional peaks in the rdf's are associated with the differences between opposite and parallel spin neighbors of Greig and Pawley as well as the crystallographic directional dependence of intermolecular distances in the first three shells of the neighbors. Also, the structural quantities from the RB calculation show considerable deviation from the FB calculation for pressures greater than 5 GPa, which suggests that the rigid body assumption for molecules may not be valid above this pressure.

Temperature-induced dynamical conformational disorder in 4-vinyl benzoic acid molecular crystals: A molecular simulation study

Extensive molecular simulations are carried out as a function of temperature to understand and quantify the conformational disorder in molecular crystals of 4-vinyl benzoic acid. The conformational disorder is found to be dynamic and associated with a flip-flop motion of vinyl groups. The population of minor conformer is less than 3% up to 300 K and is 13.2% at 350 K and these results are consistent with the experimental observations. At still higher temperatures, the population of minor conformer increases up to 25%. The evolution of structure at both molecular and unit-cell level of the molecular crystal as a function of temperature has been characterized by various quantities such as radial distribution functions, average cell parameters, volume, and interaction energies. The van't Hoff plot shows a nonlinear behavior at lower temperatures as it has been reported recently by Ogawa and co-workers in the case of stilbene, azobenzene, and N-(4-methylbenzylidene)-4-methylaniline molecular crystals. A set of rigid body simulations were also carried out to quantify the effect of conformational disorder on structural quantities such as unit-cell volume and interaction energy. The anomalous shrinkage of vinyl C=C bond length as a function of temperature has been explained by combining the results of simulations and a set of constrained optimizations using ab initio electronic structure calculations for various molecular structures differing in torsional angle.

Effect of Pressure on Pedal Motion in Stilbene Molecular Crystals and Its Dependence on the Crystallographic Site

We report computer simulation of a stilbene molecular crystal as a function of pressure up to 4 GPa. Molecular structure and the crystal structure of stilbene have been characterized by calculating the radial distribution function and dihedral angle distribution, features associated with pedal motion and cell parameters. Results suggest that the population of minor conformer at site 2 disappears altogether above 1.25 GPa. In contrast, the population of minor conformer at site 1 remains at around 12%. Pedal motion is not observed beyond a pressure of 0.8 and 1.4 GPa at site 1 and site 2, respectively. Specific heat and compressibility exhibit an anomaly around 1.25 GPa. The anomaly seems to be associated with the disappearance of pedal motion at site 2. Initially, increase in pressure leads to an increase in the magnitude of lattice energy, but beyond 0.5 GPa it decreases.